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There are numerous physical situations in which a hole or leak is introduced in an otherwise closed 
chaotic system. The leak can have a natural origin, it can mimic measurement devices, and it 
can also be used to reveal dynamical properties of the closed system. In this paper we provide an 
unified treatment of leaking systems and we review applications to different physical problems, 
both in the classical and quantum pictures. Our treatment is based on the transient chaos theory 
of open systems, which is essential because real leaks have finite size and therefore estimations 
based on the closed system differ essentially from observations. The field of applications reviewed 
is very broad, ranging from planetary astronomy and hydrodynamical flows, to plasma physics 
and quantum fidelity. The theory is expanded and adapted to the case of partial leaks (partial 
absorption/transmission) with applications to room acoustics and optical microcavities in mind. 
Simulations in the limagon family of billiards illustrate the main text. Regarding billiard dynamics, 
we emphasize that a correct discrete time representation can only be given in terms of the so- 
called true-time maps, while traditional Poincare maps lead to erroneous results. We generalize 
Perron-Frobenius-type operators so that they describe true-time maps with partial leaks. 
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I. INTRODUCTION 

A. Motivation 

Perhaps the most important distinction in the tempo- 
ral evolution of a dynamical system is between persistent 
(asymptotic) and transient (finite-time) dynamics. Dy- 
namical systems theory reflects this division and has de- 
veloped specialized methods and tools to investigate per- 
sistent (e.g., strange attractors, asymptotic Lyapunov ex- 
ponents) and trans ient (e.g., chaotic saddle , escap e rates) 
chaotic dynamics ( |Lai and Tel 2011 Ott[ |1993[ ). These 
two approaches become connected when considering the 
effect of opening up a hole (or introducing a leak) in an 
otherwise closed chaotic system, converting by this per- 
sistent into transient chaos. Transient and persistent dy- 
namics appear in both conservative and dissipative sys- 
tems, and it is important to distinguish leakage (escape 
or removal of trajectories) from dissipation (contraction 
in the phase space). Introducing a leak never generates 



an extra phase space contraction and, e.g., a conserva- 
tive system remains conservative after becoming leaky. 
More than a tool to investigate the relationship between 
different theories, problems described by a closed chaotic 
system with a leak appear nowadays in a great variety of 
fields: 

• Room acoustics: the decay of the sound en- 
ergy characterized traditionally by the so-called 
reverberation time can be considered a conse- 
quence of leaks: opening s and absorbing surfaces 
on the room's boundary (Bauer and Bertsch||l990 



Legrand and Sornette[ |1990b ~ [l991a} |Mortessa gne 



et 



le 



1993). Absorbing surfaces provide examples 



of partial leaks. 

• Hydrodynamical flows and environmental 

sciences: the fact that certain regions of flows have 
special hydrodynamical features and might there- 
fore change the properties of particles advected into 
these regions can be described by the so-called re- 
setting mechanism ( |Pierrehumbert 1994), which is 
a kind of leak from the point of view of chaotic ad- 



vection ( Neufeld et al. , 


2000 Schneider et al. , 2005 


2007 


Schneider and r . 


?el| 2003| Schneider et al. 


2002 


Tuval et al. 


2004). 



• Planetary science and cosmology: the (inelas- 
tic) collision of a small body with larger planetary 
objects leads to a drastic change in its dynamics 
compared to that in a point mass approximation 
of the larger bodies. In a first approximation the 
problem can be treat ed in the point mas s picture 
as a loss due to leaks ( |Nagler , 2004[ 2005 ). Similar 
ideas apply in cosmology ( |Motter| " 2001 ) . 



• Optical microcavities: light rays in dielectric 
materials are partially transmitted and reflected 
(with the exception of regions where total inter- 
nal reflection takes place). Chaotic cavities can 
be constructed to provide a strong directionality of 
emission through such a partial leak, a requirement 



for laser ( 


Altmann , 


2009 Dettmann et al. 


2009 


Harayama and Shinoharal 201 1| Lee et al. 


2004 



Noc kel and Stone[ |1997| |Ryu et al.\ |2006| |Schwe- 



fel et a/.[|2004| |Shinohara et _ a/.[ |2010[ |2011[ |2 009; 



Wiersig and Hentschel| |2008| |Yan et al.[|2009| ) 



• Plasma physics: particles in magnetic confine- 
ment devices are lost through collisions with sen- 
sors, antennas, or the chambers wall itself. These 



regions play therefore the role of a leak (Evans 
etal.\pM2\ |Portela et al. | [20081 [20071 |Viana etaf. 
|2011|[Wingen et al.\\2067\. 



• Wave/quantum signatures of open systems: 

features related to that of a leaking classical dy- 
namics appear in properties such as the (fractal) 
distribution of eigens tates (|Casati et al\ [19 99a; 
Ermann et ~al\ [2^091 |Keathig et al.\ |2006l |Kuhl| 



et al.\ 2005| 


Nfonnenmacher and Schenck, 2008 Pe- 


drosa et al. 


2009 ) , the survival probability in sim- 




ulations anc 


experiments (Alt et al., 1996| 


1995 




Casati et al. 


, 1999b| |Fendrik and Wisniacki[ 


1997 


Friedman et al] |2001| jKaplan et al.\ 2001), and in 



the fractal Weyl's law (Ermann and Shepelyansky 
20101 |Kopp and Schomerus| |2010| |Lu et a/.||2003a 



Nonnenmacher[ 1201 1| |Ramilowski et al. H2009| |She- 



pelyansky 2008 



Wiersig and Main 2008). 



In dynamical systems theory, the idea of leaking an 
otherwise closed chaotic systems was first proposed by 
Pianigiani and Yorke as early as 1979: 

Picture an energy conserving billiard table 
with smooth obstacles so that all trajectories 
are unstable with respect to the initial data. 
Now suppose a small hole is cut in the table so 
that the ball can fall through. We would like 
to investigate the statistical behavior of such 
phenomena (Pianigiani and Yorke, 1979). 



Their main motivation was precisely to investigate tran- 
sient chaos, as opposed to persistent chaos. The leak- 
age procedure was therefore a tool to create transiently 
chaotic systems. Interestingly, the development of the 
theory of transient chaos happened not to follow this line 
over decades. 

The importance of this mathematical approach be- 
comes apparent when one realizes the multitude of situa- 
tions in which the leak region has a well-defined physical 
interpretation. The importance of this aspect was first 
emphasized by Smilansky and coworkers, who pointed 
out that any measurement (both classical and quantum) 
leads unavoidably to a kind of leaking of the system. 
They wrote in 1992: 

A discrete spectrum is a property of a closed 
system. However, the process of measuring 
the spectrum of a bounded system consists of 
coupling the system to an external continuum. 
Thus, for the purpose of measurement, the 
closed sys tem is turned into a scatt ering sys- 
tem (|Doron and Smilansky 1992a). 



Physical realizations of the leak can thus be either the 
effect of measurement devices or intrinsic properties of 
the system, such as absorbing materials in room acoustics 
or the collisions in planetary astronomy mentioned above. 

Apart from physical leaks, there are also different the- 
oretical motivations for considering leaking systems: 

• Leakage is a tool to understand the dynamics of 
closed system s, providing thus a sort of chaoti c 

1992a|b ). 



spectroscopy (Doron and Smilansky 



More generally, systems with leaks help monitoring 



2007 ) (see also ( Nag! 



or peeping at chaos (Bunimovich and Dettmann) 



2007)): Here, as 



er et al. , 

in ( Pianigiani and Yorke| 1979f ", billiards with leaks 
were the first systems investigated because they al- 
low a natural connection between the classical and 



3 



quantum pictu res (|Alt et a/.||1996||1995||Bauer an j 

Bertsch 



1990). 



Leaking systems have been explored in the con- 



text of synchronization of chaotic oscillators (Ja- 



cobs et al. 


1998 


), and of the control of chaos ( 


Bul- 


jan and Paar, 20011 Paar and Buljan[ 2000| 


bar 


and Pavin, 


1997). 



• Leakag e reveals the foliation s inside the closed sys- 
tem ( [Aguirre ancT Sanjuai ij |2003| |Aguirre et al. 
2009 Sanjuan et al] 2003 Schneider et al. 



2002 



that lead, e.g., to fractal exit boundaries (Bleher 
eTal] [I9§8l |Portela et al\ [20071 |Ree and Relchlf 



2002) 



• The distribution of Poincare recurrences, which is 
commonly used to quan tify properties of chaotic 
Hamiltonian dynamics (Chirikov and Shepelyan- 



sky, 1984 Zaslavsky 2002), is equivalent to the 



survival proba bility in systems with leaks ( Altmann 
and Tel[[2008| . 



• Several quantifications of wave or quantum chaos, 
such as Loschmidt e cho ( |Gorin et al] 2006[ Jacquod 
and Petitjeanl 120091) or fidelity decay (IPeresl 119841), 



can be realized physically in configurations that 
are analogous to introducing a localized leak in a 
closed syste m (|Ares and Wisniacki[ |2009| |Goussev 
and Richterl |2007HGoussev et al] |2008| |Hohmann 
et a/.l 120081 Kober et al] 120111). 



The common feature in all applications and theoret- 
ical procedures listed above is that one has some free- 
dom when choosing the opening, i.e. the leak in a well- 



defined closed chaotic system ^Schneider et al. , 2002 ) . 



This should be contrasted to genuinely open systems in 
which the openness is intrinsic, and only slight paramet- 
ric changes are physically realistic, which typically do 
not allow to go to the closed system limit. Although 
both classes of systems are dynamically open, one aim 
of our paper is to emphasize the benefits of considering 
leaking systems, which are more precisely defined by two 
key elements: 

(i) the existence of a well-defined closed system which 
can be used as a comparison; 

(ii) the possibility of controlling (some) properties of 
the leak such as position, size, shape, or reflectivity. 

Property (i) guarantees that one can compare transient 
and asymptotic dynamics and can be considered as a 
particular case of (ii) if the possibility of reducing the 
leak size to zero is assured. Leaking systems can be both 
dissipative and conservative (Hamiltonian). Within this 
latter category, we consider the problem of chaotic scat- 
tering (as typically defined, e.g., in ( jGaspard , 1998| )) to 



be beyond the scope of this review because it lacks the 
properties (i) and (ii) above 1 . 

Our main approach in this Review Article is based on 
transient chaos theory, which is described for the case of 
leaky systems and connected to different recent applica- 
tions. Our aim is to be understandable by nonspecialists 
interested in learning what the implications of nonlinear- 
dynamics theories are to specific applications. At the 
same time, we will emphasize how specific applications 
pose new questions to the theory. Thus, we devote spe- 
cial attention to developing a theory consistent with the 
following two aspects required by different applications: 

• leaks are not necessarily full holes, they might be 
"semi-permeable" , i.e, the energy content of trajec- 
tories entering a partial leak is partially transmitted 
and partially reflected. 

• discrete-time maps of open flows might lead to a 
loss of information over the temporal properties, 
and therefore it is essential to use the generalized 
concept of true-time maps |Kaufmann and Lust- 
feldl f200ll), which will be denned in Sec.II.CI). 



We note here that even though our focus and numerical 
illustrations are on billiards (Hamiltonian systems), the 
theoretical framework, and many of the specific results 
can be naturally extended to systems with dissipation. 

In the remaining of this section we motivate the gen- 
eral problem through a historical example and a simple 
simulation. In Sec. [Ii] we confront the simplest theory, 
based on the properties of the closed system, with the 
appropriate transient chaos theory for open systems. A 
generalization of this theory to partial leaks is also given. 



Section III is devoted to a Perron-Frobenius-type opera- 
tor formalism that is able to describe any kind of leaking 
dynamics. The main implications of transient chaos the- 



ory are explored in Sec. IV including the case of multiple 
leaks and emission. In Sec. Owe discuss how to describe 
the generic situation of weakly chaotic Hamiltonian sys- 
tems (mixed phase space). Finally, in Sec. 



VI 



we use our 



results to give a detailed view on some of the problems 
we started this Section with. Our conclusions appear in 



Sec. VII In the Appendices (Sec. VIII) we discuss some 
important but technical aspects of open billiards (like 
e.g., different types of measures and algorithms). 



B. Classical leaking: kinetic theory and Sabine's law 

Historically, perhaps the first problem involving sys- 
tems with leaks was one related to the kinetic theory of 
gases. Consider a container filled with ideal gas. How is 



1 In some scattering cases it is possible to "close" the inside of 
the scattering region (e.g., in the three disk problem, when the 
disks touch). However, in this case the incoming trajectories are 
unable to enter the interesting, chaotic scattering region. 



4 



the container emptied after a small leak / is introduced 
on its boundary? 

The answer follows from an elementary application of 
the kinetic theory. Here we follow basically the treatment 
of ( |Bauer and Bertsch , 1990) and ( |Joyce| [l975| ). Let I 
be a disk of area AA on the surface of the container and 
/(v,t) be the phase-space density of the particles, for 
which 



/ 



f(v,t)d 3 



N(t) 
V ' 



(1) 



where N(t) is the number of particles in the container 
of volume V at time t. The number of particles with 
velocity v leaving the system over a short time interval dt 
is then dN = dtAAvnf(v,t)d 3 v, where n is the normal 
vector of the surface at the leak /. The total number 
is then obtained by carrying out an integration over all 
velocities. Thus, the time derivative of the number N(t) 
of particles inside the container is 



dN(t) 
dt 



-AA J vn/(v,t)d 3 v, 



(2) 



where the minus sign indicates that particles are escap- 
ing. 

From molecular chaos, a basic ingredient of kinetic 
theory, it follows that an equilibrium phase-space den- 
sity is homogeneous (location independent) and isotropic: 
all velocity directions are equally probable. In the 
limit of small AA we can expect that there is a quasi- 
equilibrium distribution in the open system which sets in 
on a time scale shorter than the average lifetime. This 
quasi-equilibrium distribution shares the properties of 
that of closed systems. In this case, isotropy guaran- 
tees that the phase space density depends only on the 
modulus v of the velocity, and it is, therefore, convenient 
to use spherical coordinates for the integration. With 6 
being the angle between velocity and the normal vector, 
v.n = vcosO, Eq. ([2| reads as 



dN(t) 
dt 



poo r^/Z p'Ztt 

-AA / vf(v,t)v 2 dv / cos sin OdO / d<j>. 
Jo Jo Jo 

(3) 

The spherical symmetry of the phase space density ap- 
plied to ([!]) leads to 



f(v, t)v 2 dvAiT ■ 



N(t) 
V ' 



(4) 



and implies that w(v) = An f(v,t)v 2 V /N(t) is the proba- 
bility density for the velocity modulus v in the gas. 

Substituting this into (|3|, the first integral is found to 
be proportional to the average (v) of the velocity modu- 
lus. By carrying out all integrals, we find 



dN(t) AA(v) 



dt 



N(t). 



(5) 



As long as (v) is independent of time 2 the decay of 
the particle number is thus exponential, of the form of 
exp (— Kt), with an escape rate 



AA(v) 
W 



(6) 



For simplicity we will focus here on an ensemble of identi- 
cal particles with the same velocity v colliding elastically, 
in which case (v) \-> v in Eq. (J6| 3 . In this case, the re- 
ciprocal of the escape rate, which turns out to be the 
average lifetime, can be written as 



(r) 



1 



w 

AAv' 



(7) 



This is the time needed for the decay of the survivors 
by a factor of e. Since the result is linear in AA, and 
the velocity distribution is not only isotropic but also 
homogeneous, i.e. independent of the position along the 
wall, the expression remains valid for small leaks / of any 
shape, and AA is then the total leaking area. Since AA is 
small, (r) is large, and hence the assumption of a quasi- 
equilibrium distribution becomes justified a posteriori. 

An interesting, historically independent development 
is Sabine's law, a central object of architectural acous- 
tics. This law says that the residual sound intensity 



in a room decays exponentially with time (Joyce 



1975 



Mortessagne et aZ. , 1993). The duration to decay below 
the audible intensity is called the reverberation time, T r , 
and was found experimentally by W. C. Sabine in 1898 
to be 



T r = 6 In (10) 



W 
~AAc' 



(8) 



Here c is the sound velocity, and AA is the area of 
the union of all openings of the room (or of all energy 
absorbing surfaces after proper normalization). With 
c = 340 m/s, the numerical value of T r in SI units is 
T = 0.16V/ AA. Sabine's experiments also showed that 
the reverberation time for a pleasant sound perception is 
on the order of a few seconds for a good auditorium, and 
he designed concert halls (like, e.g., the Boston Music 
Hall) according to this principle. 

A comparison of ^ and (|8| reveals that Sabine's law 
is nothing but an application of the exponential decay of 
the particle number evaluated with v = c as the particle 
velocity. What is escaping in this problem is not parti- 
cles, but the energy of the sound waves. In the geomet- 
rical limit of room acoustics, one can consider the decay 



2 In a thermodynamical system the decay of particles eventually 
leads to a reduction in the pressure and temperature inside the 
container, and thus to a reduction of v. Here we shall be inter- 
ested in systems with constant v. The exponential decay is then 
valid for any t > 0. 

3 Notice that the elastic collisions of identical particles is equivalent 
to the dynamics of independent particles, what can be seen by 
exchanging particle labels at collision. 



5 



of energy as the problem of particles which travel along 
sound rays and lose part of their energy upon hitting 
the leak or the absorbing surface. The most remarkable 
property of expression (|8| is its universality: the rever- 
beration time is independent of the location of the sound 
source and of the shape of the room, provided the absorp- 
tion is weak and sound disperses uniformly around the 
room e.g., due to roughness or irregular geometry of the 
walls ( Mortessagne et al. 1993). The pre- factor 6 In (10) 
in Eq. ([8| results from the fact that in the acoustic con- 
text the decay below the audible intensity implies 60 dB, 
i.e., a decay factor of 10 6 , instead of a factor of e in 
Sabine's law (J8|, dated back to 1898, appears thus to be 
the first application of leaking chaotic dynamical systems 
in the history of science! 

Let us now take a closer look at the assumptions in the 
derivations above from the perspective of the dynamics. 
In terms of the modern theory of dynamical systems, the 
isotropy and homogeneity of the velocity distribution are 
a consequence of the following two hypotheses: 

HI: the leak size is small, so that the phase space dis- 
tribution does not change due to the openness; and 

H2: the particle dynamics inside the room is chaotic, 
more technically, the dynamics is ergodic and 
strongly mixing (implying exponential decay of cor- 
relations in time). 

Under these assumptions, the exponential decay is 
valid also in other dimensions. For instance, the escape 
rate in two-dimensional billiards is then found to be 



AAv 



(9) 



where AA stands for the length of the leak along the 
perimeter, and V is the two-dimensional volume, the 
area, of the billiard table. The replacement of factor 
4 by 7r is due to the geometrical change from spherical to 
planar polar coordinates. 

It should be noted that in both cases the survival prob- 
ability P(t) up to time t is 



p(t) 



(10) 



as obtained from Eq. ([5]), with initial condition P(0) = 
1. The probability p(r) to leave around the escape time 
r = t is the negative derivative of P(t) and thus 



p(r) 



(11) 



and P(t) = J t °° p(r)dr. Since the exponential decay 
holds from the very beginning, the average lifetime 



(r) 



POO POO 

/ t'p(t')dt' = / P{t')dt' (12) 
Jo Jo 



is found to be (r) = 1/ft, which was used in Eq. (|7|). The 
symbol (...) can be interpreted as an ensemble average. 



Finally, it is instructive to write both expressions ([6| 
and (Tol) of the escape rate as 



(^coll) 



(13) 



where = ^ is the relative size of the leak compared 
to the full wall surface, and can therefore be considered 
as the measure of the leak (taken with respect to the 
Liouville measure). The denominator has the dimension 
of time and is given by 



it 



coll/ 



or 



(^coll) 



ttV 
~Av J 



(14) 



in the three- and two-dimensional case, respectively. 
These (t co ii)s turn out to be the precise expressions of the 
average collision time between collisions with the walls 
(or, after a multiplication by v, the mean-free-path) well 
known for three- and two-dimensional clo sed rooms or 
billiards. A s emp hasized in Ref. ( Joyce| 1975} Mortes- 



sagne et al] |1993[ ) , these results were obtained already 
in the late XIX century by Czuber and Clausius. It is 
the average collision time that sets the characteristic time 
with which the average lifetime should be compared: for 
small leaks (r) <C (t C oii), i-e. the time scales strongly 
separate. 

By definition, (t co ii) can be expressed as the average 
over the local collision times £ co ii(x) for the phase space 
coordinates x along the wall as 



coil) = J £ C oii(x)d/i, (15) 



where \i is the uniform phase space (Lebesgue) measure 
characteristic of conservative systems. All the expres- 
sions found are taken with respect to the distributions 
characteristic of the unperturbed system. This is consis- 
tent with the small leak assumption (HI above) so that 
the escape rates obtained can be considered as a leading 
order result in a perturbation expansion where averages 
can yet be taken with respect to distributions character- 
izing the unperturbed (non-leaking) system. 

In the modern applications mentioned in Sec. |I.A[ how- 
ever, conditions HI (small leaks) and H2 (strong chaos) 
are typically not met. In this paper we will discuss in 
detail what happens in such cases. For instance, in any 
practical application the leak size is not, or cannot be 
made, infinitesimally small so that HI is violated and 
perturbation expansions break down. 

We shall see that an exponential decay of the survival 
probability typically remains valid for finite leak sizes, 
at least after some initial period. The estimation of the 
escape rate can be greatly improved by considering a sim- 
ilar expression as in Eq. (13), the measure of the leak di- 



vided by the average collision time, however, both taken 
with respect to a different measure: 



fl(I) fJL c (I), (tcoll) (*coll)c 



/ 



£coii(x)d/i c 



(16) 
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The new relevant measure fi c differs from the original 
Lebesgue measure /i since many particles have left the 
system by the time of observation, and what counts is 
the set of long-lived particles. With finite leaks, the decay 
differs substantially from the naive estimate obtained by 
using the original Liouville measure, as illustrated for 
our billiard example in Fig. [I] Even if precise definitions 
and further details will appear only later, the conceptual 
difference between \i and fi c is clear (compare Figs, 
and [6] below for an illustration of the dramatic changes 
in the phase space of the billiard). 

The theory of open dynamical systems tells us that 
this new c-measure is the s o-called conditionally invari- 
ant measure introduced in (Pianigiani and Yorke, 1979), 



which shall be didactically introduced and investigated 
in Sees. fljlV The violation of hypothesis H2 of strong 
chaos leads to even more radical changes, e.g. to a de- 
viation from the exponential decay for long times. This 
case will be investigated in Sec. [V| 



C. Billiard dynamics and true-time maps 

In dynamical systems theory, the kinetic problem with 
fixed velocities and Sabine's picture of room acoustics are 
described as billiard systems, as notice already by ( |Joyce| 
1975). Billiards are defined as bounded volumes or areas 
inside which particles move in a straight line with con- 
stant velocity v between collisions at the boundary, where 
they experience specular, elastic reflection (i.e., the abso- 
lute value of the velocity v is conserved and th e angle of 
incidence is equal to t he angle of reflection) ( [Chernov' 



and Markarian, 2006). For numerical and visualization 



convenience, we illustrate our results in two-dimensional 
billiards. In this case the dynamics can be described in a 
two-dimensional phase space, achieved by replacing the 
continuous time dynamics by a corresponding discrete 
time system f that maps the position s along the bound- 
ary and angle 6 of the n-th collision into those of the 
(n+ l)-th collision at the boundary. By convention, map 
/ connects the momenta right after the collisions. This 
procedure corresponds to a Poincare surface of section. 
The dimension of the full (four-dimensional) phase space 
is reduced by two (using momentum conservation and 
the condition of collision). The shape of the billiard's 
boundary uniquely defines the dynamics of the particles, 
and system-specific properties depend sensitively on this 
shape. It is convenient to write the phase space of the 
map in terms of Birkhoff coordinates 



sine 



in which case 



f : (x n ) (x n+ i) 



(17) 



(18) 



is area preserving ( Berry[[l981| Chernov and Markarian| 



2006). 



A faithful representation of the temporal dynamics of 
billiards requires augmenting ( 18 ) by keeping track of the 



information about the time of each trajectory: 

tn+l =tn+ £coll(x n+ l), (19) 

where t n denotes the time of the n-th collision at the 
boundary of the billiard, and t co \\ stands for the time 
between two subsequent collisions. In what follows we 
associate t co \\ with the Birkhoff coordinates of the later 
collision (x n+ i) in order to be able to speak about the 
collision times within the leak when systems with leaks 
are considered. 



Map (|18|)-(|19|) is called the true-time map as coined 
by (Kaufmann and Lustfeld, 2001), which is also fre- 



quently used in the billiard context [see, e.g., (Buni- 
movich and Dettmann, 2007[)]. 



More generally, true-time maps provide a link be- 
tween discrete-time maps and continuous-time flows in 
the same spirit as described by the mathematical con- 
cepts of suspended flows, special flows, or flows under a 
function (Gaspard 1998 Katok and Hasselblatt 1995). 



A true-time map is equivalent to the continuous-time rep- 
resentation, but leads to faster and more reliable results 
than a direct integration of the billiard flow. The dif- 
ferent collision times can be taken into account also in 
the Perron-Frobenius representation of the dynamics, as 
shown in Sec. IIII.DI 

In contrast, the often used Poincare map, represented 
by (18) alone, provides a distorted image of time. It im- 



plies associating to each pair of collision the same time 
interval and thus looses contact with the time dynamics 
of the continuous time physical system (e.g., it can over- 
estimate the importance of events with short collision 
times). The Poincare map generates a measure differ- 
ent from that of the true-time map, and leads thus to 
erroneous results. When speaking about maps in the bil- 
liard context, we shall, therefore, always mean true-time 
maps. Poincare maps of billiards will only be mentioned 
again in Table [TT] and Appendix |VIII.D[ to illustrate the 
difference to true-time maps. 



D. Example in a chaotic billiard 

The main properties of two-dimensional billiards can 
be illustrated by the family of limacon billiards intro- 
duced by (jRobnik 1983) whose borders are defined in 



polar coordinates (p, <j)) by limacon-like curves 
p(<t>) = S[l + ecos($], 



(20) 



where S scales the size and e controls the shape of the 
billiard. The ratio S/v defines the unit in which time t 
is measured, which is the only effect of S and v on the 
dynamics. We therefore set S = v = 1 in what follows, 
which implies that the perimeter length is A — 8, the 
billiard's area is V = 37r/2, and the mean collision time 
(14) is thus (t co ii) = 37r 2 /16. For convenience, we use 
the convention throughout the paper that the perimeter 
coordinate s is parameterized between —1 and +1 (see 
e.g. Figs. [2] and I]) . 
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FIG. 1 (Color online) Escape of particles in a billiard with a finite leak. The survival probability P(t) for the fully chaotic 
cardioid billiard [Eq. ( |2Q[ ) with e = 1] is shown for two different configurations of the leak / (same effective size but different 
leak positions xi, and reflectivity R). The first (bottom, red) curve corresponds to a full leak [see Eq. (21)] centered at the top 
of the billiard, si — 0.5, with size 2As — 5% of the perimeter as shown in the left inset. The second (next to the bottom, blue) 
corresponds to a partial leak with R = 0.5 [see Eq. (22)] centered at si = —0.25 with a size 2 As — 10% of the perimeter as 
shown in the right inset. The observed escape rates are k = 0.03002 ± 0.00007 (full leak) and k = 0.02904 ± 0.00003 (partial 
leak), clearly different from the predictions k — 0.0270 (full line), based on Sabine's law ( |13| )-( 14 ), and_ft* = 0.0277, Eq. (24) 
with = 5%. Initial conditions were uniformly distributed in the phase space (s,p — sin#, 



Fig.§. 



For s = we recover the circular billiard, exhibiting 
regular dynamics. For e = 1, Eq. (20) defines the car- 



dioid billiard, which is ergodic and strongly mixing (Rob- 



nik| |1983| |Wojtkowski| |1986[ ) , satisfying the hypothesis of 
strong chaos H2 (previous section). For < e < 1 the 
billiard typically shows coexistence of chaotic and regu- 
lar components in the phase space (|Dullin and Baecker| 



2001), and exhibits weak chaos. The collision time t co n 
needed for the true-time map ( 18 )-(fl9|) can be deter- 



mined numerically and is shown in Fig. [21 for the cardioid 
case (e = 1). 




FIG. 2 (Color online) Collision time t co ii(x) as a function of 
the phase-space coordinates x = (s,p) in the cardioid bil- 
liard, Eq. (20) with s = 1. Consistent with the convention 
in Eq. (19), t co ii(x) is defined as the distance (or time, since 
v = 1) between x and the previous collision / _1 (x). The dis- 
continuity close to the diagonal reflects the billiard 's cusp at 
s = zbl, see Fig. [3^l. 

We now introduce a leak in such a closed billiard and 
test the limitations of Sabine's prediction. For concrete- 
ness, consider removing 5% of the top part of the perime- 



ter of the fully chaotic Cardioid billiard, as shown in 
the lower left corner of Fig. [I] Numerical simulations 
of the survival probability of trajectories in this system 
yields an escape rate n = 0.030, which differs substan- 
tially from the escape rate n = 0.0270 obtained from 



Sabine's original estimate (13) by using /i(J) = 0.05 and 
(^coii) = 37r 2 /16 . In fact, Sabine's estimates holds for 



infinitesimally small leaks only, and a naive extension for 



finite leaks will be presented in Section II.A and leads to 
Eq. (24), which is a generalization of Sabine's prediction. 
This yields n* = 0.0277 which is still about 10% below 
the observed one. Although this difference appears to 
be small, it shows up in the exponent of an exponential 
time dependence. After 500 time units, the number of 



observed survivors is a factor of 



D 0. 003x500 



1.5 



4.5 



times smaller than the one based on the closed system 
estimate. In Fig. [I] it corresponds to the difference be- 
tween the dashed (generalized Sabine's formula) and the 
solid red (direct simulations) lines. 

This very basic observation is just the simplest tempo- 
ral manifestation of a series of discrepancies that will be 
discussed in this paper and that are all originated in the 
difference between the dynamics of the closed and of the 
leaking systems (see Figs. [3] and [4] for the illustration of 
the change in the phase space). All these illustrate the 
need for a deeper theoretical understanding of systems 
with leaks, beyond the results obtained under simplify- 
ing assumptions such as those used to obtain Sabine's 
law in Sec. II. Bl 

Before exposing the theory in the next section (Sec.[n|, 
we define in full generality the problem of introducing a 
leak in an otherwise closed system. We wish to emphasize 
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that our motivation for using two-dimensional billiards is 
visual convenience and direct connection to applications. 
The idea of introducing leaks in dynamical systems ap- 
plies to a much broader class of systems where the results 
of this paper can also be applied, like e.g., non-billiard 
type Hamiltonian systems, dissipative systems, and also 
in higher dimensions. 



E. Definition of the leak 

Consider a closed system described by a map f c i OS ed (x) . 
In this article we are mainly interested in maps f c iosed 
that admit chaotic motion, but the introduction of a leak 
is independent of this requirement. Choose the leak / as 
a subset of the phase space Q. In its simplest version, 
a particle is regarded as having escaped the system af- 
ter entering the region /. The dynamics can thus be 
described by the following map: 



but the intensity of each particle will change as 



Wi+i 



f(x„) 



^closed (Xn) 

escape 



if ^ ^ 7 
if x n G 7. 



(21) 



Since escape is considered to occur one step after enter- 
ing 7, map f is defined in 7. 

In the example shown in Fig. [I] the leak 7 is centered 
at the boundary point si = 0.5 with width 2As = 5%. In 
general, a leak 7 can be centered at any phase-space po- 
sition = (si,pi) G Q. The leak mentioned above corre- 
sponds thus to 7 = [si — As, si+ As] x [— 1, 1], representing 
a rectangular strip parallel to the p-axis. A prominent 
physical example of leaks represented by strips parallel 
to the s-axis is that of dielectric cavities. In this case 
light rays coming from a medium with higher refractive 
index (n- m > n ou t) are totally reflected if they collide with 

\P\ > ^critical = ^out/^in, where ^critical = sin(<9 crit ical) is 

the critical momentum (Critical is the critical angle). The 
leak is then \p\ < ^critical? s arbitrary. 

A general leak 7 can have arbitrary shape (e.g., cir- 
cular, square, oval, etc.) and can also be composed of 
disjoint regions: 7 = U7^. In this last case, a natural 
question is that of the nature of the set of initial con- 
ditions which lead to each 7^, i.e., of the properties of 
the escape basins Bi. This problem will be discussed in 
Sec. |IV.B| For presentational convenience we will focus 
on leaks at the billiard's boundary, in which case we can 
still faithfully represent the phase space with Birkhoff co- 
ordinates. For leaks inside the billiard a representation 
in the full phase space is needed then. 

There are also physically relevant types of leaks that go 
beyond the definition in Eq. (21). For instance, in room 
acoustics, or in the above mentioned dielectric cavities, it 
is very natural to consider objects with partial reflection 
and partial absorption (or transmission). In this case we 
associate to each particle an intensity J that monotoni- 
cally decays due to collisions at the leak. The dynamics 
of particles is given by the closed map x n+ i = /closed (x n ), 



J n+l 



[ Jn 



if x n 7 
if x n G 7, 



(22) 



where the reflection coefficient < R < 1 might depend 
also on the phase space position x within the leak. The 
full leak defined in Eq. (21 ) is recovered by taking R = 
(J n = J ) in Eq. (22). Altogether, a leak 7 is defined 



by its size, position, shape, and reflectivity. In Sec. |TV 
we show that all these different characteristics of the leak 
affect the observable quantities of interest. 



II. THEORY FOR FINITE LEAKS 

A. Theory based on closed system properties 

(a) 





FIG. 3 (Color online) Dynamics in the closed cardioid bil- 
liard [Eq. ( |2Q| with e = 1]. (a) Configuration space with 
parameterization of the perimeter s £ [— 1,+1] and colli- 
sion angle 0. (b) Phase space depicted in Birkhoff coordi- 
nates x = (s,p = sin(0)) obtained at the collisions with the 
boundary. Two trajectories are shown in (a) and (b), one long 
(gray lines / black dots) and one short (red). 



The spirit behind Sabine's theory described in Sec. LB 
is to calculate the observables quantities of the open sys- 
tem based on the properties of the closed system. While 
the results of this theory are exact only for infinitesi- 
mally small leaks, it is natural to extend them to systems 
with finite leaks. As already shown above, the dynamics 
of two-dimensional billiards can be conveniently repre- 
sented by the true-time map (fl8|)-([l9|). Since for limacon 



billiards both s and p change in [—1,1], f preserves the 
measure d\i = \ cos(0)d0ds. Figure pi illustrates how this 
map is applied in the case of the cardioid billiard. Upon 
the n-th collision with the wall the length s n along the 
perimeter is determined (measured from the point lying 
furthest from the cusp) along with p n = sin# n . Any tra- 
jectory in the configuration space (like the red and gray 
curves in Fig. |3^) are thus mapped on a sequence of 
points in discrete time (red and black dots) in Fig. [3]o, 
and the time is monitored via Eq. ( fl9| ) in the knowledge 
of t co ii(x) shown in Fig. [2] 

This billiard is fully chaotic and the measure \i is the 
Lebesgue measure. This means that the predictions of 
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the theory based on the closed system are extremely sim- 
ple: trajectories are assumed to follow the natural invari- 
ant density of the closed system, p M (x) = 1/4, i.e., they 
are uniformly distributed in x = (s,p). 

Let us apply this theory to estimate the escape rate 
of a chaotic system with a finite leak. The average col- 
lision time (tcoii) f° r the closed system [Eq. ( [l5| )] is in- 
dependent of the leak size, and Eq. ( [l4| ) remains valid. 
The escape rate estimation ft* again depends only on 
the size (measure) of the leak but this time we 

do not assume to be small. For instance, a leak 

I = [si — As, si + As] x [pi — Ap,pi + Ap] has size 2 As 
along the s axis, height 2Ap in area 4AsAp and a mea- 
sure jji(I) = AsAp. This is the measure of trajectories 
escaping on the time scale (t co ii) of one collision. The 
survival probability after n = t/{t co \\) collisions can be 
estimated as 



P(t) = (I -»(!)) 



t/ (tcoll) 



<L*t 



which yields a naive estimate for the escape rate 



n = 



■In(l-Mffl) 

(*coll) 



(23) 



(24) 



This can be considered a generalization of Sabine's law 
because it is a natural extension of Eq. ( 13 ) to finite 
Form ula ([24]) is usually attribu ted to Ey r ing's work in 
1930 flMortessagne et al\ |1993| ), but see QJoyce| |1975| > 
for a detailed historical account. The naive prediction 
in Fig. [I] was determined with ( 24 ) and still consider- 
ably differs from the measured decay. It is important 
to note that while Sabine's theory is exact for infinites- 
imally small leaks, Eq. (24) is just an approximation of 
the finite-size case. Although it leads to an improved 
understanding of the problem of room acoustics, it ne- 
glects the fact that the presence of a large leak essentially 
changes the dynamics because only a small portion of the 
closed system's orbits has sufficiently long lifetime to give 
a considerable contribution to both the escape rate and 
the average collision time. A precise understanding of 
the dynamics in systems with finite leaks, including an 
explanation of the behavior observed in Fig. [T] can only 
be given if one abandons the approach based on closed 
systems and adopt s a description in t erms of the theory 
of transient chaos (Lai and Tel, 2011). 



B. Theory based on transient chaos 

The basic idea of transient chaos theory is to look at 
the invariant set of orbits that never leave the system 
for both t ±oo. A key statement of the theory is 
that there is a non-attracting chaotic set in the phase 
space t hat is responsible for the transiently chaoti c dy - 
namics flGaspard| [1998] [Ottj |1993| |Tel and Gruiz[ [2006 ) . 
This set is a chaotic saddle and is of course drastically 
different from the chaotic set of the closed system. To 
illustrate this difference we present in Fig. [4] a leaky bil- 
liard, its chaotic saddle, and two short lived trajectory. 



It is apparent that the long-lived orbits are rather excep- 
tional and the saddle is very sparse: it is a measure zero 
object (with respect to Lebesgue measure), a set that ex- 
hibits double fractal character. The difference between 
the closed system's Sabine-type theories and the ones 
based on transient chaos can pictorially best be expressed 
by comparing Figs. [3}d and[4]3. It becomes evident that 
transient chaos is supported by a strongly selected and 
extremely ordered subset of the closed system's trajecto- 
ries. Hence the measures (/i and fi c ) with which averages 
should be taken are fundamentally different in the two 
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FIG. 4 (Color online) Dynamics in a leaky cardioid billiard 
[Eq. ( |20| ) with e = 1]. (a) Configuration space with the leak / 
centered around si = 0.5 with As = 0.1 (in the momentum 
space pi = and Ap =1). One short lived (red) and one long 
lived (blue) orbit are shown, (b) Phase space of the true-time 
map with the chaotic saddle (blue dots) and the short lived 
trajectory (red). 

The saddle is responsible for the exponential decay of 
the survival probability 



p( t ) 



(25) 



where ~ indicates an asymptotic equality in t. The es- 
cape rate k is a property of the saddle and is independent 
of the initial distribution of the trajectories used to rep- 
resent an ensemble. 

The invariant set of transient chaos is called a saddle 
because it possesses a stable and an unstable manifold. 
The stable (unstable) manifold is composed by all tra- 
jectories that approach the chaotic saddle for t — » oo 
in the direct (inverted) dynamics. These manifolds of 
repelling and attracting character are extremely impor- 
tant to understand the properties of the open system. 
This is the reason why the term chaotic saddle is more 
appropriate than the often used term repeller (see, e.g., 
Ref. ( |Gaspard]|1998| )), which suggests (erroneously) that 
only unstable directions exist. 

Here we present, in brief, the so-called sprinkler 
method, that can be used to calculate not only a chaotic 



saddle but also its stable and unstable manifolds ( Lai and 



Tel 2011 Tel and Gruiz, 2006). One starts with N > 1 



trajectories distributed uniformly over the phase space. 
One then chooses a time t* ^> 1/k and follows the time 
evolution of each initial point up to t*. Only trajectories 
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that do not escape are kept, whose number is approxi- 
mately Noe(~ Kt *\ If at* is sufficiently large (but not too 
large such that only a few points remain inside), trajecto- 
ries with this long lifetime come close to the saddle in the 
course of dynamical evolution, implying that their initial 
points must be in the immediate vicinity of the stable 
manifold of the saddle (or of the saddle itself), and their 
end points be close to the unstable manifold of the sad- 
dle. The latter is so because most points still inside after 
time t* are about to leave. The points from the middle 
of these trajectories (t « t*/2) are then in the vicinity of 
the saddle. In the spirit of true-time maps, we used a 
generalization of this method which considers an inter- 
val of different times [t*,t* + At*] (see Appendix [VIII. C ) 
to generate the chaotic saddle of Fig. and the corre- 
sponding manifolds plotted in Fig. [5^i,b. 
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has been demonstrated to b e relevant in several transient- 
chaos-related phenomena ( |Lai and Tel , 2011). Fig. [6] 
shows the c-measure over the full phase space of the car- 
dioid and indicates that the distribution is rather irregu- 
lar. This should be compared with the smooth Lebesgue 
measure characterizing the closed system. 




FIG. 6 (Color online) Density of trajectories on the unstable 
manifold shown in Fig. [5] This distribution corresponds to the 
c-measure, the measure according to which averages are to be 
taken in the transient chaos context. Note that the c-measure 
is defined within the leak. (jV = 10 s , t* = 80, At* = 80, see 
Appendix IVIII.C) 



FIG. 5 (Color online) (a) Stable and (b) unstable mani- 
folds of the cardioid billiard shown in Fig. ^] obtained by the 
sprinkl er meth od (JV = 10 s , t* = 120, At* = 40, see Ap- 
pendix IvnTc ) . 



From the construction above it is clear that the parti- 
cles being in the process of escape are distributed along 
the unstable manifold. When compensating the loss due 
to escape by pumping in new particles according to an 
appropriate way, which corresponds in practice to mul- 
tiplying the density by exp(^t), we obtain an invariant 
density as the one shown in Fig. [6j This stationary dis- 
tribution is known to be the conditionally invariant mea- 
sure. 

Traditionally, a measure /i c is said to be conditionally 
invariant (c-measure for short) if for any subset E of the 
region of interest Q (Demers and Young, 2006; Pianigiani 
and Yorke[ [1979] ) 



(26) 



This means that the c-measure is not directly invari- 
ant under the map ji c (f~ l (E)) ^ /j, c (E), but it is pre- 
served under the incorporation of the compensation fac- 
tor /i c (/ _1 (^)) < 1 (the c-measure of the set remaining 
in Q in one iteration). 



Although many c-meas u res might exist ( Collet and 
Maume-Deschamps] |2000 Demers and Young| |2QQ6^ 7 
here we are interested in the natural c-measure which 
is concentrated along the saddle's unstable manifold, and 



Dimensions of the invariant sets 

Both the chaotic saddle and its manifolds are frac- 
tal sets, as can clearly be seen from Figs. [4)o and Fig. 
[5^L,b. Commonly, there are (at least) two different di- 
mensions used to quantify the fractality of these sets, 
the box-counting dimensions, .Do, and the information 
dimensions D\. The former characterizes the mere geo- 
metrical pattern, the latter also the distribution of par- 
ticles on the pattern ( Ott , 1993| ). 

The chaotic saddle of a two-dimensional map has a 
clear direct product structure: it can locally be decom- 
posed into two Cantor-set-like components, one along 
each manifold. The dimensions along the unstable (sta- 
ble) manifold is called the partial dimension in the un- 
stable (stable) direction and is marked by an upper index 
1(2). None of the partial dimensions can be larger than 
one. The dimension Do and D\ of the chaotic saddle is 
the sum of the two correspond i ng pa rtial dimensions ( |Lai| 
and Tell 120111 |Tel and Gruizl 120061): 



D =D { 



(i) 



(2) 



D 1 = D{ 1] + D\ 



)(2) 



(27) 



The saddle might also contain very rarely visited, and 
thus atypical, regions. Consequently, the information 
dimension Di cannot be larger than the box-counting 
dimension Do? which naturally holds for the partial di- 
mensions, too: < Dq \ j = 1,2. The value of the 
box-counting dimension is found often to be close to that 
of the information dimension and it is then sufficient to 
use only one of them. 
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The manifolds are also of direct product structure, but 
one component of them is a line segment, an object of 
partial dimension 1 (see Fig. [5^i,b) . Since the saddle can 
be considered as the intersection of its own manifolds, 
the dimension Dq 1 ^ of the unstable (stable) manifold 
is 1 plus the partial dimension in the stable (unstable) 
direction 



(u) 
0,1 



1 + A 



(2) 
0,1' 



W 

0,1 



1 + D x 



(i) 
o,i- 



(28) 



It should be noted that the information dimensions 

of the unstable manifold is nothing but the information 

dimension of the c-measure that sits on this manifold. 

A simplifying feature occurs in Hamiltonian systems: 
due to time reversal symmetry the partial dimensions 
coincide 4 . Therefore in Hamiltonian cases, we have 



A),: 



2D t 



(i) 



0,1' 



(u) 



0,1 



to 



0,1 



1 + 1?, 



(1) 



0,1 



1 



A),iA 
(29) 

i.e., all relevant dimensions can be expressed by the par- 
tial dimensions in the unstable direction. 

It is an important result of transient chaos theory that 
this dimension can be expressed in a simple way by the 
escape rate and the average Lyapunov exponent A on the 
chaotic saddle: 



(i) 



1 



K 

A' 



(30) 



< 



This relation, the Kantz-Grassberger relation (Kantz and 
Grassberger| |1985| ) states that the dimension observed 
along the unstable direction deviates from 1 the more, the 
larger the ratio of the escape rate (a characteristics of the 
global instability of the saddle) to the average Lyapunov 
exponent (a characteristics of the local instability on the 
saddle) is 5 . 

Implications for systems with leaks 

All results described so far are valid for open systems in 
general, they are not particular to systems with leaks. In 
the latter case, however, the chaotic saddle and its man- 
ifolds depend sensitively on the leak /. This implies that 
the unstable manifold's dimension, and the c-measure in 
general, might strongly depend on the size, location, and 
shape of the leak. 

We assume throughout this work that (a) / is not too 
large so that trajectories do not trivially escape after a 
short time (P(t) ^ 0, for any t); and (b) ergodicity and 
the chaotic properties of the closed system lead to one 



This symmetry explains why Fig. |5^i,b are mirror images of each 
other with respect to the s axis (after the leak is removed from 
the unstable manifold). The closed cardioid has an additional 
symmetry (s, p) H> (— s, —p). 
5 Relation ( |30| is valid in dissipative cases as well, but the partial 
dimension along the stable manifold is then = D^X/\X'\, 

where Y is the average negative Lyapunov exponent on the sad- 
dle. 



and only one chaotic saddle after the leak is introduced. 
We also exclude the possibility of trajectories re-entering 
the billiard after hitting a leak. This assumption is nat- 
urally satisfied in convex billiard, but should be enforced 
in the limacon billiard for leaks around the cusp 6 . 

We can take advantage of the fact that a leaky sys- 
tem is obtained from a fully chaotic closed system. For 
instance, we can construct explicitly the chaotic saddle 
by extracting from the original phase space all the f c iosed 
images and pre-images of the set used as leak / in f ope n- 
To see this, in Fig. [7] leak / of Figs. [4] and [5] is shown to- 
gether with its forward and backward images. Compared 
to Fig. [5j we see that the white regions in the stable (un- 
stable) manifold correspond to the backward (forward) 
iterates of the leak. 

Due to area preservation and ergodicity of the closed 
system, the result of removing infinitely many images can 
only be a set of measure zero. The chaotic saddle in a 
leaking system is the set of points that remain in the 
complement of the leak / and all its images for both for- 
ward and backward iterations. Indeed, the complement 
of the union of the sets in panels (a) and (b) of Fig. [7] 
already provides a good approximation to Fig. [5|a) and 
(b), respectively, and the complement of both panels is a 
good approximation to Fig. Qb). The chaotic saddle is 
a fractal subset of the original chaotic set (the full phase 
space in our strongly chaotic example). Furthermore, all 
invariant sets of the leaky system (periodic orbits, man- 
ifolds of the saddle, etc.) are subsets of those in the 
corresponding closed system. 




FIG. 7 (Color online) Phase space (s,p) of the cardioid bil- 
liard showing the (closed system) iterates of the black vertical 
stripe that corresponds to I in Figs. [4] and [5]). (a) Backward 
iterates, (b) Forward iterates. In both cases the first (blue), 
second (magenta), and third (yellow) iterates are shown. 



An improved escape rate formula 



We now notice that since escape in Eq. (21) is consid- 



ered to occur one step after entering /, the map f with 
leak is defined in / and thus the unstable manifold of the 



Physically we can imagine that an absorbing material is placed 
on the border of the closed billiard so that trajectories crossing 
the leak are immediately absorbed. 
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chaotic saddle enters /. It is also possible to compute While ( Mortessagne et al. 1993) uses a Gaussian ap- 



measures of the leak, as indicated in Eq. (16). The re- proximation for the distribution of the collision times of 



gion of interest Q in Eq. ( 26 ) is identified with the closed 
system's phase space, /i c (^j = 1, and the compensation 
factor in this equation is obtained as 



long-lived trajectories, (Ryu et al. 2006) takes into ac- 



Mc (f- 1 (^)) = Mc (^\/) = i- Mc (/). 



(31) 



We now recall that the c-measure is invariant and is 
distributed along the unstable manifold. The escape rate 
can be estimated by the same simple calculations that 
lead to Eq. (24) by replacing the natural measure /x(J) 
by the c-measure fi c (I) of the leak. Carrying out also the 
averaging of the collision times (t co ii) with respect to the 
c-measure [see Eq. (fl6|] we find that 



count only the collision times inside the leak. None of 
them is equivalent to ( 32 ) or to the exact expressions in 
Sec. UlLCl 

Another general statement we can make about systems 
with leaks is that when the size of the leak goes to zero, 
the properties of the open system tend to those of the 
closed system (Aguirre and Sanjuan 2003; de Moura and 



k(I) 



ln(l- Mc (J)) 

(^coll)c 



(32) 



Letelier| 1999[ ), i.e., the theory of Sec. II. A becomes cor- 
rect. In particular fJL c (I) — » —> 0, which implies 
that k in Eq. ([32]) tends to k* in Eq. ([24]) and both tend 
to Sabine's prediction = /i(I)/(t co n), Eq. (13). In 
terms of dimensions, Do,i —> 2 but a nontrivial closed- 
system approximation of the information dimension 
can be be obtained from the Kantz-Grassberger rela- 
tion ([30l) with = k*, X(I) = A (Lyapu nov exponent 



where subscript c stands for the c-measure of the true- 
time map and ~ indicates approximate equality. The 
validity and limitations of Eq. ( 32 ) will be carefully dis- 
cussed in Sec. |III.C"l but by now it is instructive to discuss 
the implications of this equation. It clearly shows that for 
finite leaks the escape rate cannot be obtained from the 
properties of the closed system and the c-measure should 
be used, a measure which deviates essentially from that 
of the closed system. The difference between /x(J) and 
Hc{I) is visually clear from Fig. ^jp. 

Let us make a short remark on usual open maps, in 
contrast to true-time maps. Their escape rate we denote 
by 7 (to sharply distinguish from the continuous time 
or true-time map escape rate k) implying that the dis- 
crete time survival probability P(n) decays as e~ in . The 
escape rate formula 

e- 7(/) = 1 - /x c (J) 7 (J) = - ln(l - Mc(/)) (33) 



of the closed system) ( Neufeld et al. , 2000 ) . 



Periodic orbits in maps with leaks 

We also review the dynamics in general maps with 



has been known since QPianigiani and Yorke] |1979|) and of the closed system into two sets: Tn , the orbits 



leaks. Periodic orbits analysis ( |Cvitanovic et al , 2004) 
is a powerful method to investigate chaotic systems and 
also illustrates very clearly the spirit of leaking systems. 
Generally, a dense set of unstable periodic orbits is em- 
bedded in the chaotic saddle and this set can be used to 
obtain an expression f or the escape rate of the saddle. 
As pointed out in Ref. ( |Altmann and Tel , 2009), in sys- 
tems with leaks this can be done using either the periodic 
orbits of the open system (that never hit the leak) or us- 
ing only the periodic orbits of the closed system that hit 
the leak. To illustrate the simple arguments that lead 
to this, let us consider the case of computing the escape 
rate 7 of a generic discrete mapping f . First we split the 
set T^n 11 ^ of all periodic orbits of length n (i.e., all or- 
bits that have an integer period equal to n,n/2,n/3, ...) 

(inside) 



can also be obtained directly from Eq. Q3ip fo r leaky that have at least one point inside /, and the comple- 
mentary set rfa utslde ) ? 1 e ^ oil orbits for which all points 



maps (Altmann and Tel 2008 Paar and Bu 



Ljan, 



2000). 



It expresses the fact the c-measure /i c of the leak is the 
proportion of particles escaping via the leak within an 
iteration step. Since starting from the c-measure the 
decay is exponential from the very beginning, the pro- 
portion of the surviving particles after one time unit is 
exp (—7), of those who escape is 1 — exp (—7), and thus 
H c = 1 — exp (—7) which is equivalent to ( |33| ). 

Equation ( [32] ) leads to Eq. (33) when t co \\{x) = l.It 
shows also that when using the true time of the system to- 
gether with a surface of section at the billiard's boundary, 
it is essential to take into account that the average colli- 
sion time differs from (£ C oii) given in Eq. ( [14] ). As already 
anticipated in Eq. ( [l6| ) , with finite leaks the correct aver- 
age collision time is given by (t C oii)c = J ^coii( x )Pc( x )^ x ? 
where p c is the density of the c-measure characterizing 
the system in the presence of leak /. We note that dif- 
ferent corrections for n due to the collision times were 



are outside the leak /. In the limit of large n the fol- 



lowing relation holds for hyperbolic systems (Dorfman 
1999] [Ottl (19931) 



E 



1 



| A(r (outside) )r 



(34) 



where the sum is over all points i of periodic trajectories 
in r^ olltslde ) ^ anc [ \ is the largest eigenvalue of the n-fold 
iterated map f n on the orbit. 

Next, we notice that in the closed system 7 
escape). Therefore, 



(no 



1 



E 



1 



|A(r 



(all).. 



(35) 



suggested in (jMortessagne et a/.||1993||Ryu et a/.||2006[). where the sum is over all points i in ri^ ll \ Subtract 
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ing ( 34 ) from ( 35 ) we obtain 

1 - e -«7 = 



| A(r (inside) ) 



(36) 



where the sum now is over all points i in Fn nside \ i.e., all 
points that belong to periodic orbits that have at least 
one point in /. Altogether this means that even if we 
are allowed to probe the system only through the leak, 
the identification of the periodic orbits entering / suffices 
for the computation of the main properties of the chaotic 
saddle that exists inside the system. This can be applied 
also to more efficient methods based on expansions of the 
zeta function ( [Artuso et a/.|[l990| ). 

A simple relation can be obtained in uniformly expand- 
ing (piecewise linear) maps in which leaks are selected as 
elements of a Markov partition. In this case it is possible 
to prove that for two different leaks I\ and I2 the relation 
7(^1) > 7(^2) holds if and only if the sho rtest periodic 
orbit in Ii is shorter than the one in I2 ( |Bunimovich , 
20121). This follows also from Eq. (36) with constant A 



(as in piecewise linear maps with constant slope). 



C. Initial conditions and average escape times 

Typical observable quantities in transient chaos theory, 
such as k and the fractal dimensions, are independent of 
the choice of the density of initial conditions po(x) be- 
cause they are directly related to the properties of the 
invariant chaotic saddle. More precisely, the underlying 
assumption is that po( x ) overlaps the stable manifold of 
this saddle. The stable manifold of the chaotic saddle 
typically spreads through the phase space in a filamen- 
tary pattern (e.g., as in Fig. [5|i) and therefore smooth 
Po(x)'s will typically fulfill this requirement. Even in 
this typical case, there are important quantities that do 
depend on po( x ) such as any quantity averaged over a 
large number N of trajectories. The dependence on ini- 
tial conditions and the universal asymptotic decay of the 
survival probability P(t) ~ e~ Kt can be reconciled by 
noticing that for short times, t < t s , the escape of tra- 
jectories is non-universal, and P(t) ^ e~ Kt . Even if t s is 
short, a large fraction of the trajectories may escape for 
t < t s . 

Here we discuss in more detail the simplest yet repre- 
sentative case of the average lifetime (Altmann and Tel 
[2QQ9| 



(t) = Jim 



1 



N 



N ^ 

2 = 1 



rOC nC 

— / rp(r)dr = 
Jo Jo 



P{t')dt', 

(37) 

where T{ is the lifetime of trajectory i. Here we used 
p = -dP/dt [see Eq. Q], P(0) = 1, and P(t) 
faster than \jt. 

For maps, the averaged discrete lifetime is 



(i/) = £ n'p(n') 



(38) 



where p(n) is the probability to escape in the nth step. 
Notice that for true-time maps ((18), ( p~9]) ) in general 
( r ) 7^ M(*coii) CJ with (t C oii) c given by Eq. Q16|). Instead, 



(r) = (£t coll (x^)) 



where Vi is the total number of collisions along the i- 
th trajectory, x^' J ) is the position of the j-th collision 
(j = 1 . . . vi) of trajectory i that has initial condition 

x^'°), = ^ DjLi^coii(x (i ' j) ) is the mean collision 
time of trajectory z, and the average (...) is taken over 
i = 1 ... N trajectories. The reason for this difference 
is that for short times ensemble and temporal averages 
may differ considerably. 

We would like to see if (r) and (is) can be expressed 
as a function of other easily measurable quantities. We 
also try to find a relation between (r) and (v) for the 
following particular initial densities po{x): 

C-density: p c . We take initial conditions accord- 
ing to th e c-measure po( x ) = Pc( x )- As explained in 
Sec. II. B, p c (x) describes the escaping process and is 



achieved by re-scaling the surviving trajectories from an 
arbitrary smooth initial density. Therefore, for pp(x) = 
p c (x) we findp(t) = ne~ Kt for allt > and from Eqs. (37) 
and (|32|) the simple relation: 



(r)c = 



1 



(^coll)c 



ln(l - fi c (I)) 



(39) 



follows 

For maps with escape rate 7, the normalization of p(n) 
implies Y^=iP( n ) ~ 1? an d thus p(n) = (e 7 — l)e 



-jn 7 



This leads to a different result ( |Altmann and Tel 2009) 



1 - e-T 



Mc(/)' 



(40) 



In the last equality we have used (33). It is important to 



note that for maps obtained from flows, the c-densities p c 
of the map and flow (or true-time map) are usually dif- 
ferent due to the nontrivial collision time distribution. 

Rec urrences density: p r . As pointed out in 
Refs. QAltmann and Tel| |2008| |2009| ), there is an ini- 
tial density po(x) = /v(x) connected to the problem of 
Poincare recurrences that leads to simple results for (r) 
and (v) . The Poincare recurrence theorem asserts that in 
a closed dynamical system with an invariant measure /i, 
almost any trajectory (with respect to p) chosen inside a 
region / with p(I) > will return to I an infinite num- 
ber of times. The times TVs between two consecutive re- 
turns are called Poincare recurrence times, a central con- 



cept in dynamical systems theory (|Haydn et a/., 2005). 



The (cumulative) distribution of recurrence times P r (T) 



Note that e< - 1 w 7 for 7 < 1. 



14 



Large leaks 
Finite fjL(I) / fjL c (I) 



Limit of small leaks 



Po(x) 



c- measure: p c 



recurrence: p r 



natural, smooth: 



Pr,c 



Mean time 

continuous, t 

discrete, n 



(T)c = 



k - ii c (I)) 



(T)r = ^ = W^coll} + - 
fl(I) K 



1 - e-T 



Mm,. 



1 



<T> = 



Mc(J) 



(tcoll) 

f(-D 7 



Escape rate 

continuous, t 

discrete, n 



k „ ln(l - p c (J)) _^ ln(l-M/)) = K , 

(tcoll)c (tcoll) 

7 = - In (1 - Mc(/)) # - ln(l - = 7* 



(*coll) (' 

7 = mW = 



TABLE I Dependence of the average lifetime (r) and (z/) for flows and maps, respectively, on the initial distribution po(x), 
and expressions for the escape rate k, 7. The natural invariant measure of the leak (closed system) p (I) and the conditionally 
invariant measure /jl c (I) of the leak / coincide only in the limit of small leaks. Table inspired by Ref. (Altmann and Tel 2009) 



is also used to quantify chaotic properties of specific sys- 
tems ( Chirikov and Shepelyansky , |1984[ ) . Figure [8] illus- 
trates the Poincare recurrence set-up in billiard systems. 
Using the notation introduced in Fig. [8j the average re- 
currence time T of a single long trajectory is calculated 



as 



1 N 



^total _ jVcoll(^coll) _ (£coll) 



(41) 



where Ti is the zth recurrence time along the trajectory, 
(£ C oii) is the average collision time of a typical trajectory 
starting inside the leak that, due to ergodicity, coincides 
with the closed system (t co ii) given by Eq. (15), N co \\ is 
the total number of collisions up to time ttotab N is the 
number of such collisions inside /, and = N/N co \\ is 
the fraction of collisions on /. Relation (41) is valid for 
large t totah N, N co u. 
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FIG. 8 Schematic illustration of Poincare recurrences in a 
closed billiard. The symbols indicate the times t of the n-th 
collision with the boundary of a single trajectory. The recur- 
rence times TVs are defined as the times between successive 
collisions in /. In the total time ttotai there are N co n collisions, 
out of which N collisions are inside the region /. 



For maps, analogously, we obtain that the average dis- 
crete recurrence time AT is given by 



The inverse relation between measure 
time shown in Eqs. (41) and (|42* 



(42) 

and recurrence 
known as Kac's 



lemma ( |Kac[ |1959| |Zaslavsky[ [2002 ) . 



The connection to systems with leaks is achieved by 
identifying the recurrence region and the leak /. One 
can find an appropriate initial density po( x ) = Pr(x) f° r 
the open case for which the survival probability in the 
presence of leak / coincides with the recurrence time dis- 
tribution to / in the closed system 



Pr(T) = P(t) 



(43) 



for any t = T. This can be done by using the positions 
x G / of the N recurrent points as initial conditions (see 
Fig. |8|. Due to the ergodicity of the closed chaotic sys- 
tem, in the limit of long times, the points of this single 
trajectory are distributed according to the natural den- 
sity p M (x) of the closed system, justifying the notation 
/i(J) in the equation above. In the case of the billiard sys- 
tems discussed here, p r (x) corresponds to initial condi- 
tions in the leak, uniformly distributed in x = but 
with velocities pointing inward, i.e., not escaping through 
the leak. Physically, this corresponds to shooting trajec- 
tories into the billiard through the leak. 

If time is counted discretely, p r (x) corresponds to the 
next iteration of the uniform distribution (natural mea- 
sure of closed system p^(x)) of initial conditions in x G /. 
In the general case of an invertible map f this is obtained 
by applying the Perron- Frobenius (Dorfman, 1999 Gas- 
pard[ |1998[ ) operator as flAltmann and Tel[|2009| ): 



Pr(x) = 



j(f-i(x)niH/) 



for xE/(J), (44) 



15 



where f ~ 1 (x)f)I denotes the points that come from /, J is 
the Jacobian of the map, and the factor ensures nor- 
malization. Notice that the single iteration introduced in 
the definition of p r is compensated at the end because the 
escape is considered also one time step after entering /, 
see Eq. With p (x) = p r (x), P(t) = P r (T = t) for 

all t > 0, showing that the problem of Poincare recur- 
rence can be interpreted as a specific problem of systems 
with leaks. In particular, P r (T) ~ e~ KT with the escape 
rate k of the system opened up in /. The average lifetime 
is given by Eq. (41) as 



(r) 



(t C o\\) / (^coll)c 



(^coll)c 



/*(/) ^ fi c (I) r ln(l- Mc (/))' 



For maps, Eq. (42) implies that 
{v)=M = 



til)' 

Note that both relations above reveal that for p r there is 
no difference in the dependence on the measure for maps 
and flows, in contrast to Eqs. (39) and (40) obtained with 
p c as the initial density. 

Closed system density: or any smooth p s . 
The most popular initial conditions for systems with 
leaks is by far p M . This corresponds to introducing the 
leak after the trajectories have reached equilibrium in- 
side the system. While this is sometimes implicitly as- 
sumed as natural, the results of this section show that p^ 
is not the only possibility, and certainly not the sim- 
plest one. For p M and any other smooth initial distri- 
bution po = p si the surviving trajectories relax towards 
the c-density p c ( x ), the natural densit y of the open sys- 
tem. The numerical investigations of ( jAltmann and Tel] 
2009 ) indicate that results with p^ or other smooth den- 
sities p s are usually similar to those obtained by p c , and 
more different from those obtained by p r . This can be un- 
derstood by noting that the hyperbolicity of the chaotic 
saddle leads to a fast convergence po(x) —> p c (x). While 
a convergence also holds for p r , it is made much slower 
by p r 's lack of any overlap with the chaotic saddle, with 
orbits t aking, therefore, a long er time to approach the 
saddle flAltmann and Tel[ |2009[ ) . 

The results of this section are summarized in Tab. [H 
All formulae here were confirmed within a 3% precision 
for the cardioid billiard. In the limit of small leak, the 
simple case (r) = 1/ n discussed in Eq. (12), and Sabine's 



law in the form of (13) are recovered, independent of 
the initial distribution. With a finite leak we obtain the 
values presented in Tab. [TTJ 

D. Extension to partial leaks 

So far we have restricted our theory to the case of full 
leaks that completely transmit the trajectories falling on 
them. At the end of Sec. LE we have emphasized the 



need to extend the theory to leaks that partially trans- 
mit and partially reflect trajectories. This problem can 



Po(x) 


Pc Pr Pf^,s 


Mean time 

continuous, t 
discrete, n 


(r) c = 15.24 (r) r = 18.50 (r)^ s = 14.23 
(v) c = 8.28 (i/) r = 10.0 {v)^ s = 7.78 


Escape rate 

continuous, t 
discrete, n 


k = 0.06559 ^ k* = 0.05693 
lPma P = 0.1286 ^ 7* = 0.1054 



TABLE II Numerical results for the average lifetime in the 
cardioid bi llia r d w ith a leak si = 0.5, 2As = 10% (as in 
Figs. [4| [5] [42] |43| and [44]). Other data: (t co ii) = 1.85055, 
fi(I) = 0.1, /icTO = 0.1175, (tcoii)c = 1-916. In oder to 
illustrate the case of maps, instead of the true-time map, ex- 
clusively for this simulation we have used the Poincare map 
of the billiar d. In or der to minimize the effect of sliding or- 
bits (see Sec. VIII. D ) we used in all simulations the following 
restrictions: a cut-off in the maximum collision time at 83 
collisions (t = 158 in Fig. 43), and p s taken with s £ [—1, 1], 
p £ [—0.9, 0.9], and p c was built by iterating p s . For jpmap we 
use a more restrictive cut-off, 2(t) /XjS , because it is more sen- 
sitive to the sliding orbits. Th e quantities /jl c (I) and (£ co ii)c 
were calculated using Eq. (116) 



be handled by associating to each trajectory i an inten- 
sity that decreases with time J^+i — R(^)Jn^ clue to 
collisions at the leak regions with reflectivity i?(x) < 1, 
as described by (22). Here we show how the results of 
Sec. III. Bl can be extended to this case. 

For simplicity we consider that all trajectories of 
the initial ensemble start with intensity Jo = 1. Non- 
homogeneous energy intensities can be achieved by ma- 
nipulating the density of initial trajectories po( x )> as dis- 
cussed in the previous section. For longer times the en- 
ergy density depends not only on the density of trajec- 
tories p(x, t) but also on their intensities, which can be 
thought as weights < J < 1 attributed to each tra- 
jectory. The survival probability P(t) of trajectories is 
physically re-interpreted as the fraction P(t) of the total 
energy still inside the system. In fully chaotic system, 
P(t) is expected to decay asymptotically as 



p( t ) 



with k as the energy escape rate. In practice, P(t) can 
be obtained as a sum over all N trajectories 



Pit) « * £f =1 J t w 

= In Jt{x)p{x,t)dx = J n p(x,t)dx, 



(45) 



where p(x, t) is the density and Jt(x) is the intensity of 
trajectories in x at time t. The intensity of each tra- 
jectory i at time t is given by = J^IL^f R(x.^), 
with Vi(t) as the number of collisions of trajectory i until 
time t and x^ the position of these collisions. We have 
defined the combined density 



(46) 
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FIG. 9 (Color online) Invariant densities for the Cardioid billiard considered in Figs. [4] [5] and [6] with partial leak R = 0.1 
for x G /. The different panels show the densities associated to the (a) stable manifold, (b) chaotic saddle, and (c) unstable 
manifold. The densities were obtained considering the intensity weighted trajectories at times 1/k* <C 60 < t < 80 and plotting 
at (a) their initial conditions, (b) their positions around the time t/2, and (c) their positions over the time interval 60 < t < 80. 
See Appendix |VIII.C] for details. 



which is, by convention, not normalized for t > 0. The 
appearance of p in Eq. (45) reflects the more general fact 



that the division between trajectories and intensities has 
no observable consequence, despite their natural inter- 
pretation in physical terms and their natural implemen- 
tation in ray simulations. It is only their combination, 
Eq. (46), that leads to a physically relevant density for 
the case of partial leaks, both for the c-measure and for 
estimations based on p. A mathematical description co- 
herent with this interpretation considers operators act- 
ing on densities, and is given in the next section. In the 
remaining of this section we consider how partial leaks 
modify the closed system (Sec. |ILA) and transient chaos 



(Sec. II.B) theories discussed above. 



We start with a straightforward extension of the naive 
estimate of the escape rate. When estimating the trans- 
mission happening through the leak one has to account 
for the the partial transmission because for a trajectory 
in x G / only a fraction 1 — i?(x) of its intensity is lost. 
Therefore, the measure of the leak (p(I) = j T dp) used 
in Eq. (24) to compute the naive estimate n = ft* has to 



be replaced by 



m(I) = J (1 - i?(x))d/i = J (I- R(s, 0))^ cos OdQds, 

(47) 

where i?(x) is the (position dependent) reflection coeffi- 
cient. The argument leads to 

.. = -ln(l-m(J)) 

(tco\\) 



see also (Joyce 1975 Ryu et al. 2006). This shows that 
the closed-system estimation for a leak with R = 0.5 is 
the same as the one for a leak of half of the size, which 
was tested in Fig. [I] (and proved to be inaccurate). 

The essential extension of the transient chaos theory 
is to consider a (modified) c-density p c which should be 
normalized. As p c , p c can also be achieved by re-scaling. 



Here the intensity should be multiplied by a factor e Kt 
to compensate the global decay of P(t). This re-scaling 
implies that the c-measure p c (E) of a region E G ^ in 
a partially open system is the fraction of the intensity 
at time t — >• oo that is in E. A more proper estimate 
of the escape rate k from this measure has to include, 
as i n Eq. (|47| ), the partial transmission through the leak 



_m 

iann| 



as (Altmannl 2009) 



m c (J) = J (I - R(x))dp c = J (I - R(x))p c (x)dx. (48) 



The analogue of the escape rate formula ( 32 ) for a system 
with partial leaks is thus 



ln(l - m c (/)) 

(icoll)c 



(49) 



where c in (t C oii)c indicates average taken with respect to 

Pc- 

The complete extension of the transient chaos theory 



presented in Sec. II.B| is considerably more subtle and re- 
quires the extension of the invariant sets discussed above 
(chaotic saddle and its invariant manifolds). We illus- 
trate our general considerations introducing uniform re- 
flectivity R = 0.1 in the leak of the Cardioid billiard of 
Figs. [4 [5] The procedures used previously (see Figs. 4j 
|5| and 6j were employed to obtain Fig. [9j but now the 
intensity of the trajectories were used as weights at- 
tached to each trajectory, see Appendix. |VIII.B| for de- 
tails. The support of the regions with non- negligible den- 
sities shown in Fig. [9] corresponds to the stable manifold, 
chaotic saddle, and unstable manifold, respectively. In 
the case of partial leaks it is important to distinguish 
the manifold of the trajectories from the manifold of the 
(trajectories weighted with) intensities J. Whenever the 
reflection coefficient is non-vanishing, i?(x) 7^ for all 
xG trajectories survive forever (and correspond thus 
to that of the closed problem). The chaotic set is then 
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the full phase space, its manifolds are also space-filling, 
and are not very informative. The interesting patterns 
present in all panels of Fig. [9] motivate us to interpret the 
results in terms of manifolds of the (trajectories weighted 
with) intensities. This generalizes the case of full leak 
(# = 0forxel). 

Consider first the case of the stable manifold W s . 
While in the full leak case W f ^ n is obviously outside the 
leak (e.g., in Fig. in the partial leak case W s should 
certainly include trajectories which start in the leak and 
never return. These new trajectories are nothing but the 
closed map pre-image of W f ^ n in /. Indeed, in Fig. [9^a) 
the filamentary structure extends inside I (with much 
lower intensity). These new trajectories modify the man- 
ifolds in the whole phase space because their pre-images 
can be outside the leak. More generally, it is natural to 
consider Wf nll as the set of all points that carry a nonvan- 



ishing intensity J for arbitrarily large t (Altmann, 2009) 



III. OPERATOR FORMALISM 

In the previous sections we used extended trajectories 
which depend not only on the phase space coordinate x of 
the map but also on the true time t and on the intensity 
J. These two quantities are labels which are attached 
to each trajectory and change depending on (but do not 
affect) the sequence of collisions x^, i = 0. . .n. In this 
section we develop a more elegant formalism based on op- 
erators acting on densities p(x, t) that naturally accounts 
for both of these aspects, allowing for a more rigorous 
treatment. In the spirit of systems with leaks, and simi- 
larly to Sec. [TT[ we start with the case of closed systems 
that will be used for comparison with the leaking case. 



f (s,p) 




x G W s <^> for t -> +oc, J t (x) -> J +00 (x) > 0, (50) 



where Joo(x) is the asymptotic intensity. The same rea- 
soning applies to the unstable manifold W u , in which 
case the limit t —> — oo is taken in Eq. ( |50| ). The pres- 
ence of new orbits, compared to the case with full leak, 
is clearly seen when comparing Figs. [6] and [9] (c). Fig- 
ure [9jc) corresponds to the c-density p c (x) of the partial 
leak case. The chaotic saddle (CS) is shown in Fig. (9^b) 
and should be compared to Fig. [4^b)- Again, it can be 
thought as CS = W s H W u , i.e., the points that remain 
with J > for t — >• ±oo. The natural extensions of the 
concepts of chaotic saddle and its invariant manifolds are 
expected to hold in all systems with partial leaks. 

A complete description of distributions such as those 
shown in Fig. [9] cannot be achieved by simply includ- 
ing the new trajectories mentioned above. For instance, 
trajectories that collide infinitely many times with R < 1 
but still less frequently than other trajectories might also 
contribute to the asymptotic properties. Indeed, it ap- 
pears natural that the support of the distributions such 
as those shown in Fig.[9]will fill an area of the phase space 
and therefore have a trivial fractal dimension Dq = 2. 
Given the peaked structures of the densities, we can spec- 
ulate that the information dimension Di and the gener- 
alized dimension D q w ith q > might be below 2 [see 
(Tel and Gruiz 2006) for simple examples of sets with 



trivial Dq and non-trivial D\\. It remains, however, to 
be verified whether relations (28)-(30) are applicable to 
systems with partial leaks. 

In the next section we introduce a formalism based 
on Perron-Frobenius operators, that directly calculates 
the (normalized) density p c and the artificial (non- 
observable) distinction between J and p c becomes su- 
perfluous. 




FIG. 10 (Color online) Schematic illustration of flow and map 
coordinates. We use the convention that a particle in the flow 
(■) has coordinates (s,p,r) = (x, r). Smaller dots represent 
other particles in the flow who move along the same trajec- 
tory. The true-time map coordinates are denoted by • at the 
boundary of the billiard, and are characterized by coordinates 
x (s,p). 

Before discussing how densities evolve in time, it is 
essential to clarify the relationship between the different 
invariant densities in billiard systems. Here we search for 
exact relationships between the measure fip of the flow 
and the measure \i of the true-time map (for more de- 
tails see Sec. VIII.B). In closed Hamiltonian systems the 



natural measure /i is the Lebesgue (or Liouville) measure. 

A trajectory of the billiard flow can be represented by 
the Birkhoff coordinates x = (s,p) of its next collision 
with the boundary and the time r after th e previous col- 
lision, as i l lustra ted in Fig. [To] (see Ref. (Chernov and 
Markarian, 2006) for a rigorous mathematical formula- 



tion). Since the velocity of the particle remains constant 
in a billiard, it can always be chosen to be unity. Time r 
also corresponds to a distance in the billiard, and we call 
it therefore the coordinate time. Of course, the coordi- 
nate time fulfills > r < t C oii( x ) where the collision time 
^coii( x ) is defined as the time t between the collision in x 
and the last collision 8 . 



These definitions are convenient because they assure that t co \\ (x) 
and the c-measure are defined inside the leak x £ /. 
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A. Closed system 

In the time-continuous representation, let Pf(x, r) de- 
note the probability density for finding a flow trajectory 
at (x, r). The corresponding invariant measure dpp of 
the flow can then be written as 



djip — Pf(x, r)dxdr. 



(51) 



Consider now the true-time map in standard Birkhoff 
coordinates x = (s,p). Let p M (x) denote the probabil- 
ity density of the invariant measure in this map. The 
corresponding measure d\i is 



d\i = p fI (x)dx. 



(52) 



We wish to connect p to pf- Since the dynamics between 
collisions is a uniform motion of unit velocity, pp is inde- 
pendent of r and, therefore, the flow density in variable 
x is proportional to that of the map 



(53) 



where A is a constant. Its value follows from the normal- 
ization of the measure J dpp — 1: 

/r /**coii(x) 
Pf(x, r)drdx = / / Ap fI (x.)drdx 



A / t co n(x)p AX (x)dx = 1. 



(54) 



Since the integral in the last equality equals the average 
collision time in the closed system's true-time map, we 
have 



A 



1 



(*coll) 



(55) 



From this relation, and Eqs. ([51]), (52), and (55), it fol- 



lows that the flow measure and the map measure are 
related as 



dp F = 



dr 

(^coll) 



dp. 



(56) 



This is a well known relation connecting flows and maps 



in clo sed billiards (see, e.g., (Chernov and Markarian 
|2QQ6| )). Recalling that (t coll ) can be related to simple 
geometric properties of the billiard as in Eq. (14), this 



formula shows that the mean collision time provides a 
simple and elegant connection between properties of the 
flow and map. For instance, Lyapunov exponents of the 
flow are equal to those of the true-time map divided by 

(*coll)- 



From Eq. (56) one can also estimate the escape rate 



in a system with a small leak by assuming that the ex- 
istence of such a leak does not influence the validity of 
this relation. The rate of change of pp over coordinate 
time r is in this type of perturbation approach 



dpF 
dr 



dp 

(^coll) ' 



(57) 



and this rate is independent of coordinate time r. The 
fraction of trajectories escaping through leak / over time 
unit can be computed as J 7 (I) = Jj dxT . Consider a time 
interval of the length of the average collision time (t co ii). 
The decay of surviving trajectories is exponential, P(t) = 
e~ Kt , for any time t for trajectories distributed according 
to the c-density, which is then approximately p^ . We can 
thus estimate k using Eq. ( [57| ) and writing the proportion 
of particles that escape the billiard in a time t = (£ co ii) 



as 1 



-K>(t co n) . 



J(7)M= /1 (/) = l-e- 



-«(*coll) 



This leads to the naive estimate k* = 



(58) 
ln(l - 

p(I)) / (t co \\) stated in Eq. ([24). 

Finally, the dynamics of densities in closed maps can 
be re-written in terms of the Perron-Frobenius opera- 
tors QDorfman[ |1999| |Gaspard| |1998| |Lai and Tel] |2Q11[ ) 

Pn(x) 



Pn+l(x') 



(59) 



K6f _1 (x / ) 



where is the Jacobian at point x. In Hamiltonian 

systems J7(x) = 1, and the constant Lebesgue density 
p M = 1 is a stable fixed point of Eq. ( [59] ) . Since there is 
no escape, the largest eigenvalue of the Perron-Frobenius 
operator is unity. In the next subsection we repeat the 
procedures presented above for leaking systems, i.e., we 
derive a connection between flow and map measures, an 
expression for the escape rate, and establish a Perron- 
Frobenius formalism. A direct connection between the 
Perron-Frobenius operator of the open and leaky systems 
has recentl y been investigated for cases wi th Markov par- 
titions by (Froyland and Stancevic 2010). 



B. Flow and map measures in billiards with leaks 

From the point of view of the escape process, the ana- 
logue of the natural measure is the c-measure for both 
the flow and the true-time map. It is essential to use 
the true-time map that properly takes into account the 
real time of trajectories. Indeed, it is known that av- 
erages taken with respect to flows (or true-time maps) 
and to t raditional maps differ considerable in open sys- 
tems ( |Kaufmann and Lustfeld| |2001[ ) . 

The main difference with respect to the closed case is 
that the c-measure pp c of the flow is not uniform in the 
coordinate time r, rather it decays as e~ Kr . Therefore, 



Eq. (53) has to be replaced by 



p Fc (x,r) = A c p c (-x)e 



(60) 



where pp c is the density of the flow's c-measure, and p c 
is the density of the true-time map c-measure, p c . This 
latter is the measure we have used in the main part of 
this paper 9 . Proportionality factor A c again follows from 



In view of the convention in Fig. llOl pfc( x t t — 0) is actually 
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J dfiFc — 1 as 

Pi? c (x, r)dr<ix = I I A c p c (-x)e~ Kr drdx 



= A cJ Pc(x)^(l -e~ Ktc ^)d^= 1. (61) 



This yields 



l-(e- 



' ^^coll 



(62) 



which tends to l/(t C oii) for ^ ^ as in Eq. (55). The 
map and flow c- measure fip c and fi c are related as 



djii 



A c e Kr drdja c . 



(63) 



This is the generalization to open systems of rela- 
tion ( [56] ) . In contrast to the case of closed billiard, the 
relation between flow and true-time map is not given by 
(^coii) ( or (^coii)c) alone, but includes an explicit depen- 
dence on k and the coordinate time r. In Sec. IVIIIAl we 
also connect \±f c and \i c to the flow measure projected to 
the billiard's boundary. 



C. Escape rate formula 

We can now obtain exact formulas for the escape 
rate following the arguments used for closed systems at 
the end of Sec. |III.A From Eq. (63) one can again in- 
troduce the rate of change over coordinate time r as 
T c = dfiFc/dr which is, in contrast to Eq. ( |57| ), not 
independent of coordinate time r. The fraction of tra- 



jectories that escape through leak / in time {t co \\) c is 



J^ tcou)c T c (f(I))dr. In view of (|62|) and (|63[), Eq. (|58) 
turns into 



I 



(tcon)c 



T c (f(I))dr = /i c (/(J)) 



1 



-K(t c , 



(e~ KtcoU ) c 



1-e 



-K(t c , 



(64) 

Here, in the last equality we used again that for the c- 
measure the fraction of trajectories escaping up to time t 
is 1 — exp (—nt). Thus, using also the convention that the 
escape happens immediately (i.e. t co \\ (escape) = and 
thus /i c (/(/)) = /i c (/)), we obtain an implicit relation for 
k as: 



(e-^) c = l-fi c (I). 



(65) 



the map density p c right after the collision at the boundary point 
/ — 1 (x), but this density remains unchanged until the collision at 
x. The flow density pp c decays, however, in r until the boundary 
point x is reached at r = t C oll( x )- The e~ Kr dependency of 
pF c on r guarantees the stationarity of p c in the discrete time 
dynamics. 

10 In view of the footnote after Eq. ( |60| , we use the image of the 
leak /(/) = "escape" to compute the flow at the boundary point 
x 6 I. 



This formula, which is a new result and will be derived 
more formally in the next section, establishes a relation 
between the average of an expression containing n times 
the collision time and the measure of the leak, both aver- 
ages taken with respect to the c-measure of the true-time 
map. It is thus a generalization of the Pianigiani-Yorke 
formula (33), valid for usual maps, which is recovered 
from (65) in the limit of t co u 



1. 



It is worth applying the cumulant expansion to the left 
hand side of (65): 



ln((e" 



-Kt co n \ 



=) = E 



(-«)' 



Cr(t co u) = ln(l - Mc(/)), 



(66) 

where C r (y) are the cumulants of y (C\ = (y) c ,C2 = 
cr?(?/), etc.). Keeping only the first term of the expansion, 
we find a first order approximation, k\, of the escape rate 
as 

ln(l - p c (I)) = -fti(t C oii) c . 
This is Eq. ([32]) which was obtained as an improvement 



to the naive estimate Eq. (24) based on the theory of 
transient chaos and has been used in Sec. [Til 

We can estimate the deviation between ki and k by 
including the second term of expansion (66). For the 
second order approximant a quadratic equation is ob- 
tained: 



ln(l - fjL c (I)) = -tt 2 (t C oll) C 



^2 2 



^coll,c ' 



where a? , is the second cumulant of the collision time 

kcoll,c 

distribution t co n(x) taken with respect to the c-measure 
fi c . This yields the explicit form for ^2' 



^2 



coll/ 



<tcon)? + 2<x t 2 colliC ln(l- Mc (/))]/<7 t 2 collc 



(67) 

Similar corrections due to the uneven distribution of t co \\ 



have been obtained by ( jMortessagne et al. , 1992), see 
also (Joyce 1975). 



The approximation ft 2 ~ fti is valid for 

T 2 



t coU ,c 

(icon) 2 



ln(l - Mc (/)) < 1. 



(68) 



Our numerical simulations (see Appendix VIII. B) yield 



independent estimations for n and In our typical 

configuration (e.g., Fig. [5J, the agreement between k,i 
and ^2 was on the order of 0.6% and between K2 and 
k of 0.06% (below the precision of the results reported 
in Tab. [ll|). From Eq. ( [68] ) we see that there are two ef- 
fects that can make k\ to be a good approximation of k\ 
(a) small leak: fi c (I) <C 1, and (b) cr tcoluC / (t co n) c <C 1. 
The latter corresponds to a small variance of the colli- 
sion time distribution, and is the limit under which the 
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true-time map reduces to a usual map [and Eq. (33) is 
recovered]. Note, however, that the naive estimate ft*, 
Eq. (24), does not follow in any order of the cumulant 
expansion. As previously noticed by ( |Joyce| 1975, 1978), 
ft* coincides with ft only under the physically unrealistic 
conditions that: (i) the distribution in the room is uni- 
form p c = e.g. non-specular collisions immediately 
randomize the trajecotires; and (ii) t co \\ is constant. If 
only condition (i) is satisfied, Eqs. (65) and (67) could 
be applied using p instead of p c . Joyce also pointed that 



an analogue of Eq. (67) was published by Kuttruff in the 



1970's. The ft* estimate (24) is the result of a not fully 



consistent argumentation: it intends to correct for the fi- 
nite size of the leak, but fully forgets about changing the 
measure to the c- measure, which would provide a compa- 
rable other correction. The naive estimate ft* is obtained 



in (Bunimovich and Dettmann, 2007) as the first term of 



an expansion in the leak size. 
Partial leaks 

We are now in a position to heuristically generalize 



(65) for systems with partial leaks R ^ discussed in 



Sec. II. D The key observation is that the left and right 
hand sides of Eqs. (65) can be interpreted as global and 



local quantities, respectively. The global quantity in the 
partial leak case remains unchanged: in one iteration of 
the true time map the c-density at position x decays as 
e -«tcozz(x)^ anc L gi ODa i estimation of the proportion S 
of the survivors after a time t co \\ is obtained simply as 
the average of this factor over the full phase space Q 



S = [ p c (x)e-^ (x W= {e 
Jq 



Kt co n \ 



(69) 



We denote the c-density of problems with partial leaks 
as p c . Index c refers to such c-densities. The same quan- 
tity is obtained in a local approach by considering the 
proportion of particles not crossing the leak in one iter- 
ation of the true-time map. In the full leak case this is 
1 — J / p c (x)(ix = 1 — p c (I) which appears on the right 



hand side of Eqs. (65 ). When the leak is partially reflect- 



ing, the portion JJ~R(x)p c (x)<ix of the leak's measure 
enhances the number of survivors. The full proportion S 
of the survivors is 



1 ~ Jj Pc(x)dx + Jj i?(x)/5 c (x)dx 
l-/ 7 (l-^(x))p c (x)dx 

<#(x))c, 



(70) 



where we used that the last integral [which is exactly 
m c (I) used in (48)] can be written as an integral over 
the full phase space — i?(x))p c (x)<ix, since out- 

side the leak R = 1, and we used the normalization 
J^p c (x)(ix = 1. This means that the fraction of sur- 
vivors is the full phase space average of the reflection 
coefficient (R)c- Equating the global ( [69] ) and local (70) 
expressions for the survivors we obtain the formula 



(e- kt ™") £ = (R) e 



(71) 



This is a general exact relationship for ft that applies to 
all cases discussed in this paper and is one of our main 
new results. In the next section we provide a derivation 
of this expression based on Perron-Frobenius operators. 
Equation (71) goes over into relation (65) for full leaks, 
i.e. for R = inside the leaks. 



It is worth returning, in view of (71), to Sabine's orig- 



inal problem (in the sound ray approximations) for a 
closed hall: a system with very weak partial leaks. The 
escape rate is then small, averages can be taken with re- 
spect to the closed system's measure, and ( 71 ) goes over 
into 



l-ft&oii) = (R). 



(72) 



Using ([14]) for a 3D billiard with the sound velocity as 
, (1 - (R))Ac 



speed, we find 



(73) 



Thus the effective area appearing in Sabine's formula ^ 
can be interpreted for a closed hall a posteriori as 



AA = <(1 - R))A, 



(74) 



i.e., the average of the transmission coefficient T = 1 — R 
taken with respect to the Lebesgue measure times the full 
surface area. The discovery of the concept of reverbera- 
tion time can thus be considered as the first observation 
of a closed chaotic billiard with partial leaks extending 
over the full surface. In such cases only energy is escap- 
ing the system. 

If, in addition, doors of small surface area AAq are 
also opened, sound rays also escape and this expression 
is modified as 



AA = AA + (1 — R)(A — A ) 



(75) 



It is remarkable that the problem of partial leaks shed 
now light on Sabine's result more than hundred years 
after the original work. 

It should finally be noted that the value of our es- 
cape rate relations (71) and (65), just like that of the 



Pianigiani-Yorke formula (33), is conceptual. They do 



not provide an efficient way for determining the escape 
rate (that is very well handled numerically), rather they 
illustrate how the escape rate of the flow follows from 
properties of the true-time map. 



D. Operators for true-time maps with partial leaks 

Consider an invertible open map f . Its escape rate 7 is 
known to appear as the largest eigenvalue of an operator, 
the Perron-Frobenius operator. This operator is defined 
(|Lai and Tel} |2011| [Till |1987[) by the iteration scheme of 



a density function p 

e" 7 Pn+i(x) 



Pn(x) 



(76) 
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where J7(x) is the Jacobian at point x. By considering 
the right hand side to be the result of an operator acting 
on function p, the left hand side shows that e -7 is an 
eigenvalue (in the space of positive p-s, the largest eigen- 
value) of this operator. Equation ( 76 ) expresses that the 



total probability in a small region at step n is the same 
as in the image of that region under map f , when tak- 
ing into account a factor e 7 for compensating the escape. 
The escape rate follows from the requirement that the 
integral of p n over a fixed phase space region containing 
the chaotic saddle remains finite in the limit n — » oo. The 
limit distribution p^ is then the density p c of the condi- 
tionally invariant measure concentrated on the unstable 
manifold of the chaotic saddle. This is the well-known 
picture for open maps ( |Lai and T6\\ |2Q11| [Till |1987[ ) . 

As a generalization of this idea, in a true-time prob- 
lem where the distribution of collision times t co u(x.') is 
known (remember, the collision time is associated with 
the later collision point), the continuous-time escape rate 
K is determined by the iteration scheme 



-Kt co n(x'] 



Pn+l(x') = 



Pn(x) 
I^(X)|| 



(77) 



xef-^x') 



In this equation the true-time property is incorporated in 
the left hand side in the multiplicative factor containing 
the escape rate, in a similar fashion as in Eq. (60) and 



consistent with our convention for t co \\ in Eq. ( 19 ). This is 



a higher-dimensional extension of the true-time formal 
ism used for o ne-dimensional maps by ( |Kaufmann and 
Lustfeld 2001) and, in a slightly differe nt context, has 
been previously considered by (Gaspard, 1996, 1998). 

The escape rate n again can be considered as an eigen- 
value, and its value follows from the requirement that the 
integral of p n over a fixed phase space region containing 
the chaotic saddle remains finite in the limit n — » oo. The 
limit distribution p^ is the density p c of the conditionally 
invariant measure of the true-time map. 

The problem of maps with complete or partial leaks 
can be treated as a further generalization. Since the dy- 
namics of trajectories is then closed, we write the map as 
f 'closed- The escape rate k for partial leaks follows from an 
iteration scheme in which the reflection coefficient i?(x) 
also appears, in a similar spirit as in Tanner's work (Tan- 



ner 



be seen as partial leaks with R(x) =0 for x G /. 



For traditional maps with leaks, t co ii(x) = 1, Eq. (78) 
yields the escape rate 7 of such maps: k \-> 7. It is inter- 
esting to see from (77) and (78) that in closed systems 
(R = = 0) the Perron- Frobenius operator coincides 
with the classical form (59). There is then no essential 



difference between the true-time and the map picture due 
to the simple proportionality of the flow's and the map's 
measure as expressed by ( 53 ) . 



The operator formalism for true-time maps with par- 
tial leaks developed above is a new result of our paper. It 
unifies and generalizes all previous approaches. In fact, 
Eq. ( 78 ) suggests that the physically rather different phe- 



nomena of collision times and of reflection (without tak- 
ing into account collision times) are described by essen- 
tially the same mathematical mechanism: the density 
p(x) in Eq. (78) is multiplied by a function of the phase 



space coordinates x. As a consequence, the inclusion of 
the collision time distribution has similar effects on the 
iteration as a reflection coefficients larger than unity. 

It is now straightforward to obtain an exact relation 
for k. Consider the c-density p c (x) of a system with full 
or partial leaks, i.e., the limit distribution of Eq. (78). 



By integrating both sides over the full phase space, we 
obtain 



,i?(x)p c (x) 



(79) 



= f Q dxR(x)p c (x.) = (R)c, 

where we used |j7(x)| = |dx'|/|dx|, and /dosed (^) = ^- 
This provides a proof of our fundamental formula (71 1 11 . 



Although Eqs. (71) and ( |79| ) are identical, it is worth em- 
phasizing that (71) was obtained from a qualitative ar- 



gument based on properties of billiard dynamics. Since 
Eqs. (77) and (78) are valid for any true-time map, the 



derivation presented here shows that Eq. (79) is not re- 



stricted to billiards, it holds for leaky systems in general. 
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In our notation, R shows up on the right hand side since 
there is an immediate loss of density wherever R is dif- 
ferent from unity. If collision times are also taken into 
account, we find 



e-^°" (x 'Wi(x') 



■R(x)p n (x) 

U(x) 



(78) 



In order to illustrate our formalism in a simple ex- 
ample, we consider area-preserving baker maps. This is 
motivated not only by the possibility of an analytic treat- 
ment, but also by the fact that several papers ( Ermann 
et al\ |2012| Keating et al\ |2006| |Nonnenmacher and 



I. 



The limit distribution p^ is the density p c of the con- 
ditionally invariant measure in the true-time map of the 
partially leaking system. One might also have leaks in 
a naturally open system, in which case the open map f 
should be used in the relation above. Full leaks can also 



11 Had we defined the collision time as a function of the initial 
coordinate as t C oll( x ) — ^coll(/( x )) 5 it is natural to shift the 
factor e -Ktcoll ( x ) to the right hand side of Eq. < [t8^ , and the 
same argument that led to Eq. ( [tq] ) leads to (i?e Ktcoll )c = 1. For 
the full leak case, Eq. ( pi] becomes (e Ktcoll )c = 1 + l^c(I)- 
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Zworskil 2005 Novaes et al t 2009 Pedrosa et al, 2009) 



use leaky versions of these m aps to investigate quantum 
systems (see Section VI.E.2| ). 

First we consider the triadic area-preserving baker map 
defined on the unit square (x, y) G [0, 1] x [0, 1]: 



(Xn+liVn+l) = (lx n ,3y n ) 



for y n < 1/3, 



(x n +i,y n + 1 ) = (\(x n + l),3y n - 1) for 1/3 < y n < 2/3, 
(x n +i,yn+i) = (l(x n + 2),3y n - 2) for y n > 2/3, 

(80) 

with a very simple choice of the leak: I is a band of 
height 1/3 in the expanding (y) direction in the middle 
of the square. In this model there are no partial leaks 
present. The collision times takes on two values only: 
^coii = r i if the point is mapped into the column x < 1/3 
and t co n = T2 if the particle is mapped into the column 
x > 2/3. 

Let us start from a constant distribution po = 1 on the 



unit square. Since the Jacobian is unity, Eq. (77) tells us 



that the measure (under p ) of the lower band y < 1/3 
is mapped on the column x < 1/3 with the measure on 
it multiplied by e KTl . Similarly, the measure of the right 
column 2/3 < x < 1 will be e KT2 /3. The measure from 
the mid band is mapped nowhere because it is in the leak. 

This construction extends to finer scales in a self- 
similar manner due to the simple choice of the leak. The 
measure converges to the c- measure, which therefore has 
to remain (conditionally) invariant under the above iter- 
ation. The value of the escape rate therefore follows by 
prescribing the invar iance of the c-measure. In this case 
it is sufficient to consider the one-time measure projected 
on the x axis because the density is constant along y. By 
construction, this measure is originally 1, and after one 
step it is the sum of the two values just determined, and 
thus 



3. 



(81) 



This is an (irrational) equation for the escape rate n. The 
validity of the new formula (65) can be easily verified: 



( e -«t coU ^ 



2/3, 



1, Eq. (81) 



which corresponds to 1 — p c (I) since the leak has height 
1/3 and p c is independent of y. 

In the limit of traditional maps, t\ = T2 
yields n — >• 7 = In (3/2). Since the stretching rate is 
everywhere 3 in the phase space, the maps average Lya- 
punov exponent is A = In 3. The Kantz-Grassberger for- 
mula (30) then yields D±~\l) = In 2/ In 3, i.e., the unsta- 
ble manifold of this leaky baker map carries the structure 
of the classical triadic Cantor set. 

As a more complex example, consider the dyadic baker 
map 



(x n +i,y n +i) = (^x n: 2y n ) 



for y n < 0.5, 

(xn+uyn+i) = (1 - ^~2 n \ 1 - 2(1 - y n )) for y n > 0.5, 

(82) 



with collision times taking on again two values only: 
^coii = r i if f° r the image point x' < 1/2, and T2 other- 
wise, but with partial leaks. These leaks are introduced 
with reflection coefficients R4 on four horizontal 



strips of height 1/4, as Fig. 11 illustrates. 
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Ri{ 
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FIG. 11 (Color online) Illustration of the dynamics of a 
dyadic true-time area-preserving baker map with partial 
leaks. The right square is obtained applying baker map (82) 
to the left panel. 

Let us consider again po = 1. After one iteration one 
finds that four different c-densities appear in the four 
quadrants of the unit square. A next iteration refines the 
picture, but one feature remains: the c-measure coarse 
grained on these four quadrants has 4 different densities. 
The values on the 4 rectangles change as iteration goes 
on, but the four-value structure remains unchanged. 

In order to find an analytic expression for the escape 
rate, it proves to be sufficient to deal with the course 
grained c-density p c {y) projected on the expanding (y) 
axis. As follows from above, after a large number of 
iterations, this c-density is piecewise constant and has a 
jump at y = 1/2. Let us therefore assume its form as 



Pc(y) = 



c for < y < 0.5 
for 0.5 > y < 1, 



(83) 



which fulfills normalization p c (y)dy — 1 for any c. The 
values of c and k follow from the requirement that this 
projected measure remains invariant after one more time 
step. 

The lowest horizontal strip of height 1/4 (and of c- 
measure (2-c)/4) is mapped, in view of ( |78| ), on the 
bottom left quarter of the square with a new measure 
exp(/cTi)i?i(2 — c)/4 (note that the Jacobian is unity). 
The second horizontal strip is mapped on the top left 
quarter with c-measure exp(£ri)i?2(2 — c)/4. The third 
and fourth strips come into the remaining quarters with 
weights exp(^T2)i?3c/4 and exp(&T2)i?4c/4, respectively 
(the last one representing the bottom right quarter). 

The criterion of the invariance of the projected mea- 
sure is that the total c-measure in the two bottom (top) 
quarters is the same as the integral of (83) over < y < 



1/2 (1/2 < y < 1). Thus, we find two equations 



| = I[ e - 1 i? 1 (2-c)+e^i? 3 c], 



-[e~ KTl R 2 (2 - c) + e' KT2 R 4 c}. 



(84) 
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After rearrangement, we obtain 

4 - 2e* Tl i?i - 2e kT2 R 4 + e^+^i?^ - R2R3) = 0, 

(85) 

and 



A-2e kri R 1 



2 - e kT ^Ri +e kr iR 3 



, (86) 



Equation (85) is an implicit equation for the escape rate, 



while (86) provides the value of c determining the jump 



in the c- measure (83) projected onto the y axis. We have 
verified that k obtained from such implicit relations agree 
with direct numerical simulations of the baker map in- 
cluding the intensity J n and real time t n . 

To have more analytic insight, now we focus on the 
particular case of t\ = r 2 = 1 (map with partial leaks). 
Equation (85) leads then to a quadratic expression for 



escape rate 7: 



4 - 2(R 1 + # 4 )e 7 + (R1R4 - #2#3K 7 = 0. (87) 

As a simple particular case, let us assume that there 
is a full leak over the uppermost horizontal band, and no 
partial leak anywhere: R4 = 0, Ri = R2 = R3 = 1. From 
Eq. ([871) and pfl) we obtain 



e 7 = v / 5-l, c = 3-V5. 

Since the average Lyapunov exponent in this uniform 
baker map is A = In 2, the stable manifold's information 
dimension is, in view of (T30l) 



D 



00 



ln(x/5- 1) 



In 2 



clearly below 2. The numerically generated unstable 
manifold of this map can be seen in the left panel of 
Fig. 12, and confirms the fractal property. The right 



panel shows the projected c- measure which clearly ex- 
hibits a jump at y = 1/2. The plateau values agree well 
with the theoretical predictions. 

Had we taken the leak with R\ = (R2 = R3 = R4 = 
1), the same result would have been obtained. The situ- 
ation is different, however, for R2 = or R3 = 0. In the 
first case, we see from (86) that c = 0, implying vanishing 



c- measure for the entire upper half-square. The escape 

(s) 

rate is then 7 = In 2 which implies D\ } = 1. Fractal- 
ity is then lost, the stable and unstable manifolds are 
one-dimensional. The dynamics is fully leaked, there is 
no chaos, a single unstable fixed point, the one at (0,0), 
governs the escape dynamics (hence 7 = A). For R3 = 0, 
c = 2, the c-measure vanishes in the bottom half square, 
the situation is otherwise the same. These observations 
illustrate again that the location of a leak of the same 
area is very important, even chaos can be lost if they do 
not overlap with one period- 1 orbit. 

Now, we consider the more general case of a partial 
leak of arbitrary reflection coefficient over the uppermost 
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FIG. 12 (Color online) Numerical results for the dyadic baker 
map with a full leak / over the uppermost strip of height 1/4: 
Ra = 0, Ri = R2 = R3 = 1, T\ = T2 = 1. Left panel: unstable 
manifold obtained as the endpoint of particles surviving up 
to n — 30 iterations. Right panel: Distribution of the c- 
density projected on the y axis p c (y)- The straight lines are 
the analytical results: 2— c = 1+v^ for y < 0.5 and c = 3— 
for y > 0.5. 



strip: R 4 < 1 (R x = R 2 = i? 3 =l). The solution of ([87]), 
(86) leads to 



R 4 - 3 ± y/R\ - 2R A + 5 



Ra 



1 



We took here the + root because for the — root the mea- 
sure 2 — c for y G [0, 0.5] would be negative. Using the 
result above and Eq. (f8TJ) with R\ = R2 = R3 = 1 it is 



not difficult to confirm the validity of the formula (71) 
for this partial leak case. As expected, for R4 — > 
obtain c = 2 — c = 1 and 7 = 0. 



1 we 



IV. IMPLICATIONS IN FULLY CHAOTIC SYSTEMS 

A. Dependence of the escape rate on the leak 

The main message of the previous section is that for 
systems with finite leaks the naive theory based on the 
close dynamics differs from the correct theory based on 
open systems. One of the most striking and best stud- 
ied effects arising due to this difference is the depen- 
dence of the escape rate on the position o f a fixed- 
size le ak (|Afraimovich and Bunimovich , 2010 Altmann 
et all 120041 IBunimovich and Yurchenkol 1201 II IBuni- 



movich and Dettmann[ |2007| |Demers and Wright [ |20 11 



Paar and Pavin 1997 Schneider et al. 



2002). This re- 



sult is shown in Fig. 13 for the cardioid billiard with a 
leak / = [si — As, si + As] x [pi — Ap,pi + Ap] with 
fixed size As = 0.1, Ap = 0.2,^ = 0, and different 
positions s/. Physically this type of leak, illustrated 
in Fig. 13'a,b), could be realized in optical systems by 
replacing the perfect mirror boundaries in the region 
[s/ — As, si + As] by dielectric material with refraction 
index n = 1/ sin(Ap). The results in Fig. 
the non-trivial dependence of k on s^, w: 



13 



c) confirm 



mich can take 
values both smaller and larger than the naive estima- 
tion k* = -ln(l - AsAp)/(t coll ) given by Eq. p4|. The 
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FIG. 13 (Color online) Dependence of the escape rate on the position of the leak, (a) Cardioid billiard with a leak centered 
at position si,pi = 0. (b) Phase space representation showing the leak I = [si — As, si + As] x [pi — Ap,pi + Ap] with 
Ap = 0.2, As = 0.1, pi = 0, and si £ [0, 1]. (c) The escape rate k obtained from numerical simulations as the ones in Fig. [I] 
(d) The c- measure of the leak /jl c ( I) and (e) the inverse of the mean collision time. Results in (d) and (e) were calculated from 
p c (s,p) as described in Appendix |VIII.B Using these values to compute k through Eq. (32) leads to results undistinguishable 
from those of panel (c). Horizontal red line in (c)-(e) corresponds to the estimate based on the closed billiard theory of Sec. II. A 
in (c) K* given by Eq. (24), in (d) = AsAp, and in (e) (t co ii) given by Eq. (14). The y axis on the right edge of panels 
(c)-(e) indicate the relative deviation between these values and the actual data. 



0.01 



0.008 




the l eak around short periodic orbits of the system (see 



FIG. 14 Dependence of the escape rate 7 on the location of 
the leak in the fully chaotic logistic map: x n +i = 4x n (l — x n )> 
Leaks have = 1% and are symmetrically centered at 

xi. The positions of the unstable periodic orbits (UPOs) 
of period 1 and 2 are indicated by arrows. See also (|Alt 
mann et al.\ [2 004; Bunimo vich and Yurchenko[ |2011| |Paar 
and Pavin[|1997) . 



theory developed in Sec. |ILB| tells us that n depends on 
Hc(I) and (t C oii)c through Eq. ( |32| ). Figure 13'd,e) show 
that both factors p c {I) and (£ C oii/c contribute to k ^ n* 
but the variation of p, c (I) is the stronger factor in the 
dependence of n on s\ (at least for this size of the leak). 

Apparently the first to report the dependenc e of the 
maps escape rate 7 on leak position were (Paar andl 



Pavin, 1997). The most pronounced effect shown for the 
logistic map, which is reproduced in Fig. [l4j was the 
relationship between small values of 7 and positions of 



also (Altmann et al t 2004 Bunimovich and Yurchenko 
|2011[ )|! This effect is less pronounced but also visible 
in Fig. [l3j where local minima of the n vs. s\ curve 
are obtained when the leak is placed around the low- 
est periodic orbits of the billiard: the horizontal orbit at 
(s = 0,p = 0) H> (s = — l,p = 0) and the vertical orbit 
(s = -0.5, p = 0) ^ (s = 0.5, p = 0). 

An intuitive explanation of these results is found by 
looking at the images and pre-images of leak / (Buljan 



and Paar 2001 Paar and Pavin 1997). The surviving 



trajectories at time t correspond to all trajectories that 
are not in any of the time t pre-images of /. Now, if 
pre-images overlap repeatedly, there are more surviving 
trajectories for increasing t, and therefore the escape rate 
is smaller. The overlap between the leak and its im- 
ages (or pre-images) is obviously increased when the leak 
is around periodic orbits. We have already learned in 



Sec. |II.B| that the escape rate can be calculated in terms 
of the periodic orbits inside the leak through Eq. (36). In 
this formalism, orbits with low period play an important 
role, with weights inversely proportional to their insta- 
bility. Therefore we can expect that if short orbits with 
small expansion rates fall into the leak, i.e., they appear 
in the sum (36), the value of 7 is reduced. 



This qualitative argument has been rigorously ex- 
tended for different classes of fully chaotic systems, such 
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FIG. 15 Dependence of the escape rate on the size AsAp of the leak, (a) Cardioid billiard with a leak centered at si = 0.4, pi = 
and variable size, (b) Phase space representation showing the leak / = [si — As, sj+As] x [—Ap, +Ap], with Ap — As G [0, 1/ a/5]- 
(c) Escape rate ft as a function of AsAp. (d) The c-measure of the leak fi c (I) and (e) the inverse of the mean collision time. 
Results from (d) and (e) were calculated from p c (s,p) as d escrib ed in Appendix VIII. B Black symbols: numerical results. Red 
lines: estimate base d on the closed billiard theory of Sec. II. A in (c) ft* given by Eq. (24), in (d) = AsAp, and in (e) 



(tcoii) given by Eq. (14). Blue dashed line: relative deviation between black symbols and red line, see the y axis at the right 
edge of panels (c)-(d). 



as ID expanding linear maps ( Afraimovich and Buni- 


movich 


2010; Bakhtin and Bunimovich[ 2011 




Buni- 


movich 


, 2012| Bunimovich and Yurchenko[ 2011 




^ergu- 


son and Pollicott, 2012 Keller and Liver ani, 


2009). An 



important feature in these approaches is that all systems 
admit a Markov partition and the holes are chosen to 
coincide with one element of the partition. For small 
leak sizes, a fr actal dependence of the escape rate 7 has 
been observed ( Bunimovich and Yurchenko"! 2011 Knight 



et al.\ |2012[) (see also (jAltmann and Endler||2010D). As 



noticed in (Bakhtin and Bunimovich, 201 1| | Bunimovich] 



and Yurchenko| 2011), the position dependence of 7 
might be so strong that there are cases in which a hole 
two times larger than another one can have a smaller es- 
cape rate. These results can also be understood in the 
exact expansions developed and applied to fully chaotic 



billiards with finite but small leaks in (Bunimovich and 



Dettmann[ |2007| |Dettmann| |2011a[ ). 



Beyond the Markovian approach, an alternative expla- 
nation based on recurrence times and Kac's lemma ap- 



pears in Ref. ( Altmann et al. , 2004). While a great num- 



ber of analytical results can be obtained in ID Markov 
systems, the results shown above in the Cardioid billiard 
(area-preserving map) appear to show a smoother depen- 
dence on position. 

Finally we would like to mention that similar strong 
dependencies are observed also by varying other parame- 
ters of the leaks such as the orientation of (asymmetric) 



leaks, as observed for Hamiltonian systems in (Schneider 



et 



2002). Interestingly, similarly complicated even 



fractal dependencies with the position of the leak have 



been observed in the diffusion coefficient (Klages 2007 



Knight et al. 2012). Exponential decay and leaks with 
different shapes were investigated in detail in the periodic 



Lorentz gas in (Demers et al. 2010) 



Complementary to the position-dependence of ft, the 

dAir 



dependence on the leak size is illustrated in Fig. 15 



mann et al. , 2004 Bunimovich and Dettmannl 



2007F 



Schneider et al. , 2002 ) . We change the size of the leak 
I = [si - As, s] + As] x \pi - Ap,pi + Ap] by chang- 
ing Ap = As at a fixed si = 0.4,^ = 0, as illustrated 
in Fig. [l5^ a,b). The dependence of ft on As depicted 
in Fig. |15[c) follows roughly the dependence of ft* (red 
line in panel c) but a non-trivial behavior is observed 
apart from this trend. Again, both the measure /i c (^) 
shown in Fig. [l5^d) and the mean collision time (t co \\) c 
shown in Fig. |15[e) are clearly different from the closed 
system prediction (red lines). The oscillations in /i c (i) 
are stronger than the ones in (t co \\) c (at least for the leak 
sizes considered here) and provide the strongest contri- 
bution to ft 7^ ft* . Notice that in view of the dependency 
of ft on the leak position discussed above, one can easily 
find situations in which larger leaks have smaller escape 
rates not only in the case of Markov leaks. 

One particularly important limit is the case of vanish- 
ingly small leaks As, Ap — >• 0. This is the traditional 
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limit mathematicians are interested in (e.g., in the con- 
text of Poincare recurrences) (Haydn et ah, 2005). It is 
interesting to see what happens with the position depen- 
dency of 7 in the limit — >• 0. For the case of the 
doubling map with Markov leaks, the escape rate 7 de- 
pends only on the periodic orbit of lowest period p inside 
the leak and is given in leading order by (see e.g. (Buni- 

" 20091 )) 



movich 2012 Keller and Liver ani 



7=^)(i-n, 

for arbitrary small /i(I). In this limit almost every leak 
position will have p — >• 00 so that Sabine's result is re- 
covered in /i-almost every case. Indeed, also the results 
in Fig. 15'c-e) show that all quantities converge to the 
closed systems prediction (red lines) in this limit. The 
blue line in Fig. [l5^c) shows that the relative difference 
(k-k*)/k* is of the order of 10% for /i(J) « 0.1 but that 
it also consistently decays for small As. Altogether, these 
convergences are physically relevant, particular mani- 
festations of the more general converge nce p c (s,p) — >• 
p fl (s : p) for n(I) — > 0, discussed in Sec. II. B An opti- 
mal j3xpam^on^^ leak s ize was developed 
in ( |Bunimovich and Dettmann 2007| ). Very recently, a 
devil's s taricase dependence of k ( and 7) has also been 
noticed (Demers and Wright, 2011). 



B. Multiple leaks and basins of escape 

The idea of placing more than one leak to the system 
is ver y natural in numerous circumstances (Ble her et al.\ 
1988). Quantum systems often have more than one leak 
due to input /output /transmission/reflect ion/antennas. 
Multiple leaks in chaotic systems have also been consid- 



ered (jBuljan and Paar||2001 Bunimovich and Dettmann| 


2005, 


2007 


Dettmann and Georgiou[ |2011b| 


Portela 


et al. 


2007 


). The results of the previous sections show 



that the c-measure depends sensitively on the position of 
the leak. As a consequence, the results for multiple leaks 
are not only different from the closed system approxima- 
tion, but also cannot be easily obtained from the results 
for each leak alone. For instance, with two leaks 1\ and 
^2, the escape rate is, in general, different from the sum 
of the single leak case: 



K{h+I 2 )^K{h)+K{I 2 ). 



This results is not surprising in view of the non-trivial 
dependence of n on leak size, reported in Fig. [l5j in 
which case it was clear that doubling the leak size 
does not imply doubling k. Indeed, the difference re- 



and Paar 



2001 



ported in Eq. (|88|) has been understood in t erms of 
the overlap of t 



re pre-images of the tw o leak s ( Buljan 
Pikovsky and Popovich| [2003 ), and has 



been developed more systematically by Bunimovich and 
Dettmann. They found that k{I\ + h) can be expressed 
as k{I\)+k{I2) plus a series of correlation terms (with de - 
creasing importance) (Bunimovich and Dettmann, 2007). 



In contrast to the escape rate, the dimensions of the 
invariant sets of the system opened with multiple leaks 
can be estimated from the dimensions of the invariant 
sets with single leaks. Consider the case of two leaks I\ 
and I 2 . The saddle (or the manifolds) obtained for the 
case in which I\ and I2 are simultaneously opened corre- 
sponds to the intersection of the saddle (manifold) when 
only I\ is opened with the saddle (manifold) obtained 
when only I 2 is opened. Very generally, the fracta l di- 
mension Dq(I\ + h) of the intersection is given by (Fal 



1985) 



Do(Ii + 1 2) — Dq(Ii) + A) (^2) — ^embedding, 

where Dq(Ii) (Dq(I 2 )) is the dimension of the set ob- 
tained when only 1\ (I 2) is opened and ^embedding is the 
dimension of the embedding space. For full leaks, the di- 
mension and the escape rate are related to the Lyapunov 
exponent through relation (30). The argument above can 



trivially be extended to more then two leaks. It is not 
valid, however, in cases when the saddles of the two leaks 
are trivially connected to each other (e.g., I2 C h or I2 
is an image of 

The interesting phenomenology of multiple leaks is bet- 
ter illustrated through an example. We consider the car- 
dioid billiard with two leaks, Ii and I2, of the same size 
/i(Ii) = ^(^2) = 0.04 but at different positions, as de- 
picted in Fig. [16] Let us first consider the effect of each 
of these leaks separately. Following the procedures de- 
scribed in Fig. [l3j we obtain the results reported in the 
two first columns of Tab. IIIII 




FIG. 16 (Color online) Escape basins for two leaks in the 
cardioid billiard. Leak 1\ is centered at s\ — —0.5, pi = 
with As — Ap — 0.2, and h at si — 0.6, pi — with 
As = 0.04 and no restriction in the collision angle is applied 
(i.e., Ap = 1). Both leaks have the same area in the phase 
space, nevertheless, 41% of the trajectories escape through 
leak ii (blue region) while 59% of trajectories escape through 
leak I2 (yellow region) when both leaks are open. Different 
characteristic values obtained for this system are reported in 
Tab. EH 



27 



Leaks open during experiment: 


Only h 


Only I 2 


Both I\ and I2 


Measurement (\) applied in leak Ij = (— >) 


h 


h 




h 


h 




0.04 


0.04 


0.08 


0.04 


0.04 


K 


0.01607 


0.02360 


0.03776 


0.03780 


0.03773 




0.02856 


0.04332 


0.06697 


0.02766 


0.03931 


(tcoll)c 


1.8059 


1.8799 


1.8424 


1.8424 


1.8424 




0.02856 


0.04332 


0.06697 


0.06999 


0.06500 


(*coll>^ 


1.8059 


1.8799 


1.8424 


1.9273 


1.7869 



TABLE III Measurements in the system depicted in Fig. 16 While the standard c- measure (in \i c and (£ co ii) c ) considers 



normalization in the full phase space, the superscript f, j in the last two lines indicate that the computation of the c-measure 
was restricted to the set of points (x) which escape through the leak Ij. Obviously, [i\ ,j (Ij) — f^c(Ij) whenever only Ij is open 
(Ij = {ii, I2, I\ + ^2}, compare lines 5 and 7). Note that the measured k through leak ii, I2, and I\ + 12 are consistent to each 
other when both leaks are simultaneously open (three columns on the right). These values are consistent with those obtained 
from the relationship ^1 = — ln(l — fi c ) / {t co i\) c with (Ij), {t co \\)i^ (the values in the last two rows of these columns). For the 
computation procedures, see Fig. [42] (for the f,j cases the set S(t*) was divided in two subsets according to the leak through 
which trajectories escape). Error bars are of order 5 in the last digit. 



Let us now consider the case when both leaks are si- 
multaneously open. The result reported in Tab. |III| con- 
firms inequality ( [88] ) . Similarly, we compute /i c (^)> the 
c-measure for both leaks open calculated at Ii. The mean 
collision time is not defined with respect to a given leak, 
as it depends on p c in the whole phase space, hence the 
corresponding dashes in Tab. |III| 



A natural question is that of the nature of the sets of 
initial conditions which escape through each Zj, i.e., of the 
properties of the escape basins Bi, as the ones depicted in 



Fig. [16} The border between the two escape regi ons con- 
tains the stable manifold of the chaotic saddle dLai and 



Tel, 2011). Notice that the total area of the basins cor- 



respond to the amount of initial conditions that escape 
through each leak, while p c (h) and /i c (^2) are propor- 
tional to the rate of escape through each leak for large 
times. Note that the measured n through leak Ji, I 2l 
and I\ + I2 are the same within the numerical precision 
when both leaks are simultaneously open (three columns 
on the right). In addition, it is interesting to note that 
all values are consistent with those obtained from the re- 
lationship k ~ -ln(l-/i c )/(t co ii) c with /4 ,J (/j), {tco\\)l )J 
(the values in the last two rows of these columns). 



It is shown in (Bunimovich and Yurchenko 2011 ) that 



it is possible to construct examples in which arbitrary 
small escape rates are achieved even in the presence of 
leaks of arbitrary large sizes. This surprising claim can 
be understood intuitively from the results of this sec- 
tion. Starting from a system with an arbitrarily small 
leak consider expanding the leak in such a way that the 
new leak also contains many images (and/or pre-images) 
of the original leak. This new leak can take an arbitrarily 
large proportion of the phase space without affecting the 
saddle and thus the escape dynamics which is by con- 
struction slow. 



C. Emission 
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FIG. 17 (Color online) Far-field emission for the Cardioid 
billiard with a leak centered at (si,pi) = (0, 0), with As = 0.5 
and Ap = 0.2. The far- field intensity distribution is computed 
by collecting the number of trajectories emitted under angle 
in the direction <3>, measured from the center of the billiard 
over the time interval (to . . . 00) with to ^> 1. The initial 
density was uniform in the full phase space po(x) = yO M (x), 
1000 bins in $ £ [— 7r,7r] were used. 

The most natural observable quantity in the configu- 
ration space of leaking systems is the emission of trajec- 
tories through the leak. Here we provide a representative 
configuration in which emission plays an important role, 
and connect the observed quantities to our theoretical 
formulation. 

Consider that detectors are placed around a circle far 
away from the cardioid billiard, as usually considered in 
optical microcavities. We introduce a leak on the right 
hand side of the Cardioid billiard with si = 0,pi = 0, 
As = 0.5, Ap = 0.2. The detectors collect the intensity 
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FIG. 18 (Color online) Emission from the Cardioid billiard 
described in Fig. |TT] over the full configuration space. The 
color code indicates the density of trajectories outside the 
billiard at large times t > 50. Fig[l7|is obtained by summing 
up all the intensities in a given direction. 10 8 initial conditions 
were used and data are presented on a grid of 200 x 200. 



of light emitted through the leak under different polar an- 
gles measured from the center of the billiard. Numer- 
ically, we distribute a large number of trajectories with 
an initial density po ( x ) m the phase space and record the 
emission angle of escaped trajectories in the configura- 
tion space. The proportion of trajectories escaping in a 
small interval around <3> is measured. For short times, the 
distributions are very irregular and depend strongly on 
the initial distribution po( x ) of the rays inside the cavity. 
After a transient period of the time, p(x, t) —> p c (x), the 
total intensity decays exponentially and the shape of the 
spatial distribution remains the same. Figure [TT] shows 
the spatial distribution as a function of <£, the so-called 
far-field emission. A multi-peaked non-uniform emission 
is observed. The spatial density of light rays for t ^> 1 is 
shown in Fig. [18) 




FIG. 19 (Color online) Invariant density p c (s,p) inside the 
leak of Figs. |TT| [18] The emission patterns are ultimately de- 
termined by p c (s,p), e.g., the zero density around (s,p) = 
(±0.5,0) leads to the zero emission around cj) — ±tt/2 in 

Fig. m 



For long times, the emission is fully described by the 
c-density p c (x) inside the leak (x G I). For instance, the 
emission angle $ is a geometrical function of x = (s,p) 
in the leak, $ = 0(x), and the far- field emission p(&) is 
proportional to 



p c (x)(5($ - <£(x))dx. 



(89) 



Figure 19 shows p c (x) in / for our example. The leak 
placed only on the right side of the billiard s G [—0.5, 0.5] 
and the non-uniform distribution of p c ( x ) are responsi- 
ble for the peaked and non-uniform emission in Figs. [TT] 
, The filamentary structure of p c (x) inside / seen 
19 reflects the filamentary pattern of the unstable 



and 18 



in Fig. 

manifold of the chaotic saddle, as discussed in Sec. II. B 
After performing the projection by p c (s,p) as indicated 
in Eq. (89), these filaments give rise to the zigzaged far- 



field emission shown in Fig. 17 In Sec. VI.D we shall see 
that similar emission properties can be experimentally 
observed in lasing microcavities, in which case the factor 
(1 — i?(x)) has also to be included in Eq. (89) in order to 
account for the partial reflection/transmission property 
of the leak (see also Sec. |II.P| ). 



V. EXTENSION TO WEAKLY CHAOTIC SYSTEMS 
A. Closed-system phase space 

So far we have focused on the case of fully-chaotic sys- 
tems and used the cardioid billiard to illustrate the the- 
ory. While these results apply to a broad class of (Hamil- 
tonian and dissipative) systems, there is an evident need 
to expand them to the larger class of weakly chaotic 
systems. For instance, in closed Hamiltonian systems 
stable periodic orbits and, around them, quasi-periodic 
Kolmogorov-Arnold-Moser (KAM) tori form regions of 
regular motion in the phase space. These regions co- 
exist with regions of chaotic motion. The situation is 



illustrated in Fig. 20 for the limacon billiard (20) with 
e = 0.46. Such mixed-phase- space systems are generic 
among all Hamiltonian systems and we focus on this type 
of nonlinearity. This is the typical case for billiards whose 
boundaries are defined by arbitrary (smooth) curves. Fa- 
mous examples are the limacon ( |Robnik , |1983 ) b illiard 
for any e 1 in Eq. ( 20 ) and the annular billiard ( |Saito| 



et al ,1982). 



The arguments above emphasize that the closed system 
is not ergodic: the phase-space is divided in multiple in- 
dependent components. One could imagine that the dy- 
namics inside any chaotic component would be described 
by the theory above. In reality, the situation is more in- 
volved because KAM tori are sticky surfaces that can be 
thought to affect trajectories in the surrounding chaotic 
component. Due to the smoothness of the dynamics, the 
local finite-time Lyapunov exponent close to the border 
of the KAM islands approaches zero (the Lyapunov ex- 
ponent of the tori). A trajectory in the chaotic region 
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(a) 



(b) 1 





FIG. 20 (Color online) Limacon billiard denned by Eq. (20) 
with e = 0.46. (a) Configuration space and (b) phase space of 
the closed billiard. The blue trajectory (•) is a stable period-4 
orbit around which a KAM island exists. A chaotic trajectory 
is shown as gray lines in (a) and as black dots in (b). The 



dashed line in (b) indicates the leak used in Figs. 21 and 22 



that come near some KAM surface, wanders close to that 
surface for a long time before leaving it and showing in- 
termittent bursts of chaos. This effect is called stickiness 
and the dynamics is said to be weakly chaotic (refer- 
ences are given in the next Section). Stickiness is typical 
in non-hyperbolic dynamical systems and can be thought 
as a consequence of the vanishing local Lyapunov expo- 
nent in the sticky region. The independence of this effect 
from non-ergodicity becomes evident by noting that even 
zero measure sets can lead to stickiness and weak chaos 
(e.g., the bouncing ball orbits discussed below). 

When a leak is introduced in a mixed-phase-space sys- 
tem, only the ergodic components that intersect the leak 
will be affected. Here we focus on leaks placed in chaotic 
components and, accordingly, the estimations of the mea- 
sure of the leak has to be normalized by the mea- 
sure of the chaotic sea /i(^ c haos)- This usually involves 
estimating the measure of the regions of regular motion 
M^reguiar) ano ^ subtracting from fi(ft) = 1 because typi- 
cally = /i(^chaos) + ^(^reguiar). Beyond this trivial 
correction accounting for the non-ergodicity of the sys- 
tem, weak chaos and stickiness is manifested in the sur- 
vival probability P(t). 



B. Decay of the survival probability in open systems 

The most important effect of the presence of sticky 
regions on the survival probability P(t) is that it modi- 
fies the asymptotic decay from exponential to power-law 
P(t) ~ t~ z , where z is the stickiness exponent related 
to the properties of chaotic regions close to regular ones. 



observables such as Poincare recurrences (Chirikov and 


Shepelyansky, 1984| |Zaslavsky[ 2002|) (as in Sec. II.CD, 


long-term correlations (|Chirikov and Shepelyansky 


1984 


1999 Karney 1983), 1/f spectrum JGeisel et al.\ 


1987D, 


Lyapunov exponents (|Artuso and Manchein 2009| 


Kantz 


and Grassberger 1987), and anomalous transport 


Geisel 


et al.\ 1988| Karney[ 1983). The value z = 2 can 


je ob- 



or marginally unstable (parabolic) periodic orb its (Alt 
mann et al. [[2008a Gaspard and Dorfman[|1995[), such as 



1995 



nette 



Kokshenev and Nemes[ 2000| |Legrand and Sor- 



the Sinai (Bauer and Bertsch, 1990j Fendrik and Sanchez 



1990a ), stadium ( Alt et a/.||l996||Armstead et al. 
2004f|Dettmann and GeorgiouH2Q09| |Nagler et a/.||2007 
1983[), the mus hroom (jAltmann et al. 



Vivaldi et 

2005} |Dettmann and Georgiou[ |2011a| |MiyaguchT 



Tanaka and Shudo 2006), and other billiards (Altmann 



et al, 2008a Fendrik and Wisniacki, 1997), and also 



2007 



for area preserving maps with sharply divided phase- 



space (Akaishi and Shudo 2009 Altmann et al. , 2006 



Fendrik and Wisniacki, 1997). Even a single marginally 



unsta ble point in area-preserving ma p can lead to stick- 
iness (lArtuso and Prampolini, 1998[), allowing for a di- 



rect connection to one-di mensional (Pomeau-M anneville 



type) intermittent maps (Artuso et al. |2008 ). Stick- 



iness also appears in higher dimensional Hamiltonian 



systems (|Altmann and Kantz , |2007 Ding et 



1990 



Fendrik and Sanchez, 1995). In area-preserving maps 
different stickiness s cenarios can be distinguished ( |Za-| 



slavsky, 2002 2005), but there are examples of bil- 
liards with divided phase space for which stickiness is 
absent (Bunimovich, 2008). For the generic KAM sce- 



nario strong fluctuations are observed due to the pres- 



ence of Cantori (Meiss 1992) acting as partial barrier 
to the transport of particles. The universality of z in 
the KAM scenario is an old problem that has not been 
completely solved despite different approa ches and sub- 
stantial advance in the last thirt y years (jChirikov and 
I Shepelyanskyl |l984l |l999l [2002] | Crist adoro and Ketzm- 
erick||2008||Meiss and Qttj [19851 |1986||Venegeroles||2009| 



Weiss et al.\ 2002^7 Th e most recent results ( |Cristadoro 
and Ketzmerick[ |2008 ) indicate the universal exponent 
to be z ~ 1.57. 

Even if for long times the decay is expected to be 
power-law due to the tori, there are interesting pre- 
asymptotic regimes of P(t) for systems with mixed phase 
space. For very short times t < t s: initial-condition and 
system-dependent fluctuations are typical and may play 



a predominant role in specific applications (Altmann and 
Tell [20081 IDietz et al.\ |2006a| |Grete and Markus[ |2007| 



Tanaka and Shudo, 2006[ ). Furthermore, when the leak 



is far away from any KAM tori and the chaotic compo- 
nent is large 12 , typical trajectories will exit before hav- 
ing the chance of approaching the KAM islands. These 
trajectories will experience an effective hyperbolic system 
and render an intermediate-time exponential decay, as re- 



ported in different systems ( Alt et al. , 


1996 Altmann and 


Tel[ |2008; Dettmann and Georgiou, 2009, 201 la| Fendrik 


et al. 


19941 |Gaspard and Dorfman| |1995| |Jung et al.\ 


1993 


Kokshenev and Nemes 2000 Zaslavsky, 2002 


)• 



tained analytically in leaking billiards with bouncing ball 



12 This situation is particularly important in high-dimensional sys- 
tems for which the measure of the regular regions decreases (but 
is still different from zero). 
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FIG. 21 (Color online) Survival probability P(t) inside the limacon billiard depicted in Fig. 20 with a leak / centered around 
si = 0.5, pi = with 2As = 1% of the perimeter. Two different initial ensembles po(x) are taken: p r (x) on the image of / 
(continuous, black); and p M (x) uniform in the chaotic component of the closed billiard (dashed, green). The red curve fits P(t) 
between times 1/k and 2/k. (a) Log-linear representation; Inset: magnification for short times, (b) Log- log represent atio n. 
P(t) decays with intermediate escape rate k ~ 0.0057(3), and with a power-law P(t) ~ t~ z for t > t cr oss ~ 1600 (see (90)). 
Thick blue lines with scalings z — 1.75 and Zt r — 0.75 are shown for comparison. 



Based on the different decay regimes di scussed above 
we can w rite the survival probability as ( |Altmann and 



Tel 



2008): 



P(t) 



irregular 



ae 



■b(Kt)- 



for < t < t s 

for t s <t<t z (90) 

for t z < t 



where t z corresponds to the time needed for the first tra- 
jectories to approach the sticky region, a/b is propor- 
tional to ratio of the measure of the chaotic and of the 
regular components of the phase space, and t 3 ,t z depend 
on the initial conditions po( x )- Even if trajectories 
are started in the sticky region (see next section), both 
power-law and exponential regimes are seen because for 
large chaotic components ae~ K ^> bn~ z . These regimes 
are clearly observable in the example of the limacon bil- 



liard, as Fig. [21] illustrates. 

An actual crossover time t C ross between th e exponential 
and the algebraic decay can be defined as ( Akaishi and 
Shudo| |2QQ9i |Altmann and Tel| |2008[ ) 



/^cross) 



(91) 



which is the time when the contributions from the hy- 
perbolic and the nonhyperbolic components are of equal 
importance. 

We can easily estimate the dependence of the crossover 
time on the size of the leak assuming that the ratio 
b/a depends at most wea kly on k for smal l From 
Eq. ^ we obtain that QAltmann and Tel[ |2008[ ) 



(92) 



A logarithmic correction to this relation was found 



2009) as t c 



by ( |Akaishi and Shudo 

(log(//(/)) + 1)M J )- See also. QAltmann et al\ |2012 



Dettmann and Georgiou J [20091 |2012| IKokshenev" ana 



Nemes 2000) for exact calculations and simulations in 



specific systems. Taking this correction into account, af- 
ter the usual re-scaling of time t \-> tn and in the limit 
— » we find that in the rescaled units k — >> 1 and 
Across —> i-e., the exponential decay always dominates 
P(t). This provides an example of precise statements for 
infinitely small leaks that mask an interesting dynamical 
phenomenon (the power-law) because for any finite leak 
a transition to power-law decay exists. Instead, here we 
do not apply any re-scaling and discuss the dependence 
of the intermediate time regimes in Eq. ( 90 ) on leak size 



and other parameters. Alternatively, one could take a 
different rescaling of time, e.g. t \-> tfi(I)/ log/i(J), that 
would not suppress the crossover. 



C. Dependence on the initial distribution 

The connection between escape time in open sys- 
tems and recurrence time in closed system described 
in Sec. II. C remains valid for weakly chaotic system as 
well. In particular, the intermediate-time decay of P(t) in 
Eq. (90) appears in the (cumulative) distribution P r (T) 
of the Poincare recurrence times T: P r (T) ~ P(t) t=T . 
With initial condition p r (x) [in the image of the leak / 
taken with respect to f c i OS ed as in Sec. II.C )] the entire dis- 
tributions coincide P r (T) = P(t)t=T- I n area-preserving 
cases p r is a uniform distribution over the image f c iosed (I) 
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of the leak. 




FIG. 22 (Color online) Chaotic saddles (a,c) and the corre- 
sponding unstable manifolds (b,d) characterizing the dynam- 
ics in the billiard depicted in Fig.|20]at different time instants 
t* . In (a,b) t* ~ 1/k « 175 < tcross ~ 1600 and the hyperbolic 
component of the saddle is dominant. In (c,d) t* ~ 2t croS s and 
the nonhyperbolic component is dominant. This is responsi- 
ble for the power law decay in (90). We considered initial 
conditions uniformly distributed on the leak, see Sec. |VIII.C| 
for details. 

The power-law exponent shows an important depen- 
dence on the initial conditions po( x )- In fact, z in 
Eq. (90) characterizes "scattering" cases, i.e., situations 
in which the support of po( x ) is far away from KAM 
surfaces. For po(x) inside the sticky region, arbitrarily 
close to KAM islands, the escape process - called tran- 
sient chaos situation ( |Pikovsky , 1992) - can be shown to 
be characterized by a survival probability P tr for which 



P tr (t)~t- 



for large £, 



(93) 



with a different decay exponent z tr . As shown in ( |Meiss| 
1997||Pikovsky[ [T992) and explained below, Zt r is smaller 
than z with unit difference: 



Ztr — Z — 1. 



(94) 



The impact of such a slower decay of the survival proba- 



bility can be clearly seen in Fig. 21 since the support of 
p r remains far away from KAM islands, but this is not 
the case for p M . 

The similarity between a properly opened up dynam- 
ics and the Poincare recurrences in the closed system can 
be used to explain the difference between the algebraic- 



decay exponents z and z r (Altmann et al, 2006). Con 



sider first initial conditions touching the sticky region, 
the case of transient chaos, and examine the time a 
trajectory takes to escape to a region far away from 



the sticky region. The survival probability distribution 
Pt r (r) is proportional to the natural measure /x(r) of the 
region of the phase space to which the trajectories stick 
for a time longer than r. Due to ergodicity, we can write 



P tr (r) ~ h(t) = 



^total 



(95) 



where t T is the total time spent in the sticky regions (in 
events with recurrences time T > r) within the total 
observation time ttotai? see Fig. [8|. 

For the recurrence problem with a single trajectory of 
length ttotai initialized far away from the sticky region, 
the cumulative probability P r (r) to find recurrence times 
T larger than r can be expressed as 



P r (r) 



Nr 

N 



P(r), 



(96) 



where N T is the number of recurrences with recurrence 
times larger than r and TV is the total number of recur- 
rences observed in the time interval t t otai- The right- 
hand side expresses the connection between recurrence 
times and escapes, Eq. (43). Since the total observation 



time can be estimated as N times the mean recurrence 
time f [Eq. [4l|], we have t total ~ NT. Similarly, the 
total time t T spent inside the sticky region is approxi- 
mately the number N T of recurrences with times longer 
than r multiplied by r: t T ~ N T r. As T is a constant, 
independent of r, these allow us to write 



P(r) 



N 



tr/l 



t T T P tr (r) 



(97) 



from which the shift of algebraic decay exponent by one 
immediately follows from the asymptotics of Eqs. 
and (1931). 



D. Hyperbolic and non-hyperbolic components of chaotic 
saddles 

A theory based on invariant sets has to consider that in 
this case the chaotic saddle responsible for transient cases 
encircle KAM tori, moreover, it comes arbitrarily close to 
the to ri ( |Christiansen and Grassberger[ 1993 Lau et al , 
|1991|) . The survival probability in Eq. (1901) suggests that 



the invariant sets governing the temporal decay of weakly 
chaotic systems can be divided in hyperbolic and non- 



hyper bolic components (Fendrik et al. t 1994 Jung et al 
1993 ) . The intermediate-time (t s < t < t z ) exponential 



decay is the manifestation of the hyperbolic component 
of the chaotic saddle, and the power-low decay becoming 
observable for long times t > t cross of the nonhyperbolic 
component. 



In Fig. 22 we show numerical approximations of the 
chaotic saddle and its unstable manifold in the hyper- 
bolic and non-hyperbolic regimes. We have employed the 
same procedures used before for the case of fully chaotic 
systems, see Appendix |VIII.C[ but varied the effective 
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time t* used in the simulation. For intermediate values 
of t* in which the exponential decay dominates P(t) in 



Eq. (90), the obtained saddle and its unstable manifolds 
are disjoint from the region containing the KAM island 
and show the structure typical of hyperbolic systems. For 
large times t* > t cross the asymptotic dynamics is gov- 
erned by the non- hyperbolic regions, and the saddle and 
unstable manifold concentrate around the islands. 



We assumed, as in Sec. |II.B[ that a single chaotic sad- 
dle exists. An interesting and nontrivial example that 
violates this case has been reported by Dettmann and 
Georgiou in the stadium billiard with two full leaks 1{ at 
the boundary: I\ in the flat and I2 in the circular com- 
ponent of th e billiard's boundary (Dettmann and Geor- 
|giou[ |2011b[ ). Wit h onl y one of these leaks, P(t) ~ t~ z 
as reviewed in Sec. V.B Consider now the survival prob- 



ability P- (t) of particles that start in leak i and leave 
through leak j (i = j corresponds to a reflection survival 
probability and i =^ j to a transmission survival proba- 
bility). From the four P-(t) in the example above, only 
Pi(t) shows the expected power-law tail. P x 2 (t), P 2 1 (t), 
and P| (t) decay exponentially. This surprising result is a 
consequence of the sticky region due to the parallel walls 
in the Stadium billiard, which cannot be approached from 
all positions of the phase space if a leak is placed along 
these parallel walls. More generally, this example shows 
that the ergodicity of the closed system is not automati- 
cally transferred to the leaking system. 

A curious effect tha t is beyond the result s of th is pa- 
per was reported by (Custodio and Beims, 2011): the 



shape of the leak (rounded, squared, etc.) introduced at 
the border of the billiard can modify the dynamics of the 
billiard (e.g., create sticky or chaotic motion). Similarly, 
openness in optical systems can lead to a modification of 
the dynamics of the reflected rays due to non-specular re- 
flection close to the angle of total internal reflection J Alt- 



mann et al. 2008b Schomerus and Hentschel 2006D. In 



this paper we only consider leaks for which the trajectory 
either escapes or is reflected as in the closed system. 

Finally it is worth comparing the results of this Sec- 
tion to other types of billiards with leaks. In integrable 
cases, leaks (or recurrence regions) in billiards (|Bauer and| 



Bertsch, 1990; Bunimovich and Dettmann, 2005; 


Vicen- 


tini and Kokshenev, 2001 


) and maps (Buric et al. 


, 2003) 



due to stickiness caused by families of marginally unsta- 
ble periodic orbits (e.g., parallel walls) and is consistent 



with Eq. (94) since there is no chaotic region and hence 
all initial conditions are close to sticky regions. Leaks 
placed inside regular islands will lead to similar observa- 
tions. Leaks (or recurrence regions) centered in periodic 
orbits have been considered in Ref. ( Hu et al\ 2004). 
For leaks centered at the border of regular and chaotic 
regions, a weighted sum o f exponential and power law 
escapes has been found by ( Buric et al. , 2003). A similar 
composition was also reported in (|Kokshenev and Vicen- 
[tinil|2003|[Vicentini and Kokshenev||200lD for the case of 
leaking polygonal billiards, which are neither integrable 



nor chaotic. (Dettmann and Leonel, 2012| |Leonel and 



Dettmann, 2012) report a decay of survival probability 
with a stretched exponential P(t) ~ exp(— at@) when a 
leak is introduced in the open bouncer model (vibrating 
billiard), which exhibits mixed phase space and Fermi 
acceleration in the closed version. 



VI. APPLICATIONS 

We now present different applications for which the 
results of the previous sections are relevant for the un- 
derstanding of specific phenomena and observations. Ac- 
cordingly, the models appearing below intend to capture 
the essential dynamical aspects of the system and not to 
describe all details of a given experimental configuration. 
The extend into which dynamical systems provide an ap- 
propriate description of reality has been subject of dis- 
cussion in the context of wave dynamics for a long time. 
For instance, in acoustics the dynamical systems descrip- 
tion corresponds to the ray approach, which has been re- 
peatedly debated ( | Joyce , 1975) but systematically used 
nevertheless. 



Two recent advances that confirm the va- 
lidity of the ray picture are: (i) the derivation of Sabine's 
law in the wave picture (Dennis, 2010 Legrand and Sor- 



nette 



1991b) and (ii) the recognition of the relevance of 



2010) 



classical periodic orbits in room acoustics (Berry 
Similar considerations involving ray and wave pictures 
apply to the quantum- and optics- applications below, 
while the astronomy-, fluids-, and plasma- applications 
have also simplifying assumptions and regimes of valid- 
ity of their own. 



A. Planetary astronomy 

Chaotic systems with leaks might play a role in one of 
the most traditional problem of celestial mechanics, the 
three body problem (three masses interacting through 
gra vitational forces). This has been noted by Bleher et 
al. ( |Bleher et al. |1988[ ) and exp l ored in detailed fo r par- 
ticular cases by Nagler flNaglerl 120021 120041 120051). The 



most natural leakage mechanism corresponds to collisions 
of the finite-size cel estial bodies, a real possibi lity even 
in our solar system (Laskar and Gastineau 2009). 



The simplest case discussed in Ref. ( Nagler 2004 ) cor- 
responds to two main bodies with equal mass M moving 
along the same circle around their center of mass, and a 
test body of mass m <C M (circular restricted three-body 
problem). In the co-rotating reference frame centered at 
the center of mass, this problem can be described by a 
four-dimensional time-independent system. By using the 
conservation of energy and employing a suitable Poincare 
surface of section, one can reduce the dynamics to a two- 
dimensional discrete-time system, as the billiard systems 
considered in this paper. For a fixed energy, confined tra- 
jectories coexist with trajectories that escape to infinity. 
The system is non-integrable, with regular and chaotic 
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FIG. 23 (Color online) The configuration space of the small body in the restricted three body problem in a frame co-rotating 
with the main bodies. This section of the phase space is obtained at a fixed energy under the condition that the particle's 
velocity toward the origin vanishes, and the angular velocity is negative. The colors correspond to trajectories colliding with 
main body 1 (white), colliding with main body 2 (gray, red), and escape to infinity (dark - black to light blue, from short to 
long escape times). From left to right the diameter of the two main bodies is increased, wh at correspond s to the size of the 
leaks (not visible in this section) and to an increase of the probability of the collisions, from ( |Nagler[ |2004| ). 



trajectories coexisting in the phase space. 

Leakage is introduced by considering that the masses 
of the two main bodies are not concentrated in a point 
and therefore collisions occur whenever the point-like test 
particle approaches one of the bodies to a distance r < R, 
where R is the radius of the main bodies (the velocity 
can be arbitrary). The size of the leak in the configu- 
ration space is proportional to the radius R of the main 
bodies, and the position of the leak is given by the posi- 
tions of the main bodies, which is fixed in the co-rotating 
frame. This situation corresponds to the problem of sys- 
tems with more than one coexisting leaks discussed in 
Sec. |IV.B| In fact, at a given energy, particle can also 
go to infinity, therefore there are here two leaks in an 



otherwise open system present. Figure 23 illustrates how 
the basin of escape and the basin of collision with each 
of the main bodies vary as a function of the leak size R. 
The typical fractal-like structures observed in Fig. [16] are 
clearly seen, and coexist with smooth boundaries ( |Na-| 
gler, 2QQ4[ ). ( |Nagler| [2005 ) extended this analysis to the 
case of different masses and different circular orbits for 
the main bodies. 



B. Hydrodynamical flows 

The advection of tracer particles in hydrodynamical 
flows represents an important application of dynamical 
system theory. Since molecular diffusion is typically neg- 
ligible on the relevant time and length scales, the equa- 
tion of motion for an idealized particle of zero size and 
zero mass expresses the fact that the particle velocity, r, 
coincides, at any instant of time, with the flow velocity 
u(r, t). The velocity field is assumed to be known and 
the advective dynamics is thus described by 



r(t)=u[r(t),t]. 



(98) 



The solution to this differential equation is the path r(i) 
of the particle. 

Chaos is typical in two-dimensional time-dependent 
flows and in any kind of three-dimensional flows, lead- 
ing to the phenomenon of chaotic advection ( |Aref[ |1984[ ) . 
Indeed, the main physical mechanism for fluid stirring is 



advection, whose efficiency can greatly be enhanced by 
chaotic dynamics. The spreading of pollutants on large 
scales is also dominated by advection. Potential appli- 
cations of chaotic advection thus range from laboratory 
investigations of fluid dynamics to the study of large-scale 
environmental flows, and these aspects are well reviewed 



in the literature (Aref and Naschie 1995 Lai and Tel 
20lT| lOttinol 119891). 



The existence of a leak in a flow implies a sink for 
fluid elements. In such cases the amount of fluid in a 
finite container is decreasing in time, and this leads to a 
qualitative change of the dynamics sooner or later. Here 
we focus on a less evident, but physically more appeal- 
ing, realization of leaking dynamics, which is related to 
advected particles and to the first arrival to certain re- 
gions of the flow. The problem of reactions in fluid flows 
provides an interesting example. Reactive particles have 
typically no considerable influence on flow, therefore the 
velocity field remains the same u(r, t) as without reac- 
tions. The reactive dynamics, represented, e.g., by a 
temporal change of certain particle properties, is thus su- 
perimposed on the advection problem (98). The change 



might happen upon entering a region of the flow. The 
most important examples of activities in flows involve 
chemical (e.g., A + B — > C + D) and biological (e.g., 
A + B — >• 2A) reactions. The combination of reaction 
and advectio n provides a re alistic model for a pletho r a of 
applications (Neufeld and Hernandez-Garcia 2009 Tel 



et al. 2005). Below we discuss two particular problems 



in which leaking is introduced in a closed hydrodynam- 
ical flow by allowing particles to enter or also react in 
pre-assigned fluid regions (that play the role of leaks) . 



1. Spreading of pollutants in the environment 

Imagine that a pollutant is released in an observation 
region within a water basin. A typical problem in the 
prevention of environmental pollution is to determine 
which coastal region the pollutant will be advected to 
so that one can estimate which parts are most likely to 
be affected by the pollution release. This is one par- 
ticular example of the general problem of partitioning 
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the initial conditions in an observation region of a closed 
flow according to the first arrival of advected particles 
to pre-determined sub-regions of interest. The boundary 
between the different partitions typically shows fractal 
patterns that correspond to the stable manifold of the 
chaotic saddle formed by tracers that never reach any of 
the target regions. The target regions act as leaks for the 
dynamics of the tracers [following Eq. (98)], but not for 
the velocity field u(r,t). 




FIG. 24 Spreading of pollutants in a model of a wind-driven 
lake. The wind field is periodic with period T. Positions 
x, y are colored according to which of the four coastal regions 
(bands along the boundaries, also colored) a tracer starting 
with initial condition x, y [at time t = (mod T)] is first 
advected to. Tracers that do not reach any of the bands along 
the shores over 60 hours of observation are colored white. 
Panels a)-d) represent layers at depth 0.7, 1.0, 1.3, and 1.6 m. 
(Picture by S. Szanyi.) 



Figure [24] shows an example in a square-shaped wind- 
driven lake of 2km by 2km. The boundaries are verti- 
cal walls of height 2m, and below this depth, the lake 
has a pyramidal form with the deepest point at a depth 
of 2.5m in the middle. The water flow is generated by 
wind stress. A similar problem was studied in (Karolyi 



et al. , 2010) by applying a shallow- water approximation 



in which layers of different depths are assumed to move 
in a synchronized manner. In Fig.[24]the flow is obtained 
from a numerical solution of the hydro dynamical equa- 



tions in three-dimensions (Cioffi et al. 2005| ). A wind of 



strength 12m/s changes periodically in time, with a pe- 
riod of T = 8 hours. It blows from south-west (left lower 
corner) for T/2, changes abruptly to south-east, and after 
an interval of length T/2 it changes back again, etc. 

The region of observation is the full lake area outside 
narrow bands along the costs. Layers at different depths 
are investigated. Each point in a layer is colored accord- 



ing to which of the coastal bands (of width 100m and 
of height 2m) along the four vertical walls at the shores 
will be reached by the tracer first (the vertical coordi- 
nate of the first arrival is not recorded). The results in 
Fig. [24] indicate a strong height-dependence. In the up- 
permost layer, pollutants released in the rectangl e ar e 
most dangerous for the northern shore (Szanyi] 2012). 
The western shore is somewhat less polluted than the 
eastern one. The southern coast is hardly affected. This 
is consistent with the fact that the average wind direction 
is southerly, and particles are thus pushed mainly north- 
ward. Considerable deviations from this pattern occur 
when going to deeper and deeper layers. The pollution 
of the northern coast is becoming less and less strong, 
while the southern one becomes heavily polluted. This 
is due to the development of an overturning circulation 
that has an overall southward component in the deeper 
regions. The hazard for the four shores thus strongly de- 
pends on the level in which pollution is released. The 
boundaries between different colors contain the stable 
manifold of a three-dimensional chaotic saddle residing 
in the region of observation. The main result obtained 
from this example is that a two-dimensional shallow wa- 
ter approximation might lead to an oversimplification of 
the three-dimensional flow and also of the advection pat- 
terns since the observed height-dependence in the pollu- 
tant distribution cannot appear then. 



2. Reactivity in flows, resetting 

Perhaps the simplest model of reactions is provide d by 
what is invented and called by R. Pierrehumber t ( |Ngan| 



and Pierrehumbert, 2000| Pierrehumbert, 1994) as the 



resetting mechanism. Whenever a tracer enters a prese- 
lected region of the flow, a given property of the tracer 
such as concentration (or color) is reset to a value asso- 
ciated with tha t regio n, regardless of its previous value 
( jNeufeld et al. [2000). This mimics a situation where 
dye is introduced by diffusion from a solid surface and is 
maintained at the saturation concentration in a diffusive 
boundary layer. An atmospheric example of the reset- 
ting mechanism is provided by the dynamics of water 
vapor. Water vapor is removed from an air parcel when- 
ever it enters a region where the local humidity is lower 
than that of the parcel. Its humidity is then reset to this 
lower value and the difference is rained out. When a par- 
cel of low humidity comes into a humid region (typically 
close to the Earth's surface), its water vapor content is 
reset to a high value. The water vapor distribution in the 
atmosphere is similar indeed to the ones obtained from 
resetting models (Pierrehumbert et al., 2007). As in the 



example in the previous Section, the leaking mechanism 
always refers to the particles of a given type (color), and 
not to the fluid, as fluid is never lost from the system. 

In two-dimensional incompressible flows Vu=0, and 
the advective dynamics (98) is area-preserving. The re- 



setting problem provides thus a close analog of a leaky 
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FIG. 25 (Color online) Tracer distributions with resetting in a closed time-periodic flow. The particles' color (concentration) is 
set to blue (dark) and yellow (light) at two resetting bands with total area 0.32 placed (a,d) vertically, (b,e) horizontally, and 
(c,f) as a triangle at the corners. In the full phase space the two colors come arbitrarily close to each other along a filamentary 
pattern. In panels (a)-(c) the bands foliate a fully chaotic flow, while in (d)-(f) the fluid flow is slower and KAM tori appear 
outside the resetting regions. (Data by G. Drotos). 



2D billiard. Chaotic advection in closed flows is char- 
acterized by space-filling chaos, and correspondingly by 
space-filling stable and unstable manifolds. Nonhyper- 
bolic regions might also exists around KAM tori. The 
resetting mechanism reveals the foliations of the dynam- 
ics, in a similar spirit as a leak, because in both cases 
a particle can be be considered to be lost after entering 
a preselected region. The advective dynamics with the 
leak is typically transiently chaotic. Particles never es- 
caping the complement of the leak(s), both forward and 
backward in time, form a chaotic saddle. The manifolds 
of this saddle are subsets of the closed system's manifolds 
since the advective dynamics outside the leak is exactly 
the same as in the closed system. 

With two or more resetting regions, tracers of differ- 
ent concentrations or different colors come close to each 
other along fractal-like boundaries. An example is shown 



in Fig. 25 where the flow is chosen as the alternating si- 
nusoidal shear flow model ( |Pierrehumbert| |1994[ ). The 
dark (light) dots are obtained as initial conditions of tra- 
jectories reaching the dark (light) region when iterated 
backward in time. Particles close to the boundary have 
long lifetimes outside the resetting regions and therefore 
they come close to the saddle in the time-reversed dy- 
namics. The boundaries trace thus out here the unstable 
manifold of the chaotic saddle. The resetting pattern is 
thus a fingerprint of the chaotic saddle underlying the 
leaking advection dynamics, where the leak is the union 
of the resetting regions. Figure [25] shows that resetting 
regions of the same area lead to different patterns and 
also to different average lifetimes (panels a)-c)) in full 
agreement with what we saw for fully chaotic billiards in 
Sec. IV Panels d)-e) represent a case where the average 



fluid velocity is smaller and four elliptic islands appear, 
with a nonhyperbolic component of the chaotic saddle 
around the outermost KAM tori. 



C. Magnetic field lines in plasma confinement machines 

The research on plasma confinement devices in the last 
sixty years has been driven by the expectation of de- 
signing controlled fusion reactors. One of the best stud- 
ied and most promising strategy is to use magnetic con- 
finement machines with toroidal shape, as illustrated in 
Fig. [26| These machines are called tokamaks ( |Wesson[ 
1987). In lowest order approximation the charged parti- 
cles in a tokamak follow its magnetic field lines. These 
lines are divergence free. They thus correspond to trajec- 
tories of volume preserving flows, and slices of them can 
be described by area-preserving maps ( Morrison] 2000). 



A chaotic layer at the border of the plasma column is 
long known to enhance the confinement of particles to 
the core. This enhancement is accomplished by control- 
ling plasma- wall interactions in the tokamak ( Engelhardt 



and Feneberg 1978[). Chaos is present whenever the ideal 



toroidal symmetry is broken which can be achieved by 
generating weak electric currents along the toroidal ves- 
sel via the so-called ergodic magnetic limiters shown in 
Fig. [26] (left). 

Here we are interested in the dynamics of the magnetic 
field lines which, in contrast to the complicated dynamics 
of the particles, allow for a simple description in terms of 
a near-integrable Hamiltonian systems. Figure [26] illus- 
trates the steps that map the magnetic field line dynam- 
ics into a discrete 13 map between two successive inter- 
sections with the Poincare surface of section with fixed 



toroidal angle <£, see (Portela et al t 2008) for a recent 
review. The connection to systems with leaks is estab- 
lished by noting that the magnetic field lines (and the 
evolution of the corresponding Hamiltonian trajectories) 
are interrupted at the vessel's wall and at any obstacle 
inside the tokamak chamber, like probes and antennas. 
The main loss of particles is due to collisions with the 



Resetting- like methods have been applied to vi sualiz 
ing the foliations of three-dimensional model flows ( |Tuval 
\et a/.[|2004]) and they also hav e applications in geophysics 
dSchneider et a/.l 120001 120071. 



13 The length of the field line along the toroidal angle <E> plays the 
role of continuous time in this description of the (stationary) 
magnetic field configuration. 
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FIG. 26 Construction of an area-preserving map of the magnetic field lines in a Tokamak. From left to right: main geometrical 
parameters of a Tokamak with an ergodic magnetic limiter; torus in which the magnetic field lines are confined and the surface 
of section at fixed toroidal angle $ used to construct the discrete-time map; different trajectories of the map shown in polar 
coordinates; and the same trajectories of the map shown in planar coordinates. 




FIG. 27 The saddle (blue), stable (black), and unstable (red) 
manifolds in the Ullmann-Caldas tokamak map. The (stable) 
unstable manifold was obtained as the 80-th (pre-) image of a 
small ball of initial conditions around the unstable periodic 
point indicated as a square at ~ (0.49, 0.78). The intersection 



of both manifolds approximates the chaotic saddle (Portela 
et al 



2007J). 



the heat flux is not unifor m in the poloidal angle (pro 



portional to x in Fig. 26) (Shen et al. 1989 Takamura 



et al.\ |1989| ). Apart from the disp ersion caused by colli- 



sional effects (Schelin et al. 



leave the system following c 



2011b, the charged particles 



osely the unstable manifold 



of the chaotic saddle and are deposited in the positions 
in which this manifold intersects the wall. This argumen- 
tation predicts that the heat flux should be proportional 
to the length (or escape time) of the magnetic field lines. 
Lines approaching the stable and unstable manifolds have 
increasingly large lengths. Figure [28| shows experimental 
and numerical results that confirm this prediction. In the 
upper panel the concentration of the heat flux in specific 
poloidal angles is shown for different edge safety factors 
q a , a control parameter that is inversely proportional to 
the plasma current. In the lower panel we observe that 
the same patterns are observed in the length of the field 
lines. 

Here, again, similarly as in the resetting problem, we 
have a chance to see the direct fingerprint of transient 
chaos in the configuration space since in both cases the 
phase space and the configuration space coincide. 



chamber's wall which act as a leak and is located in the 
large-?/ region in the rightmost panel of Fig. [26] As we 
saw, the leakage strongly affects the dynamics of the area- 
preserving map. 

The asymptotic dynamics of the leaking system is gov- 
erned by the chaotic saddle and its manifolds, which de- 
pend sensitively on the choice of the leak and are thus 
different from the invariant sets of the naturally open 



system. Figure 27 shows an approximation of these in- 



variant sets for the Ullmann-Caldas map (Ullmann and 



Caldas, 2000), based on the manifolds of a single fixed 



point embedded into the chaotic region ( jPortela et al 
|2007| ). Viana et al. provide a recent detailed review of 
the role of the chaotic saddle and fractal structures in 



plasma confinement devices (IViana et al. 2011). We em- 



phasize that the chaotic saddle is non-hyperbolic due to 
the presence of KAM islands, as discussed in Sec. [V) 

Experimental signatures of the invariant sets of the 
leaking dynamics appear in the heat flux into the toka- 



mak's wall ( Evans et al.||2002||da Silva et a;.[|2QQ2l|Viana| 



et al. ,2011 



Wingen et a/.||2007 ). It has been noticed that 



D. Optical microcavities 

Optical microcavities are used in applications ranging 
from dynamic filters in optica l com munications to quan- 
tum electrodynamics (see (Vahala, 2003) for a review). 
Here we focus on the case of optical microcavities used 
as lasers. Beside the practical applications, these cavi- 
ties allow for fundamental scientific investigations (e.g., 
in quantum chaos) in which systems partial leaks play 
a prominent role. An example of an optical microcavity 
is given in Fig. 29 taken from (iShinohara et al. 2010 



201l[ ), a system used here to illustrate t he achievements 



of the last 15 years of intense research (|Altmann 


2009 


Dettmann et al.\ |2009| |Harayama and Shinohara 


2011 


Lee et al.\ 2004 


Nockel and Stone| 1997 Ryu et al.\ 


2006 


Schwefel et al. 


2004; Tanaka et al.\ |2007| |Wiersig and 


Hentschel 2008 Yan et al, 2009). 



The connection to leaking chaotic systems is based on 
the observation that the shape of the microlaser mat- 
ters: while the regular trajectories present in spherical 
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FIG. 28 (Color online) Comparison between experimental 
and numerical results for the heat flux pattern in tokamaks. 
Upper panel: experimental measurements of the heat flux 
(color scale) through different poloidal coordinates (propor- 
tional to x in Fig. 26 ) shown for different edge safety factors 
q a . Lower panel: numerical simulations indicating the length 



of the field lines that hit the tokamak wall vs. q a . From ( Win- 
gen et al.\ [20071. 




FIG. 29 Scanning electron microscope image of a microlaser 
with deformed-disk shape. From (Shi nohara et al.\ 2010). 



and circular shaped cavities provide the good confine- 
ment necessary for lasing (high Q modes), deformation 
of these geometries leading to chaotic dynamics (|Mekis 



et al. 19951 iNockel and Stonel 119971) can add to the good 



confinement other desired properties such as the direc- 



tionality of the emission ( Dettmann et al. 


2009 Gmachl 


et al.\ 1998| Lebental et al.\ 2007| Liu and Lai[ |2002 


Nockel and Stone] 1997| Schwefel et al. 


[20C 
1 He 


)4j|Song et al. 


20091 Tanaka et al.\ |2007| Wiersig anc 


ntschel[ |2008 
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FIG. 30 Ray model of the microlaser depicted in Fig. 29 (a) 
Definitions of the coordinates, (b) A stable periodic orbit, 
(c) Phase space of the closed billiard. The upper edge of the 
partial leak is indicated by th e dashed (red) line in panel (c). 
From (Shinohara et al. 2011). 
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FIG. 31 Far-field emission observed experimentally (solid 
line) and predicted by the ray-model (dashed line). 



From (Shinohara et al 2011) 



|Yan et aT , 2009). Even if these systems are often far 
from the limit of vanishing wave-length (geometrical op- 
tics), the so-called ray- wave correspondence applies in a 
surprising large n umber of cases, see Ref. ( |Harayama and 
jShinohara |2011 ) for a recent review. A description based 
on rays makes robust numerical simulations possible and 
provides intuitive interpretations of experimental obser- 
vations that can be directly connected to orbits and the 
geometry of the corresponding billiard. 

The ray-model in optical cavities is equivalent to the 
dynamics in closed billiards with partial leaks. Figure 30 
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shows the phase space of t he system used in Refs. ( Shi- 
nohara et al\ 120101 |2011|). As mentioned in Sees. II. E 



and II. D[ the partial leak has a natural physical origin 
in the collisions inside the cavity (with refractive index 
ni n > n out = 1). Rays with angles smaller than the crit- 
ical angle of total internal reflection [p < p c = sin# c = 
1/riin, dashed line in Fig. |30fc) are partially transmitted 
(with angles given by Snell's law) and partially reflected. 
The intensities of the re flected and tran s mitted rays are 
given by Fresnel's law ( Altmann] |2QQ9 Harayama and 
Shinohara|[2QTT||Lee et a/.[|2004HNockel and Stone[|1997| 



Ryu at 2006; Sc hwefel et al\ |2004| |Tanaka et al 
2007). For rays with p < p c no transmission takes place. 



Interestingly, limacon shaped microlasers as the ones 
used in the main part of this pape r have also been used 
in experiments a nd simulations (|Shinohara et al , 2009 



Song et al. | [20091 1 Wiersig and Hentschel||2008||Yan et al 
2009). An important additional aspect of lasing cavi- 
ties is the gain medium that continuously pumps energy 
into the system enhancing the intensity of long-living tra- 
jectories (a precise modeling of this effect is beyond the 
scope of this paper) . For corrections to the ray model due 
to wave effects see Refs. (Altmann et al 2008b Schome- 



rus and HentschelH2006l ) . 

The main relevant observable is the far-field emission 
intensity, as computed in Fig. 17 (see also Sec. IV. C). 
A comparison between experimental results and ray- 



dynamics simulations is presented in Fig. 31 The agree- 
ment is remarkable in particular if one takes into account 
that the system is far from the formal semi- classical limit 
of small wavelength, see (Redding et al 2012 Shinohara 



et al, 20101 2011) for further evidence of the robustness 



of ray results. The important and somehow surprising 
aspect of the emission in Fig. |3l] is the concentration of 
the emission in specific directions, even if the leak is in- 
troduced in a region where the closed system is more or 
less uniformly chaotic. 

The specific directions of emission have been explained 
using the unstable manifold of a single unstable periodic 
orbits close to the border of the leak region, first shown 
by (|Schwefel et al 2004). Another important theoreti- 



cal development was the proposal of the existence of a so- 
called quasi-statio nary energy distr i bution in the clas sical 
ray dynamics, see ( Lee et al , 2004 Ryu et al 2006). As 
noticed in Ref. (Altmannl 2009), and better explained in 



this paper, these two important concepts can more gen- 
erally be expressed, respectively, as: the unstable mani- 
fold of the chaotic saddle (which aligns with the unstable 
manifold of specific periodic orbits) and the c-measure 
distributed along it. In as much the ray- model provides 
a good description of the laser, these two key concepts of 
the theory of dynamical systems have direct observable 
consequences for the far-field emission of lasing micro- 
cavities. Here, the results of Sec. [V] have to be taken 
into account bec ause microlasing cavities are typically 
not fully chaotic ( |Liu and Lai , |2QQ2| ) . The main effect of 
weak chaos is that the long-living modes (high-Q factors) 
are usually concentrated inside KAM islands and regions 



of regular motion. It remains to be shown in which ex- 
tent the effect of weak chaos described in Sec. [V|can be 
detected experimentally. 

Recent developments show links to many other results 
discussed in our paper. For instance, while dielectric mi- 
crocavities have leaks restricted in the p direction, recent 
experiments use waveguides that restrict the leak also 
in the other direction. Leaks are constructed around 
a position s\ that showed t o be efficient in channeling 
the rays out of the cavity (Redding et al 2012 Song 



et al 



2012). Such coexistence of different leaks was dis- 



cussed in Sec. IV. B Another example is the periodic- 
orbit formula for the resonance spectrum of dielectric 
cavities proposed by (Bogomolny et al\ 2008) and ex- 



perimentally veriff^oM^ and 



chaotic ( jBittner et al 2012[ ) systems. The main novelty 
in this formula is the incorporation of partial leakage as 
a multiplicative term proportional to the product of the 
reflection coefficients at collisions along periodic orbits, 
as the reflection coefficient appears in the operator pro- 
posed in Eq. ([78). 



E. Quantum and wave chaos in systems with leaks 

The quantum and wave analogues of classically chaotic 
systems with leaks appear in theoretical analyses, con- 
trolled experiments, and real- world systems. Quan- 
tum mechanically, the major differences to closed sys- 
tems originate from the non-Hermitian Hamilton oper- 
ator (non-unitarian Scattering matrix). The intrinsic 
openness of any quantum mechanical experiment, high- 
lighted in a quote from the early 1990's in Sec. \\\ is al- 



11 is a] 
" [1999 



ready fully discussed in textbooks ( Stockmann| 11999). 
Still, an important question influencing the research in 
quantum and wave chaos in the last decade was the ef- 
fect of openness on results known for closed systems. 
A driving force are experiments in quantum dots, mi- 
crowaves, optics, acoustics etc. These more recent devel- 
opmen ts appear in recen t essays on open quantum sys- 
tems: (Kuhl et al, 2005| ) reviews spectral and scattering 
properties with focus on classica l wave s (microwaves and 
sound waves); ( Fyodorov et al\ 2005) uses a formalism 
based on correlation and d istribution fun ctions that goes 
beyond random matrices; ([Rotter 2009 ) uses a Feshbach 



projection operator fo rmalism to non- Hermitian Hamil- 
tonian operators; and ( | Nonnenmacher 2011[ ) reviews the 
mathematical results and methods of semiclassical theo- 
ries of wave operators in scattering systems. 

The results of these review papers are in a great extent 
valid for scattering systems in general, with the main dis- 
tinction being between weak and strong absorption (cou- 
pling to the environment). Throughout this paper we 
have distinguished leaking from genuinely open systems 
based on the key elements: (i) the possibility of com- 
paring the results to a closed-system and (ii) the con- 
trol over properties of the leak (see Sec. LA). These two 



elements affect spectral properties, which have been re- 
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cently observed in microwave experiments ( Barthelemy| 
et al. , 2QQ5| Dietz et aT\ 2QQ6b|) and fully an alyzed theo- 



retically (Poli et al., 2009; Savin et al. 2006). Point (i) is 



used, for instance, when comparing the deviations from 
the universal distribution of nearest neighbor levels ob - 
served in fully chaotic closed systems ( Stockmann| 1999). 
Point (ii) is related to the importance of the localization 
of the leak or perturbation in the phase space, a recent 
observation that is nicely summarized by Savin, Legrand, 
and Mortessagne as: 



Open wave-chaotic systems in the presence of 
energy losses (absorption) are nowadays un- 
der intense experimental and theoretical in- 
vestigations . . . Most of the works concern 
the case of uniform absorption which is re- 
sponsible for homogeneous broadening T/ lom 
of all the modes (resonance states). How- 
ever, in some experimentally relevant situ- 
ations . . . one should take into account also 
localized-in- space losses which lead to an in- 
homogeneous part Tinh of the widths which 
varies from mode to mode. ( |Savin et al 



2006) 



The use of antennas and measurement devices typically 
fulfills properties (i) and (ii) mentioned above. 

Below we focus on three measurable properties of quan- 
tum systems which can be directly connected to classical 
properties and the main results of this paper. 



1. Loschmidt echo (fidelity) 

The absence of well-defined trajectories in quan- 
tum mechanics make it difficult to precisely define 



chaos (Berry, 1987), a concept usually defined based on 
the exponential instability of classical trajectories un- 
der small perturbations of initial conditions (one posi- 
tive Lyapunov exponent). The effects of classical chaos 
in quantum and wave systems can be observed following 
an idea by Peres, who proposed to compare the evolu- 
tion of a wave packet with its evolution in a perturbed 
Hamiltonian (Peres, 1984). A central concept here is the 



Loschmidt echo, also known as fidelity. It is defined as 
the overlap between two quantum states: the first state 
is obtained from an initial state |<3>) in the course of its 
evolution up to time t under a Hamiltonian H; and the 
second state results from the same initial state by evolv- 
ing up to the same time with a perturbed Hamiltonian 
H' . The measure of this overlap is 

M(t) = \(^\e lH, ^ h e- lHt ^\^}\ 2 , 

where h is Planck's constant. This quantity can also be 
interpreted as the overlap of the initial state and the state 
obtained by first propagating this state up to time t with 
H, and then backward in time, up to —t under H' . M(t) 
equals unity at t = 0, and typically decays in time. 



The term Loschmidt echo refers to the debates in the 
XIX. century about the ireversibility of thermodynami- 
cal systems and the foundations of statistical mechanics 
(in which Loschmidt participated). Similar echo con cepts 
are important in classical systems (Eckhardt, 2003), with 



applications in sensing techniques (Taddese et al., 2010), 



and have recently been extended to relativistic quantum 
systems (Sad urm and Seligman| |2008[ ). The term fi- 
delity is used mainly in the field of quantum information. 
The overlap M{t) is known to quantify the robustness of 
systems to perturbations, and has nowadays numerous 
applications in quantum information, statistical physics, 



and quantum chaos (Gorin et al. 2006). 



The first and most natural investigations in quantum 
chaos considered glo bal perturbations (e.g., H' = H + sV 
with V acting globally). Global perturbations affect all 
(or a dominant part of) the phase-space accessible to the 
quantum particle in the course of its time evolution. It 
was shown that the decay of M(t) with t has a variety of 
regimes, many of them reflecting classical properties of 
the closed system (e.g., Lyapunov exponents). See (Gorin 



et al. 



2006; Jacquod and Petitjean 



for 



. 2009D 
2011) for limitations 



and ( |Garcfa-Mata and Wisniacki| 
of the Lyapunov regime. 

More recently, experimental realizations in bil- 



liards ( jHohmann et al. , 2008) emphasized the impor- 
tance of the case of quantum mechanically local pertur- 
bations (Ares and Wisniacki[ |2009l |Goussev and Richter 



2007} |Goussev et al.\ |2008[ ]kober et al.\ |2011| ). Be- 



low we review recent results on local perturbations that 
show that in this case the decay of M(t) has an impor- 
tant regime dominated by properties of the corresponding 
classical system with leak (e.g., the escape rate). 



Figure [32] illustrates the concept of local perturbations 
on the example of a chaotic billiard of area A. The per- 
turbation consists of a deformation of width w localized 
in a region Bi of the boundary. The rest of the boundary, 
£?o, is unaffected by the perturbation. Three trajectories 
starting from the same location i*o (and reaching point 
r afterwards) are shown. The two nearby trajectories 
Si and Si marked by continuous and dashed lines corre- 
spond, respectively, to an evolution with the unperturbed 
and perturbed Hamiltonian. In a semiclassical approxi- 
mation, the action difference between these trajectories 
should be determined. The third trajectory (So) hits the 
boundary only at Bq. For this case the perturbed and 
unperturbed trajectories coincide, hence the action dif- 
ference of this pair is zero. 

The quantum mechanical implementation and inter- 
pretation of the Loschmidt echo with a localized pertur- 



bation is illustrated by Fig. [33] The evolution of a wave 
packet is shown in the original billiard up to time t. At 
this time the perturbation along the boundary is intro- 
duced, and the dynamics is followed in the time-reversed 
evolution up to the same time. Gray shading represents 
the probability distribution. The initial Gaussian wave 
packet appears thus as a dot in the upper left panel. The 
arrow in this panel marks the momentum direction of the 
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FIG. 32 Schematic diagram of a billiard with a local pertur- 
bation B\ along its bounda ry and three c lassical trajectories 
between 1*0 and r, from (Go ussev and Richter[ |2007| . 



wave packet. The lower left panel shows that the original 
form is only partially recovered after the full process, the 
overlap is thus incomplete M(t) < 1. 




FIG. 33 Forward-time wave packet evolution in an unper- 
turbed diamond billiard up to a time t corresponding to about 
lO(tcoii), followed by the time-reversed evolution in the per- 
turbed billiard. The width of the localized perturbation is 
w = 60, 15 percent of the length of the horizontal edge, 
and the de Broglie wavelength is A = 4.8 in the same units. 
From (iGoussev and Richterl 12007). 



Goussev and Richter considered a classically weak but 
quantum mechanically strong perturbatio n in a re gion 
of the boundary of a fully chaotic billiard (Goussev and 



Richter, 2007). More precisely, the length w and depth r 
of the perturbed boundary should be much smaller than 
the length of the billiard boundary, and the perturbation 
is considered quantum mechanically strong if w and r 
are much larger than the de Broglie wavelength A. Their 
main result is that in such cases there is a large time 
interval for which the Loschmidt echo decays as: 



M(t) ~ exp(-2«t), 



(99) 



where n is the classical escape rate of the billiard with a 
leak / extending precisely through the perturbation re- 
gion of length w along the original boundary. The time 
interval over which decay (99) is valid extends from the 
mean collision time (flight time) (t co ii) to a saturation 



time t s -l/(2«)lnMoo where = X(j/{2ttA), with 
a as the wa ve packet size at t = 0, and A as the area of 
the billiard ( Gutierrez and Goussev| 2009). Note that t s 
can be even larger than the Heisenberg time tn ~ l/h. 
This result is obtained by showing that - in this regime 
- the main contribution to M(t) is due to trajectories 
that do not collide with the perturbation region (colli- 
sions there would radically change the dynamics of the 
trajectories). In our formulation, this means that such 
trajectories start form the stable manifold of the saddle 
characterizing the billiard with this leak. Fig. [34] shows 
results of numerical simulations in the billiard depicted in 
Fig. 33 The quantity M(t) is plotted for different depths 
r, and its decay is governed by the classical escape rate 
k of the billiard with a leak of length w, independent of 
the value of r. Time is normalized by the classical mean 
collision time (£ C oii)- 
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FIG. 34 The decay of M(t) in the diamond billiard of Fig (33 
for four different values of the depth r at w = 60. Time is 
scaled with (t co ii) = £/. The solid straight line corresponds 
to exp (— 2 k£) with the escape rate given by Q. The inset 
presents the decay of M(t) , the Loschmidt echo a veraged over 
several values of r, from ( |Goussev and Richter[ [2007 ) . 



These results have been refined and generalized to per- 
turbations which are quantum mechanically not necessar- 
ily strong. The decay of M(t) remains then exponential, 
but the decay rate is a factor a times the classical escape 
rate: 

M(t) ~ exp(-otKt). (100) 
For piston-like perturbations, semiclassical expressions 



have been derived for th e coefficient a (Goussev et al 
2008 Kober et al.\ |2011[ ) . As a function of the perturba- 



tion strength, a turns out to be an oscillatory function 
converging to 2 for large perturbations. 

These results have been confirmed in recent experi- 
ments in m icrowave cavit ies performed by Kober and 
co-workers (Kober et al. 2011). The correspondence 
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between the Helmholtz equation (describing 2D mi- 
crowave cavities) and the Schrodinger equation make mi- 
crowave cavities a unique tool for the experimental in- 
vestigation of quantum chaos ( Richter 1999 Stockmann| 
|1999[ ), including the exploration of the effec t of classical 
trajectories in chaotic systems with leaks (Dembowski 



in Sec. II.B (Casati et al. 1999a). Open quantum sys- 



et al. 



Ref. (Kober et al. 



2QQ4|) . Figur e |35| shows the configuration used in 
2Q11[) for the experimental investi- 



gation of Loschmidt echo decays. For not very strong 
perturbations, very good agreement has been found with 



the semiclassical coefficient a in Eq. ( 100). On the exper- 
imental side, the investigation of the strong perturbation 
case remains a challenge since the signals are then ex- 
pected to be rather tiny. 

Altogether these results on Loschmidt echoes show 
that there is an increasing recent interest in systems with 
leaks, not only in applications but also on the side of fun- 
damental properties of quantum chaos. In the semiclas- 
sical limit such systems reflect the properties of the un- 
derlying classical systems, including semi- classical correc- 



tions to the classical escape rate (Gutierrez et al. 2009 



Sieber and Richter, 2001 Waltner et al, 2008). 




FIG. 35 Geometry of the chaotic billiard used in the exper- 
iments of Kober et al. The local perturbation is piston-like, 
of width w and of depth h. The measuring antenna is placed 
at location a. The rot at able ellipsis is used to perform en- 
semble averages, and the additional elements are introduced 
to red uce the influence of bouncing balls. From ( |Kober et al.\ 
^OlT). 



2. Fractal distribution of eigenstates 

Another signature of chaos in open quantum systems 
is provided by the fractality of certain invariant sets of 
the corresponding classical systems. The eigenstates of 
the quantum system turn out to be distributed along 
the manifolds of the classical chaotic saddle 14 discussed 



14 In spite of the fact that repellers (i.e., unstable sets with only 
expanding directions) cannot exist at all in Hamiltonian systems, 
this misleading term is very often used also in the quantum liter- 
ature to refer to the chaotic sad dle |Eberspacher et a/T]|2010||Lu 
et al.||2003b||Nov aes et al. ||2009||Pedrosa et al.||2009||Ramilowski 



terns are characterized by non-unitary evolution opera- 
tors having a set of right and left decaying nonorthogonal 
eigenfunctions. The left and right eigenstates of the non- 
unitary propagator concentrate in the limit of h — >• on 
the stable and unstable manifold of the chaotic saddle, 



respectively (Keating et al. 



tistical properties of energy 



2006). This affects the star 



evels discussed below. 



1979). 



Classical maps can be quantized ( Berry et al. 
Quantized area preserving baker maps with leaks (of- 
ten called open baker maps in the quantum chaos com- 
munity) nicely exemplify the results mentioned above, 
and can b e d irectly related to the findings reported in 
Sec. | III.E| ( Keating et al , 2006) studied the ternary 



baker map ( 80 ) with a leak taken as the middle horizon- 
tal strip, exactly as in Sec jIII.E Figure 36 shows their 
results: the long-lived (right) eigenstates concentrate in 
the semiclassical limit on the unstable manifold (horizon- 
tal filaments parallel to the q axis). Correspondingly, the 
momentum representation of any long-lived right eigen- 
state is supported on a Cantor set which is the projection 
of the unstable manifold on the momentum (p) axis. This 
is apparent in Fig. [37| which also shows a magnification 
that reveals the self-similar character of this distribution. 
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FIG. 36 Representations of the long-lived right eigenstates 
in the ternary baker map (with a leak). The classical map 
is described by (80) but the variables are here p = x and 
q = y. Here only 1/9 of the phase space is shown, outside 
the leak (a central vertical strip of area 1/3). Left panel: 
The average of the so-called Husimi function. Right panel: 
the corresponding Wigner function average. The scale of h is 
0.005, from (iKeating et al. 112006). 



In another study, Pedrosa et al. considered the quan- 
tized version of the binary baker map (82) with leaks 



as s tripes along the p ax is, centered at qi and of width 
Aq ( Pedrosa et al.\ 2009). In the classical version, they 
found a similar strong dependence of the escape rate on 
the leak position q\ (for any fixed finite Ag), as discussed 
in Sec. IV. A The influence of the shortest periodic orbit 
falling into the leak has also been noticed. Looking for 
a quantum mechanical analog of this, they determined 
the distribution of eigenvalues of the quantum evolu- 
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et al. 



[2009]). 



tion operator in the complex plane, as shown in Fig. 
Since the map is open, moduli are less than one, and the 
eigenvalues fall thus all inside the unit circle. For a leak 
with smaller escape rate, the distribution of eigenvalues is 
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FIG. 37 Probability density of the leaky ternary baker map 
in momentum space from right eigenstates averaged over the 
longest-lived states. The magnification (right pan el) illus- 
trates the s elf-similar Cantor set character, from ( |Keating 
et al. \ 12006). 



found to be dense at the outer ring, while the distribution 
is characterized by an increase of density near the origin 
for a case with faster escape. The strong dependence of 
the classical escape rate on the position of the leak is thus 
found to be reflected in the eigenvalue distribution of the 
quantum evolution operator. 
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tic maps fiw = d/2). In open chaotic maps the num- 
ber N(k) of resonances with wave nu mbers of real par t 
smaller than k scales with an exponent ( Lu et al. 2003a ) : 



(102) 



where Do is the fractal dimension of the underlying clas- 
sical chaotic saddle, and = is the partial box- 



counting dimension (see (29). Equation (102) shows that 



the original Weyl law is converted into a fractal Weyl 



law (Eberspacher et al. 


2010 Ermann and Shepelyan- 


sky[ |2010l Kopp and Sc 


iomerus[ 2010| |Lu et al.\ 


2003a 


Nonnenmacher[|2011||Nonnenmacher and Zworski 
Ramilowski et a/.||2009||Shepelyansky||2008; Spina 


2005 

, et al. , 


2010 Wiersig and Main, 2008). An investigation of 



this p roperty has also been carried out (Pedrosa et al. 
2009). As Fig. 39 shows, they determined the fraction 
of eigenvalues N v with modulus v > 0.3 as a function 
of the Hilbert's space dimension, N, for different leaks. 
The prediction of the fractal Weyl law for this case is 
log N v ~ log N, which is fulfilled for all graphs. 




FIG. 39 Logarithmic plot of the fraction of eigenstates N u for 
v > 0.3 as a function of the dimension N of the Hilbert space 
in a quantum baker map. Lines correspond to the prediction 
of the fractal Weyl law. In the upper (lower) panel qi = 0.3 
(qi = 0.5), and circles correspond to Aq = 0.05, squares to 
Aq = 0.1, and triangles to Aq = 0.2 (Ap = 0.5 in all cases). 
From (IPedrosa et al. 1120091). 



FIG. 38 Eigenvalues of a leaky baker map in the complex 
plane. In the upper panels the dimension of the Hilbert space 
is N = 602, while in the lower ones N = 2048. In the left 
(right) column qi = 0.3 (qi = 0.5) with classical escape rate 
7 0.Q90 (7 0.165). Aq 0.1, Ap 0.5 in all cases, 



from (Pedrosa et al. 2009) 



Weyl's law states that in closed systems the number 
N(k) of energy levels with wave number smaller than k 
grows as 



N(k) 



(101) 



where the exponent /3w is an integer proportional to the 
dimension of the system (e.g., for <i-dimensional symplec- 



3. Survival Probability and quantum Poincare recurrences 



In Sec. II. C we have related the escape of trajectories 
in systems with leaks to the problem of Poincare recur- 
rences. In the same spirit, Casati, Maspero, and Shep- 
elianski used the quantum survival probability P q (t) to 
investigate what they define as Quantum Poincare re- 



currences ( |Casati et al. , [l999b| ) . This is performed in- 
troducing absorbing boundary conditions (the leak) in 
quantized area-preserving maps in the weakly chaotic 
regime. The quantum decay P q (t) is identical to the clas- 
sical decay P q (t) = P(t) up to a time t q ~ y/l/ti. For 
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times longer than the Heisenberg time t > tn ~ one 
finds P q (t) ~ 1/t. This is valid up to a maximum time 
tmax ~ exp(l/ft) after which an exponential decay sets 
in, whose rate is, however, different from the classical 
escape rate. 

The origin of the n ew quantum decay regim es can be 
understood following (Wimberger et al. 2002) (see also 
references therein) and writing 



(103) 



where Tj > is the decay rate of the eigenstates \ifjj) 
and uoj is the expansion coefficient uoj = K^ol^j)! 2 f° r 
the initial state \<po). The asymptotic exponential de- 
cay corresponds to the smallest Tj with ujj ^ 0. In a 
mixed phase space this longest living state is localized 
in the center of a KAM island and has cuj =^ due to 
chaotic tunneling. KAM islands (or sticky regions) are 
not the only source of power-law decays in P q (T), which 
can be originated through genuine quantum effects such 
as the the quantum mechan ical (or dynamical) localiza- 
tion of states due to chaos ( |Stockmann| [1999 ) . In par- 
ticular, even classical systems showing exponential tails 



in P(t) have power-law tails in P q (t) (|Alt et al. 


1996 


1995 


Casati et al. 


1999b 


Fendrik and WisniackiJ | 


1997J). 



The full description of P q (t) has to consider the effect of 
the localization of many states, not only in the chaotic 
and regular components of the phase space but also at 
the hierarchical border of the KAM islands (Ketzmerick 



et al. , 2000). In general, P q de pend s both on the distribu- 
tion of the decay rates T's in ( |103 ) and also on the initial 



condition 



as seen from Eq. (103) and emphasized 



in (Wimberger et al> 2002). The existence of regimes 



of power-law decay is generic, but exponents different 
from Pgjt) ~ l/t ar e also typically found ( |Skipetrov and 
|van Tiggelen , 2006). Altogether, the decay of P q {t) and 
the previously mentioned Loschmidt echo (fidelity) decay 
show that the classical decay is still observed in quan- 
tum systems, but quantum effects dominate the large 
and asymptotic behavior. These deviations are typically 
not accessible for experiments [for more discussions on 
the time when P q (t) deviates from P (t) see (Puhlmann 
et a/.||2005| |Savin and Sokolov[ [l997| )]. 



The most detailed experiments of the survival proba- 
bility are on atom-optics billiards as r eported in ( |Fried- 



man et al. 2001 Kaplan et al. 2001). In these exper- 



iments the survival probability of atoms inside cavities 
were measured for different boundaries and leak posi- 
tions. The connection to classical billiards becomes evi- 
dent when comparing the escape from integrable (circu- 
lar) and fully chaotic billiards (tilted Stadium billiard 15 ). 
The measurements shown in the upper panel of Fig. |40| 



15 The usual parallel walls of the Stadium billiard were tilted in 
order to avoid the influence of the bouncing ball orbits discussed 
in Sec.lvTB] 
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FIG. 40 Survival probability of ultracold rubidium atoms 
confined by an optical potential. Upper panel: diamonds 
(circles) correspond to the experimental results for a circu- 
lar (Stadium) billiard with a leak at the bottom, as shown 
in the upper (lower) insets. Lower panels: different poten- 
tials along the boundary of the stadium billiard, with two 
different positions (a) and (b) of the leak. In both panels re- 
sults with harder (o) and softer (+) billiard walls are shown. 
Full and doted lines correspond to numerical simulations and 
the dashed line indicates Sabine prediction k — / (t C oii) , 



Eq. ( 


13) 


2001 





From ( Friedman et al.\ |2001| and ( Kaplan et al. 



confirm that a much faster (exponential-like) decay is 
observed for the chaotic case than for the i ntegrable 
one. A subsequent study QKaplan et al. , 2001[ ) consid- 
ered the effect of soft walls along the boundaries of the 
(tilted) Stadium billiard, which generically lead to the 
creation of KAM islands. The lower panels of Fig. |40| 
show the experimental survival probabilities. They con- 
firm a slower decay for soft walls, which was shown to 
be related to t he stic kiness of trajectories around the is- 
lands (see Sec. V.B). A strong dependence on the leak 



position is also found. If the leak corresponds to a region 
of the phase space that fully contains the island(s), no 
stickiness effect can be observed (panel b). 



VII. SUMMARY AND OUTLOOK 

There are many different configurations in which a leak 
is introduced in a chaotic systems, but in all cases the 
simplest theoretical approaches rely on the same princi- 
ple: the properties of the leaky system can be compared 
to the properties of the closed system. Simple estimates 
based on this principle led to Sabine's law for acoustical 
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reverberation at the end of the XIX century and, more 
than a century later, were used to expl ain the emission 
patterns of microlasers (see, e.g., Ref. ( Lebental et al.\ 
2QQ7|)) and the survival probability of c old atoms ( |Fried-| 
man et al t 2001 Kaplan et all 12001 1). In this paper 



we have presented a transient chaos based theory of this 
problem that goes beyond the (naive) closed system ap- 
proximation (Sec. II. A ) and leads to a proper treat ment 
of finite leaks (Sec. |ILB| ), partial reflection (Sees. II. D 
and III), and weak chaos (SecM). 



Another aim of this paper has been to illustrate the 
abu nda nce of applications of leaking systems (Sees. I] 
and VI). In the regimes in which dynamical-system mod- 



els apply, the results reviewed show that all important 
quanti ties of the t ransi ent chaos theory of leaky systems 
(Sees. II.B, III.D, and IV). have direct consequences in 



experiments and applications of contemporary physics re- 
search: 



• the escape rate (Sec. III.C| ) dominates some 
regimes of the Loschmidt echo (Sec. VI.E.l) and 
was directly measured in cold-atoms experiments 
(Sec. IVIE.3D; 



• the chaotic saddle and its invariant mani- 



folds (Sec. II.B) leave di rect fi ngerprints on as- 
tronomical problems (Sec. VI. A), hydrodynamical 
flows (Sec. VLB), magnetic field lines i n toka maks 
(Sec. |VI.C| ), optical microcavities (Sec. VI.D), and 
on the dist ribution of eigenstates in quantum sys- 
tems (Sec. VI.E.2). Many of these examples in- 



clude weakly chaotic systems, where the observa- 
tions often become easy to interpret in view of the 
division of the chaotic saddle in hyperbolic and 
non-hyperbolic components (Sec. V.D). 



• the stickine ss du e to the no n- hyperbolic com- 
ponent (Sec. V.B) becom es observable in trans- 
port properti es (Sec. V.D) and in quantum systems 
(Sec.|VLE.3). 



the dependence on leak position, size, and 



shape (Sec. IV. A| ) is clea rly ob servable in hydro- 
dynamical problems (Sec. VLB), in optical micro- 
cavit ies (Sec . VI.D), and in quantum experiments 
(Sec.|VI.E.3|). 



• basins of escape (Sec. IV.B ), shaped by an under- 



lying chaotic saddle, prove to be concepts usefull y 
applicable to the three body pr oblem (Sec. VI.A), 
in asymmetric trans port ( Sec. V.D), and in the 
plasma problem (Sec. |VI.C ); 



• far-field emission (Sec. |IV.C[ ) and partial reflec- 
tion (Sec. II. D) play a crucial role in experiments 
in optical mic rocav ities (Sec. VI.D), and in room 
acoustics (Sec. LB). 



Our review paper contains also new results. They 
mainly arise within the framework of our operator de- 
scription of true-time maps with partial leaks (Sec. III). 



The general formula (|71| relate averages taken with the c- 
measure to the escape rate (e~ Ktcoll )c = (R)c- The idea of 
a partial reflection/absorption is not incorporated in tra- 
ditional descriptions of dynamical systems, despite being 
physically very natural, and also of easy implementation 
in ray simulations. Here we have shown that reflectivity 
can be naturally incorporated into the operator formal- 
ism that avoids the artificial factorization of the total 
energy density p in ray density p and intensity J. More 
generally, the approach makes superfluous the attaching 
of labels (e.g., intensity or true time) to trajectories, that 
are modifie d at eac h collision. Also the numerical proce- 
dures (Sec. VIII. B ) to estimate the different measures of 
leaky billiards should be relevant not only for simulations 
but also in experiments for which data can be collected 
within a snapshot (as in the S set) or during some time 
interval (as in the Q set). The significance of these re- 
sults is not restricted to billiards or even to Hamiltonian 
systems, they apply to any true-time map obtained from 
a Poincare surface of section. The collision times is then 
to be replaced by the return times to the Poincare sec- 
tion. It is a widespread practice to assume that results 
for maps extend immediately to flows or true-time maps. 
While this holds for stroboscopic maps of periodically 
driven systems, autonomous problems, or more gener- 
ally, problems with Poincare maps are different. Our 



results further emphasize the message of ( |Kaufmann and 
Lustfeld 2001) that the connection between true-time 



and discrete-time maps in open systems is different from 
the case of closed systems, and often involves surprising 
non-trivial results. 

Finally, we discuss unsolved problems and future re- 
search directions. First we should note that leaky sys- 
tems are intimately related and provide further motiva- 
tion to fundamental problems in dynamical systems. In 
particular, in the case of Hamiltonian systems we can 
mention the understanding of the mechanisms of sticki- 
ness in higher dimensions and its connection to Arnold 
diffusion (Bunimovich 2008) and a formal treatment 



of the division of the chaotic saddle in hyperbolic and 



non- hyperbolic components (see e.g. (Altmann and Tel 
2008)). Among questions directly related t o leaky sys- 
tems, and in addition to ( jDettmann , 2011b), the follow- 
ing general problems claim for further investigations: 

• Adapt and extend current results on sensitivity on 
leak position to more generic situations. For in- 
stance, most rigorous mathematical results about 
the escape rate (see Sec. IV. A) are demonstrated 
for one-dimensional systems with leaks in Markov 
partitions. A similarly non- generic feature is the 
standard explanation of directional emission in 2D- 
optical microcavities (see Sec. VI.D), which relies 



on the existence of an unstable periodic orbit close 
to the leak. There is an evident need for extending 
these results to generic leak positions and chaotic 
saddles, respectively. 

• Obtain the conditions for which a single chaotic 
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saddle (with a single physically relevant c-measure) 
is created after a leak is introduced (iCollet andl 



Maume-Deschamps , 2000J Demers and Young[ 
2006|. The hypothesis of the existe nce of a single 



saddle was used in our theory (Sec. II. B| ). To ex 
plore the cases in which this is violated ( Dettmann 



and Georgiou, 2011b) remains a task for the future. 



• Consider novel configurations of leakage such as 
introducing a temporal dependence or a random 
choice on the properties of the leak (Bahsoun and 
Vaienti| [20121 |Georgiou et al.\ |2012j |NandorTand 
Szaszl 120121). 



• Partial leaks: clarify the nature of the spectrum 
of dimensions of the invariant sets (as mentioned 
in Sec. II.D), related these dimensions to the es- 



cape rate and to the Lyapunov exponent. Extend 
the results to non-invertible systems, and develop 
efficient algorithms for the computation of chaos 
characteristics (e.g., escape rate, dimensions, and 
Lyapunov exponents). 

• Recent works suggest that certain non-traditional 
(e.g., location dependent) boundary conditions in 
problems of wave cha os can be inte rpreted as gen- 
eralization of leakage ( |Berry| |2009[ ) . It remains to 
be explored how this can be related to the tradi- 
tional leakage discussed here. 

Investigate further the connection to transport 



properties ( Dettmann and Georgiou 2011b Knight 


et al. 


2012j~Kuhl et al\ 2005} and the effect 




of stochastic perturbations ( 


Altmann and Endler, 




2010 


|Altmann et a/.[|2012||Dettmann and Howard, 


2009 


). 



Besides these general problems, a clear future research 
line is to adapt the models to make them more realis- 
tic to specific applications. This involves including new 
features in the dynamics, e.g., gain medium in lasing cav- 
ities, reactions of particles, or properties of the wave sys- 
tems (e.g., modified collision laws). The challenge here is 
to show how to connect experiments and observations to 
the theory proposed in this paper. Also important is to 
show in which extent and into which level of detail the 
models match the experiments, e.g., in which extent non- 
hyperbolic properties of chaotic dynamical systems play 
a role in the experimental results. In turn, these adap- 
tations towards applications bring new problems and in- 
spiration for theoreticians, a virtuous circle illustrated 
in this paper by, e.g., the case of partial leaks. We are 
convinced that the collection of problems that can be 
modeled as leaky chaotic system will keep growing and 
pushing the research forward also in unpredictable direc- 
tions. 



VIII. APPENDICES 

A. Projected measure and averages 

Let us turn now to the c-measure jip c of the dynamics 
obtained by projecting the flow on the boundary. This de- 
scribes in an ensemble of escaping trajectories the prob- 
ability of the occurrence of collisions with Birkhoff co- 
ordinates x irrespective of the time of collision. Let the 
corresponding density be denoted as pp c (x), then 

dj^pc = pp c (x)(ix. (104) 

The projected measure is obtained by integrating the flow 
measure d/npc in Eq. ( [63] ) over coordinate time up to its 
maximum, t C oii( x ) (see Fig. 10). This measure can thus 
be written as 

dfiPc = Jo co " (x) dfi Fc = / * coii(x) Ace-^drdfic = 



-Ktcoll(x) 



-A c d/i c . 



The relation between the densities is then 



p Pc (x) = p c (x) : 



-Ktcoll(x) 



-A r . 



icoii(x) 



(105) 



, (106) 



K ' (^coll)/i 

where the arrow indicates the limit of closed systems 



(k — » 0). Equivalently, 
p Fc (x,r) = pp c (x): 



, > (107) 

1 _ e -«*coii(x) t coll (x) v ' 

Using these relations, we can write average values in 
the projected and flow measures as a function of those 
in the true- time map measure. Particularly interesting is 
the case of the average collision time. In the flow measure 
it is obtained as the average of the coordinate time r 

/Mo 00,l(x) 



(^coll) 



Fc 



dr re Kr A c p c (x.) 



I _ e -Kt co n(x) _^ coll ( x )e- Ktcoll ( x ) A , N 

J dx ^cpc(x) 

l-(e- KU )o-(«tcoiie-^")c, , <i c 2 oll ) 



-A a -> 



k- 2(t co n) 
(108) 

In the projected measure the average collision time is 
obtained from ( |106| ) as 

(t co \\}pc = Jdxt co ii(x)/9p c (x) 

= A c J dx p c (x) 

(^coll)c - (£colie"^ coll )c . 



(^coll) _ ey/j. \ 
— r - Z{t co n) F . 

\^coll/ 



(109) 
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The difference between the average collision time in the 
projection and in the true-time map is then 



(tcoll) Pc — (^coll)c = 

((icoii)c(e- K *->)c - (tcoiie- Kf »"( x )) c ) — -»• 

« (*coll) 



(110) 

where of coll is the variance squared of the distribution of 
the collisions times in the closed map. The difference be- 
tween the two averages is thus due to the inhomogeneity 
of the collision time distribution. 

An even simpler relation can be obtained between the 
averages taken in the true-time map and the projected 
representation. Rearranging Eq. ( 106 ) and using the fact 
that p c is normalized, we obtain 

J p c (x)dx = 1 = J p Pc (x) ~ _ e _^ tcou(x) rfx, (111) 



from which it follows that 
1 



( e -«£coll) c 



-K*coll 



Pc 



For closed systems this relation goes over into 



1 



(^coll) 



1 

^coll / p c 



(112) 



(113) 



B. Algorithms for open billiards 

Here we describe the numerical simulations of open bil- 
liards. For the limacon billiard considered in this paper 
a map was obtained using the standard geometric tech- 
niques (Back er and Dullirij |1997 Robnik, 1983): given a 
collision position r = (x,y) and velocity ^"[correspond- 
ingly, coordinates the next collision position is 
obtained as the first intersection of the ray r + vt with 
the billiard boundary. That is achieved by imposing that 
the point r + vt satisfies Eq. p0| ) of the cardioid, which 
can be written for S = 1 as (x 2 -\-y 2 — ex) 2 — (x 2 -\-y 2 ) = 0. 
The smallest solution of the resulting forth order polyno- 
mial in t yields the collision time, £ C oii( r/ )> and the next 
collision position, r' = r + v£ co ii. 

We search for numerically efficient procedures to com- 
pute observable quantities in true-time open maps. In 



Sec. |II.B| we learned that asymptotic quantities and in- 
variant properties of the escape process are given by the 
density p c ( x ) of the c- measure. Below we describe how 
this density can be approximated in numerical simula- 
tions, related to the different representations of invariant 
measures. 

For the billiard flow, p c is approximated by taking all 
trajectories in an initial ensemble that survive up to a 
large time t*. Mathematically, t* — » oo, but in practice, 
for t* > 2/ n a good approximation of p c is obtained. For 
true-time maps, we have only the Birkhoff coordinates 
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FIG. 41 Schematic illustration of the numerical procedure to 
obtain the true-time c-measure in open billiards. The time 
instants of collisions at the boundary are plotted, with black 
dots, for different trajectories. Trajectories i = 2 and i — 3 
belong to S(t* — 35) (inset) and have n*(i = 2) = 17 and 
n*(i = 3) = 13. S(t*) denotes the set of trajectories that 
survive inside the billiard up to time t* at least. Trajectories 
i = 1, 4, N do not belong to S(t* = 35). 



x n = ( s mPn) and the time t n of collisions. It is natu- 
ral to take for each trajectory surviving up to time t* 
the values s n * and p n * , where n* is the smallest n such 
that t n * > t*. Figure 41 illustrates this procedure. Note 



that n* is different for each trajectory and depends on the 
whole history of each trajectory [through Eq. (19)]. The 
points x* = (s n *,p n *) obtained through this procedure 
group together all points that at t = t* have coordinates 
x* irrespective of their coordinate time r after the last 
collision (small dots in Fig. [To| Therefore, these points 
are then sampled according to pp c (x), the projected den- 
sity 16 . This means that the average (B)p c of any ob- 
servable B{pc) (e.g., -B(x) = t co ii (x)) can be numerically 
obtained as 



(B) 



Pc 



N- 



where S(t*) denotes the set of trajectories that survive 
inside the billiard up to time t*. Qualitatively speaking, 
the use of the x* coordinates belonging to set S(t*) cor- 
responds to an instantaneous observation of the system. 
In a simulation using N initial conditions, let the size of 
S(t*) be denoted by N* . The survival probability can 
then be estimated as P(t *) = N* /N. Furthermore, it is 
possible to use Eq. ( |1Q6 ) and this numerical procedure 
to efficiently compute the true-time map average value of 
any observable B(x) as 



(B) c = J dx p c (x)E(x) 



16 Equation jl06| ) shows that pp c (x) is insensitive to the sliding 
trajectories with vanishing £ C oll( x )- 
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/ 



= / dx 



1 - (e~ KtcoU ) c 



— t - 

TV* ^ 1 

S(t*) 



-/Ctcoll(x) 



l-(e- 



Pp c (x)B(x) 



:S(x*), 



l(x* 



(114) 



where the last relation is valid for 7V*,t* ^> 1. In partic- 
ular, with the c hoice B(jc) = 1, one recovers Eq. (112). 



Equation (109) is recovered with the choice B(x) = 
*coii(x)(l - e~ ^ coll(x) ) /^. Once the value of (e-^ coll ) c 
is known, Eq. ( |VIII.B ) can be directly applied to an ar- 
bitrary observable. The measure of the leak p c (I) is com- 
puted, e.g., by choosing B — 1 for (x) G / and B = for 

(x)^I. 

For partial leaks discussed in Sec. |II.D Eq. (VIII.B) is 
generalized to 



(B)c 



r* S Jt 'h 

S(f) 



-Kt co ll(x*) 



^(x*), (H5) 



where Jt* is the intensity of the trajectory over which 
summation runs, and J* = J t * is the total intensity 
of the rays in S(t*) (i.e., of all trajectories in the system 
at time t*), which are all trajectories if i?(x) ^ 0. The 
advantage of this procedure, and of any application us- 
ing the S(t*) set, or of projected densities, is that we can 
have good control on the time t = t* at which the mea- 
sures are approximated numerically. This is particularly 
important for systems with mixed phase space discussed 
in Sec. [V| where a crossover from exponential to power- 
law decay is observed in P(t). 

An alternative procedure for generating averages in the 
true-time map is based on sampling a set of collision 
points according to the density p c of the map. This is 
done by collecting in a set Q(t*,At*) all collisions in a 
time interval At* = [£*,£* + At*] for At* > max{t coll (x)} 
and over all N trajectories. The use of the Q set assumes 
an observation of the system over an extended time pe- 
riod. 

To see how this procedure generates the c- measure, 
consider first q different time instants t* G At 
1, ...,(/, and let S q (t*) = {S(t*)} q j=v 
variant in time, points in all £(£!•) 's and in S q (t*) are 
distributed according to pp c . However, the coordinates 
x* of a single trajectory i will typically appear multiple 
times in S q (t*) because different time instants t* G At* 
may have the same n*. This can be seen by considering 
two closeby values of t* in Fig. (41). (For t* - t* +1 < 
(^coii)c? most trajectories will leaaTo repeated values of 
x*.) For q — >• oo, the number of times a given x* appears 
in S q (t*) is proportional to J o tcoll(x * } e~ Kt dt. 

A point x* of a single trajectory is included only once 
in Q(t*,At*) and possibly many times in S q (t*). There- 
fore, noting that the S set is sampled according to pp Cl 
we get that points in Q(t*,At*) are sampled according 
to 

PPc(x) 

4/o COll(x) e-^dt ~ M± ~ e-«*coii(x)) 



"3 ^ i J ~ 

Since pp c is in- 



pq( x ) = 



Ppc(x) 



where A c is the normalization constant (62). Comparing 
to Eq. (106), we obtain Pq(x) = p c (x). 



Numerically it is efficient to consider the case Q(t*) = 
Q(t*,At* oo), i.e, to consider all collisions until 
esca ping the system. Using this sampling procedure, 
Eq. ( | VIII.B ) is re- written as 

Wc4 E 5 ( x *)> ( 116 ) 



Nc 



where Nq ^> 1 is the nu mber of points in the set Q(t*). 
The case of partial leak (|115|) is in this case simply 



(B)c 



1 



Jq,t 



E 

Q(t*) 



J t .B(x*), 



(117) 



where Jqj* is the sum of the intensities of all collisions 
in the set Q. In our simulations we have used the set Q 



for calculating p c (I) and (t C oii)c? as illustrated in Fig. 42 
While the set Q provides better statistics and should typ- 
ically be used, the results for the set S are conceptually 
interesting and may find applications in different simula- 
tions or experimental applications. 
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FIG. 42 Convergence of the numerical estimation of j i c {I) 
and (tcoii)c These results were obtained using Eq. (116) 
for the Cardioid billiard [e = 1 in Eq. (20)] with a leak in 
(s,p) G [0.4,0.6] x [—1, 1] using 2.5 10 7 initial conditions uni- 
formly distributed in the full phase space [—1,1] x [—1,1]. 
Short time fluctuations are present during the convergence to 
the c-measure. The results /x c (J) = 0.11745 ± 0.00006 and 
(tcoii)c = 1.9161 =b 0.0004 are the mean and the standard 
deviation of the values obtained for integers t* in the inter- 
val t* e [1/k*,4/k*] with k* = 0.05693 obtained by ([24 ) 
using fj,(I) = 0.10 and (t co ii) = 1.85055. The Q sets at 
each time t* include all collisions for t* < t < t max with 
tmax = 7 / K* ~ 125. The actual decay rate k was esti- 
mated through a direct fit of P(t) (not shown) and resulted 
in k = 0.065590 ± 0.000007. The procedures described here 
were used in all computations throughout the paper. 



C. Computation of invariant manifolds and densities 

Numerical approximations of the invariant sets of tran- 
siently chaotic maps can be obtained using the sprinkler 
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method (|Lai and Tell 120111 |Tel and Gruizl |2006|). The 



idea is to use initial conditions distributed uniformly in 
the phase space of the closed system and keep track of the 
trajectories that never escape up to a time t*. Coherent 
with the notation of the previous section, t* corresponds 
to the transient time needed for p(x, t) to approach p c (x). 
For t* ^> the surviving trajectories necessarily start 
close to the stable manifold of the saddle (at t = 0), ap- 
proach the saddle (at t « t*/2), and most of them will be 
about to leave the system through the unstable manifold 
of the saddle (at t = t*). The invariant sets are approx- 
imated by the position of these surviving trajectories at 
the times t s = {0,t*/2,t*}. More precisely, the stable 
manifold, the chaotic saddle, and the unstable manifolds 
build the support of the densities p(x, t s ) obtained from 
the position of the trajectories x(t s ) at the different times 
t s = {0,t*/2,t*}, respectively. For t s = t* we recover the 
c-density p c (x). 

Modifications of the methods above are needed to ad- 
dress true-time maps with partial leaks. In order to ad- 
dress the true-time aspect of the map, the modifications 
discussed in the previous section are needed. In the sim- 
ulations shown in this paper we have considered an ap- 
proach based on the Q-set of the previous section, with 
a fixed t* ^> 1/n and a finite At* > max{t co ii(x)}. From 
trajectories surviving up to t*, the positions x of all col- 
lisions (until escape) in the time interval t G [t s ,t a + At*] 
with t s = {0,t*/2,t*} generate the densities p(x, t s ). 
This procedure ensures that the approximation criteria 
of the method (e.g., surviving for a time t > t*, for 
the stable manifold) is satisfied by all considered points. 
It is worth noting that when initial distributions reach 
the saddle relatively fast within the chosen interval At*, 
the initial positions are a better approximation to the 
stable manifold than the positions in the time interval 
t G [0, At*]. For the case of partial leaks the same pro- 
cedure is employed, and the position x of each collision 
with t G [t s ,t 3 + At*] is counted with a weight given by 
the intensity J(t* + 1 — t s ) (the same intensity at a time 
t > t* is used for the computation of the densities along 
the saddle, stable, and unstable manifolds). Notice that 
for the case of full leak (R — 0) and fixed collision time 
(t co n = 1) the usual definitions are practically recovered. 

The general procedure described above was employed 
to compute Figs. [6] and [9] For the case of Figs. [22] and [44] 
the system shows asymptotic non-exponential decay of 
P(t) and the densities around some finite time t* were 
estimated using a small region At*. 



D. Difference between Poincare and true-time maps 

The distinction between the usual Poincare map and 
the true-time map is crucial in all open systems because 
the c-measure (or the m easure on the chaotic saddle) in 
both maps are different (Kaufmann and Lustfeld, 2001[ ). 
This difference becomes dramatic in the case of billiards 
with concave borders such as the limacon billiard consid- 



ered here, because of the existence of trajectories with 
^coii(x) = sliding along the boundaries (whispering 
gallery) characterized by p « ±1. In the Poincare map, 
these orbits build one-parameter families (the boundary) 
of nonhyperbolic trajectories. However, these orbits have 
little influence on the flow, where time is counted by t and 
not by the number n of collisions. 




1000 



FIG. 43 Difference between the true-time and the Poincare 
map of the cardioid billiard. Survival probability P(t) as a 
function of the physical time t (black) and the discrete-time 
survival probability P(n) given in terms of iteration num- 
ber n of the Poincare map (red, multiplied by the average 
collision time (t co u)c for comparison, n — n(t co n)c)- The 
leak is as in Fig. [5] and Tab. |ll] 10 8 initial conditions have 
been taken, uniformly distributed in [—0.9,0.9] x [—0.9,0.9], 
to avoiding sliding initial conditions. The log-log inset high- 
lights the power-law tail of the Poincare map. 

The flow of the closed cardioid billiard is proved to 
be hyperbolic, and we expect a well-defined c-measure in 
the leaky case both in the flow and in the true-time map. 
Figure [43] shows that the usual Poincare map is not fully 
(uniformly) hyperbolic due to the sliding trajectories. 
They convert the long time decay to non-exponential in 
the leaky case. The c-measure of the Poincare map is thus 
ill-defined. From this example it is clear that the simple 
correspondence n(t co \\ c ) \-> t between the physical time t 
and the number of iterations n of the map is not able to 
explain the qualitative difference between the (physically 
meaningful) survival probability P(t) and the discrete- 
time survival probability P(n), defined as the fraction of 
initial conditions that survive inside the system up to n 
collisions at least. From Fig. [43] one can have the impres- 
sion that the differences appear only for long times. In 
Fig. [44] we present, however, p c of the Poincare map nu- 
merically obtained as usual, which should be compared 
to the true-time map result shown in Fig. [6] There is 
a disagreement not only in the region of the sliding or- 
bits. Indeed, these orbits represent just an extreme case 
of discrepancy. 

The c-measure of the Poincare map and of the true- 
time map are generically different (even in the absence 
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of sliding orbits) k ^ 7p m ap/(^coii)c 7 and there is no sim- 
ple relationship between k of the flow and 7p ma p- For the 
cardioid billiard 7p m ap is n °t even defined, as P(n) de- 
cays asymptotically as a power law. If, nevertheless, we 
consider an effective 7p ma p extracted from the short time 
behavior in Fig. 43 , we obtain that 7p ma p and K,(t co \\) c 



differ by about 2%7as can be seen from Tab. [TTJ and the 
data reported in the caption. 




FIG. 44 Conditionally invariant measure for the Poincare 
map of the cardioid billiard obtained for initial conditions uni- 
formly distributed in [-0.9,0.9] x [-0.9,0.9]. This measure 
was obtained from trajectories which survive at least up to 
n* = 10 iterations. For trajectories with life time longer than 
n* all iterations up to 40 are plotted. In comparison with the 
corresponding true-time map measure, Fig. [6] obtained with 
t* = At* = 80, orbits in the sliding-orbit region (\p\ « 1) have 
a bigger weight, although in the range 10 < n < 40, escape is 
essentially exponential (cf. Fig. 43). 
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